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ABSTRACT 


Experimental  and  analytical  investigations  on  jet  noise  generation  and 
suppression  have  resulted  in  a  verification  of  the  concept  that  jet  noiss.  ie 
uniquely  related  to  mean  flow  velocity  distribution  in  the  jet  wake.  Advances 
have  been  made  in  analytical  techniques  for  calculating  aerodynamic  properties 
of  free  jet  flow  and  the  related  acoustic  properties  of  the  flow  and  verification 
of  the  assumption  that  negligible  attenuation  of  sound  occurs  during  propagation 
through  multi -element  jet  flows  has  bean  achieved. 

Investigation  of  selective  water  injection  as  a  jet  roiae  suppression  scheme 
has  indicated  little  promise  for  in-flight  noise  reduction.  Studies  of  noise 
generation  by  high  temperature  and  high  preasure  ratio  jets  have  extended  the 
state-of-the-art  in  the  area  of  high  temperature  flow  noiae  suppression,  and  have 
resulted  in  a  means  for  predicting  total  acoustic  power  beyond  the  limits  of 
applicability  of  the  Lighthill  parameter. 
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This  Is  tbe  final  report  covering  the  work  performed  under  Phase  I,  "Jet 
Boise  Suppression  Research",  Department  of  the  Havy,  Bureau  of  Ksval  Weapons, 
Contract  Bo.  JRv  6S-C687-d,  in  accordance  vith  the  Contractor's  paropoeal 
entitled,  "Jet  Boise  Suppression  Research  and  System  Studies",  P&L-86,  dated 
March  1961. 

The  objectives  of  the  Phase  I  program  were:  to  extend  the  theoxy  of  jet 
noise  suppression)  to  high  pressure  ratio  and  high  temperature  flow  by  demon¬ 
strating  the  generality  of  the  concept  that  jet  noise  is  uniquely  related  to 
the  floe  velocity  distribution  in  the  wake;  to  further  refine  and  develop 
analytical  techniques  applicable  to  computer  solutions  far  jet  flows  and  noise; 
to  investigate  significant  variables  Influencing  Jet  noise  generation  and 
suppressor  performance;  and  to  evaluate  the  feasibility  of  selective  water 
injection  as  an  in-flight  suppression  technique.  The  advanced  knowledge 
achieved  through  the  Successful  accomplishment  of  this  research  program,  along 
with  the  findings  of  the  Phase  II  program  on  operational  noise  problems  and 
effects  of  suppression  on  mission  performance,  can  be  implemented  toward  the 
goal  of  developing  a  practical  suppressor  for  military  jet  aircraft. 

The  problems  dealt  with  in  this  research  investigation  have  been  divided 
into  the  five  sections  of  thin  report,  each  of  which  la  essentially  self- 
contained  far  the  convenience  of  the  reader.  However,  the  interrelation  of  the 
various  aspects  of  the  program  Is  maintained  by  appropriate  cross-references 
between  the  various  sections,  and  by  the  condensation  of  the  overall  effort  in 
the  Suamaxy. 
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Experimental  end  analytical  inve  s  tigat ions  on  Jot  noise  generation  nod  suppression 
ham  been  carried  out  for  the  purpose  of*  1)  d«onatratiug  and  axUnding  the 
generality  of  the  concept  that  jet  noiee  la  uniquely  related  to  jet  wake  mean 
flow  velocity  distributions;  2)  further  refining  and  developing  analytical  com¬ 
puter  techniques  for  solution  of  jet  flow  and  noise  problems;  3)  defining  the 
effects  of  jet  teaperature  and  pressure  ratio  on  noise  characteristics  for  nosxles 
and  nossls  suppressors;  4)  determining  the  effect  of  the  turbulent  jet  wake  on  the 
radiation  of  sound  energy  generated  by  Multiple  element  suppressors;  and  5) 
evaluating  the  feasibility  of  selective  water  injection  as  an  in-flight  suppression 
dsvioe. 

The  previously  established  technique  for  calculation  of  sound  power  spectra  from 
aerodynamic  properties  of  the  jet  wake  la  outlined  in  Section  1  of  this  report. 

It  is  demonstrated  that  the  method  la  quite  general,  being  applicable  to  noxale- 
suppresaors  of  various  geometric  designs.  Furthermore,  experimental  evidence  is 
presented  which  indicates  that  jet  noise  sound  power  spectra  can  be  adequately 
predicted  for  the  various  nosale  configurations  at  flow  pressure  ratios  up  to 
3.0,  and  flow  temperatures  up  to  1800°R.  Comparison  between  measured  and  pre¬ 
dicted  overall  sound  power  shows  that  the  prediction  method  used  provides  calculated 
levels  within  £  2  db  of  acoustically  measured  levels  for  nearly  all  nozzle  shapes 
and  flow  conditions  tested. 

In  Section  2,  the  development  of  analytical  techniques  for  jet  flaw  field  pre¬ 
diction  Is  undertaken  in  order  to  provide  a  working  tool  for  jet  noise  suppressor 
research.  Starting  from  the  approach  suggested  by  the  theory  of  Kelchardt, 
computer  programs  are  developed  which  represent  an  advancement  of  the  state-of-the- 
art  in  the  analysis  of  jet  flow  from  complex  nozzles.  While  it  is  shown  that 
analytically  derived  mean  flow  velocity  profiles  for  the  jet  mixing  region  in 

ix 


JL 


general  agrees  satisfactorily  with  experjasnt,  farther  work  is  required  to  bring 
computer  techniques  to  their  full  state  of  uaefullness. 

In  Section  3 ,  the  results  of  a  study  on  jet  noise  generation  and  suppression  at 
high  flow  temperatures,  up  to  those  representative  of  afterburning  jet  engine 
operation,  art  presented.  Experimental  data  on  noise  generation  by  jets  up  to 
flew  pressure,  ratio  of  3.4  and  tamperatare  up  to  3300°R  has  been  obtained,  in¬ 
dicating  that  for  low  pressure  ratios  and  moderate  temperatures  little  flow 
density  effect  exists,  while  over  the  full  range  of  observation  a  temperature 
dependent  flow  density  effect  is  found  to  occur.  This  effect  is  such  that  sound 
power  varies  with  (€j/ 6A]P  »  ''here  n  is  approximately  *  linear  function  of 

temperature.  The  results  of  seals  model  experiments  have  been  used  to  establish 
an  empirical  relationsnio  between  noksle  energy  flux  and  sound  power  generation, 
applicable  to  concentric  flows,  as  well  as  simple  oonioal  no  sale  a.  Studies  of 
suppressor  noaale  acoustical  performance  at  high  flow  temperature  have  bean 
carried  out  experimentally,  with  results  indicating  that  suppression  is  temperature 
as  well  as  pressure  ratio  dependent.  Jet  noise  suppression  at  high  flow  temperatures 
requires  further  study  for  the  purpose  of  deriving  generalised  temperature  scaling 
laws  applicable  to  all  nozzle  geometries,  since  the  experimental  results  presented 
indicate  that  temperature  effect  on  suppression  may  be  nozzle  geometry  dependent. 

The  results  of  the  high  temperature  jet  noise  investigation  point  out  the  need  for 
flow  temperature  simulation  in  suppressor  scale  model  development  directed  toward 
afterburning  jet  engines. 

Section  4  presents  ths  results  of  experimental  investigations  on  disr’.pation  of 
sound,  generated  by  a  jet,  in  the  turbulent  wake  of  the  jet  itself.  No  evidence 
of  sound  absorption  has  been  foutd,  and  thus  no  modification  of  the  correlation  of 
aerodynamic  and  acoustic  jet  properties  is  r squired  from  this  consideration.  The 
effect  previously  referred  to  by  other  researchers  as  jet  noise  "shielding"  is 
shown  to  be  the  result  of  directivity  rather  than  sound  power  change. 
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Selective  vater  injection  as  a  jet  noise  suppression  technique  has  been  ex¬ 
perimentally  sod  analytically  investigated,  with  results  presented  in  Section  5. 
Studies  indicate  that  tee  desired  alteration  of  jet  mixing  cannot  be  achieved 
through  United  vater  use,  and  that  tee  suppression  thus  accomplished  is  less 
than  that  expected  for  uniform  mixing  of  air  and  water*  Adverse  effects  of  high 
flow  temperature  and  velocity  on  noise  reduction  capability  of  vater  injection 
have  been  experimentally  demonstrated,  and  indicate  that  this  suppression  scheme 
is  not  practical  for  in-flight  noisa  reduction. 

The  knowledge  gained  in  tee  successful  completion  of  this  program  on  jet  noise 
suppression  research  can  be  Implemented,  along  with  other  existing  technology, 
toward  tee  development  of  a  jet  noise  suppressor  for  application  to  military 
aircraft  in  flight.  Thus,  noise  reduction  for  military  jet  aircraft  might  indeed 
be  a  practical  reality. 
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1.0  ’  Kxpqrlaggtal  Investigation  m  the  BeJatlor shin  Between  kexoSsm^R 

Properties  of  the  Jet,  Plow  Field  and  Kcsultant  Boise  Generation. 

1.1  Introduction 

As  s  result  of  a  previous  Kavy  sponsored  Geaersl  Electric  research  progr— 
on  Jet  noise  generation  and  suppression  (XOas  59-&L60-c)  an  empirical  f emula¬ 
tion  vas  developed  for  relating  the  aerodynamic  properties  of  a  free  Jet  to  its 
noise  characteristics  (Kef.  l).  The  hypothesis  that  the  sound  poser  generated 
et  various  axial  positions  in  the  Jet  nixing  region  depends  only  on  the  mean 
velocity  distributions  at  these  locations  led  to  a  general  quantitative  rela¬ 
tionship  for  determining  sound  power  spectra  for  siaple  conical  and  complex 
suppressor  nozzles .  It  vas  demonstrated  that,  vithin  certain  llaita  at  fl m 
pressure  and  temperature,  the  manner  in  which  the  flew  velocity  distributions 
very  vithin  the  Jet  because  of  nozzle  geoaetry  Variation  does  not  Influence  the 
generality  of  the  aero -acoustic  correlation. 

The  usefulness  of  the  previously  established  correlation  between  aero¬ 
dynamic  and  acoustic  Jet  properties  depends  to  a  great  extent  an  its  applica¬ 
bility  to  realistic  flew  temperatures  and  pressures,  and  to  various  nozzle 
configurations .  While  past  verification  of  the  theory  was  United  to  cold  flow, 
sub-critical  pressure  ratios  and  relatively  simple  nozzle  designs  consisting  of 
a  rusher  of  identical  elements,  the  present  study  considered  flew  temperatures 
up  to  lB00*R,  pressure  ratios  up  to  3*0,  and  large  variety  of  nozzle  designs, 
in  accordance  with  the  program  work  statement  in  Reference  2.  The  generality 
of  previously  established  empirical  constants  has  been  re-evaluated  through 
correlation  of  aerodynamic  and  acoustic  data,  indicating  that  the  velocity 
profile  formulation  can  indeed  be  applied  to  heated,  high  pressure  ratio  Jets, 
and  essentially  all  nozzle  geometries. 


* 
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1*2  Prediction  of  Sound  Foyer  Spectra  from  Hew  Profile  Data 
*  .  '* 

The  previously  established  relationship  between  aerodynamic  and  acoustic 

properties  of  a  free  jet  is  restated  frost  Reference  1  : 
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where; 


P  (f )  -m  Sound  power  per  cycle  at  frequency  f . 
X  ■  Empirical  constant 
-LdLgbfchill  parameter 


o  o  ao  o 
f 

n 

a 
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■  Frequency 

«  Local  Been  flow  velocity 

■  Flow  cross-section  area 

m  Axial  location  along  the  jet 
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The  assumptions  employed  are:  (l)  Ho  noise  is  generated  within  the  flow 
potential  core,  and  thus  the  integral 


/  ft/  -ft) 


must  be  evaluated  over  the  region  a  (x)  external  to  the  core  flew;  (2)  the 
axial  distribution  of  scale  of  turbulence  within  the  jet  is  linear,  that  is, 
L-^X;  (3)  the  relationship  between  axial  location  and  frequency  of  the  sound 
source  is  given  by:  f  d  Wl.25  x\  (2) 

u-*  V 

o 

where  d  ,  nozzle  equivalent  diameter,  is  defined  by: 
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and  a  it  the  number  of  floe  elements  la  the  nozzle;  (4)  the  total  generated 
sound  power  is  radiated. 

Prediction  of  the  sound  power  spectrum  is  accomplished  by  experimentally 
or  analytically  evaluating  the  integral  in  equation  (l)  at  a  sufficient  number 
of  x  -  locations  along  the  jet  so  that  the  frequency  range  of  interest  is  fully 
defined  (Equation  (2)).  In  this  section  of  the  report,  the  experimental  evalu¬ 
ation  at  jet  flow  properties  Is  discussed,  while  the  analytical  approach  is 
developed  in  Section  2.0.  — 

Under  a  previous  Navy-sponsored  research  program  (Reference  l)  the 
validity  of  this  method  of  sound  power  spectrum  prediction  was  demonstrated  for 

w 

an  eight-lobe  suppressor  nestle,  under  cold  flow  and  low  pressure  ratio  condi¬ 
tions.  For  coqplex  nozzle  geometries  employing  elements  of  non-uniform  size 
the  noxxle  equivalent  diameter,  d  ,  in  the  sound  source  location  equation 
(Eq.  (2))  must  be  suitably  defined.  The  definition  of  nozzles  equivalent 
diameter  as  previously  stated  in  Eq.  (3)  is  a  logical  choice  in  that  it  is  an 
average  elemental  diameter,  and  the  total  noise  is  generated  essentially  in  the 
region  where  the  jet  flow  exhibits  an  individual  element  identity.  Furthermore, 
the  jet  potential  core  length,  which  is  related  to  the  location  of  the  region 
of  maximum  noise  generation,  correlates  with  the  equivalent  diameter  as  herein 
defined. 

The  assumption  that  the  total  sound  power  generated  is  radiated  into  the 
far  field  is  considered  in  Section  4.0  of  this  report,  and  it  is  shown  that  no 
significant  dissipation  of  acoustic  energy  in  the  audio  frequency  range  occurs 
in  the  turbulent  jet. 
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In  order  to  determine  the  limits  of  applicability  and  to  further  develop 

the  correlation  described  above,  experimental  investigation  of  flew  velocity 

* 

profiles  was  carried  out  for  seven  nozzle  configurations,  flew  temperatures 

cp  to  l800*R  and  nozzle  pressure  ratios  up  to  3*0  were  tested,  and  aerodynamic 

data  vas  correlated  with  measured  sound  pewer  spectra,  as  described  in  the 

«• 

following  section. 

1.3  Experimental  Determination  of  flew  Velocity  Fri  files  and  Correlation  with 
Measured  £ound  Power  Spectra 

Aerodynamic  measurements  vers  mads  in  the  free  Jet  wake  for  a  number  of 
nossle  geometries  and  flow  conditions  for  the  purpose  of  defining  the  mean  flow 
properties  of  the  Jet.  flow  velocity  distribution  vas  determined  throughout 
tire  Jet  wake,  and  acoustic  properties  predicted  from  the  flow  measurements  were 
compared  with  measured  sound  pover  spectra. 

1.3*1  Aerodynamic  Measurements 

Detailed  flow  velocity  data  vas  obtained  throughout  the  significant 
noise  generating  regions  of  the  Jet  vake  for  the  various  nozzles  tested  using 
the  research  facility  shown  in  Figure  1.1 .  The  facility  consists  of  a  burner 
cf  the  type  conmonly  employed  in  aircraft  gas  turbines,  followed  by  a  flew  mixer 
and  pressure  and  temperature  instrumentation  for  accurately  determining  inlet 
conditions  to  the  nozzle  under  test.  Air  supply  capacity  for  the  facility  is 
up  to  9  lb/sec.,  and  flew  temperatures  up  to  1800*R  are  achieved  by  the  burner. 
Sozzle  pressure  ratios  up  to  3.0  were  possible  for  the  size  nozzles  tested  (up 
to  13  sq.  in.  exit  area)  vith  heated  flov.  For  high  pressure  ratio  cold  flew 
testing,  an  indoor  facility  vith  higher  airflow  capability  vas  enployed,  and 
test  instrumentation  used  vas  similar  to  that  for  hot  flov  measurement .  The 


test  nozzle  was  mounted  am  the  facility,  shown  in  Figures  1.1  and  1.2,  and  the  wake 
was  surveyed  with  a  combination  total  pressure  -  total  temperature  probe,  the 
position  of  which  was  remotely  controlled  by  means  of  probe  actuators  mounted 
from  a  lathe  bed.  During  the  research  pregram,  an  afterburner  was  added  to  the 
facility  (Figure  1*2)  increasing  flow  temperature  capability  to  3500*R,  but 
instrumentation  difficulty  precluded  velocity  profile  measurements  above  1800*B, 
and  the  addition  to  the  facility  was  used  far  studies  of  noise  generation  at 
high  tenperature  as  described  in  Section  3  of  this  report. 

Temperature  and  pressure  data  was  recorded  automatically  for  probe 
traverses  at  each  axial  location  in  the  Jet  wake,  and  point  by  point  calcula¬ 
tions  of  flow  velocity  were  accomplished  by  computer.  Since  the  velocity  in  the 
core  region  of  the  Jet  is  not  evaluated  in  the  prediction  of  sound  power  spectra, 
the  problem  of  measuring  supersonic  velocities  and  static  pressure  gradients  in 
the  care  region  of  high  pressure  ratio  Jets  was  bypassed.  For  the  pressure 
ratio*  under  consideration,  flow  velocity  in  the  noise  producing  mixing  region 
is  essentially  always  subsonic,  while  the  expansion  region  with  static  pressure 
gradients  and  shock  structure  is  assigned  to  the  Jet  core  In  which  no  noise  is 
generated.  Thus,  the  computation  of  flow  velocity  was  simplified  in  that  the 
assumption  of  static  pressure  equal  to  ambient  pressure  could  be  made. 

The  prediction  of  sound  power  spectra  from  velocity  data  does  not 
require  knowledge  of  the  velocity  distributions  at  a  given  axial  location  in 
the  Jet;  only  the  flew  velocity-flow  area  relationship  is  required  (Eq.  l). 
Therefore,  the  flow  velocity  data  was  reduced  to  the  form  of  equivalent  velocity 
prefiles.  That  is,  the  flow  was  artificially  reoriented  to  simulate  flow  from 
a  conical,  nozzle  for  the  purpose  of  simplifying  the  integration  process  and 


providing  a  comaon  means  of  presentation  for  all  of  the  flow  data  for  the 
various  configurations .  This  equivalent  data  is  presented  in  the  next  section, 
having  been  obtained  from  flow  velocity  naps  at  each  location  5n  the  Jet,  by 
graphically  accumulating  flov  area  in  the  order  of  decreasing  flow  velocity  to 
artificially  create  an  equivalent  conical  nozzle  profile  for  each  nozzle 
geometry  tested.  Fro*  the  equivalent  profiles,  the  integration  required  for  . 
evaluation  of  Equation  1  was  perforated  graphically. 

I.3.2  Acoustical  Measurements 

Aerodynaal  cally  predicted  sound  power  spectra  were  correlated  with 
actual  noise  measurements  obtained  for  each  nozzle  configuration.  The  outdoor 
test  facility  used  for  flow  velocity  measurements  was  also  employed  for  sound 
power  evaluation.  Boise  measurements  were  made  under  free -field  condition*.  A 
microphone  attached  to  a  remotely  controlled  boom  (Figure  1.3  )  was  used  to 
survey  the  sound  field  at  a  radius  of  15  feet  from  the  nozzle.  The  microphone 
a  Bruel  A  KJaer  Type  4l33»  has  flat  frequency  response  to  40,000  cps-  A 
frequency  analysis  of  the  Jet  noise  was  obtained  for  a  number  of  angular  loca¬ 
tions  from  the  Jet  axis  utilizing  the  instrumentation  shown  In  Figure  l4i  . 

From  the  directivity  data  thus  obtained,  so*”*!  power  spectra  were  computed  by 
the  process  of  integration  of  sound  pressures  over  the  radiation  area.  The 
computation  procedure  assumes  that  the  ground  plane  is  semi-absorbing  so  that 
m>v<nii  error  from  this  source  is  about  ldb,  as  the  difference  between  assuming 
a  totally  absorbing  and  totally  reflecting  surface  is  3  db.  The  location  of  the 
microphone  at  15  feet  from  the  Jet  permits  measurements  to  a  lew  frequency  limit 
at  250  cps.  The  measurement  radius  was  selected  so  as  to  provide  far -field 
conditions  for  lew  frequency  noise  and  yet  not  be  so  far  from  the  source  as  to 
oe  influenced  by  reflections  from  nearby  structures. 
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The  aerodynamic  and  acoustical  data  obtained  for  the  acres  nozzle 
geometries  tested  are  presented  is  this  section,  and  correlation  between 
predicted  and  Measured  sound  power  spectra  la  mad*. 


l.V.l 


In  order  to  determine  the  applicability  of  the  aerodynamic -acoustic 


correlation  to  various  nozzle  geometries,  a  number  of  configurations  were 
tested.  A  list  of  configurations  and  test  conditions  is  presented  In  Table  1*1. 
The  nozzles  are  all  of  approximately  the  same  scale,  with  exit  area  ranging 
from  9.6  to  13  la2.  The  geometric  variations  (as  shown  In  Figures  1.9  thru  1.8  ) 
include  simple  circular  exit,  annular  exit  (plug  nozzle),  shrouded  annular 
exit,  multiple  exit  (segmented),  and  multiple  exit  with  elements  of  unequal 
size.  The  variety  of  nozzle  configurations  provided  a  means  of  evaluating  the 
applicability  of  the  prediction  technique  to  essentially  ell  types^pf  suppressors. 


1.A.2  Correlation  of  Aerc 


lc  and  Acoustic  Data 


From  the  measurements  of  total  pressure  and  temperature  in  the  jet 
wake  far  the  various  configurations  and  test  points  (as  summarized  In  Table  l.l), 
mean  flow  velocities  were  calculated  and  flow  velocity  contours  were  plotted  for 


each  axial  location  within  the  jet  wake.  Velocity  contours  for  the  19-tube 
and  lB -segment  nozzles  are  presented  In  Figures  1.9  and  1.10  .  The  merging  of 
the  individual  flows,  which  are  quite  distinct  near  the  nozzle  exit,  can  be 
readily  observed  as  the  flow  progresses  downstream.  Equivalent  flew  profiles 
are  derived  from  the  velocity  contours  by  accumulating  the  flow  area  graphically, 
in  order  of  decreasing  flew  velocity,  resulting  in  the  profiles  as  shown  in 
Figures  1.11  through  1.16  •  The  equivalent  flow  profile  curves  enable  the  various 
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nozzles  to  be  considered  it  tens  of  a  fictitious  simple  circular  exit  geometry, 
but  with  flow  profiles  developing  differently  than  for  a  circular  exit  nozzle. 
The  normalized  equivalent  profiles  and  the  curves  showing  the  decay  of  maximum 
flow  velocity  as  a  function  of  axial  location  (Figure  1*1?  )  provide  the  data 
required  for  aerodynamic  prediction  of  sound  power  spectra  by  means  of  equations 
land  2. 

The  spectra  calculated  from  jet  flow  properties  are  compared  with 
acoustically  measured  spectra  In  Figures  1.18 thru  1,23.  Since  the  levels 
calculated  ore  power  spectrum  levels  (sound  power  per  unit  bandwidth),  the 
values  which  were  measured  in  percentage  handwiftths  have  been  converted  to  the 
same  constant  bandwidth  reference  for  ease  of  comparison.  Predicted  and 
measured  total  acoustic  power  levels  far  all  of  the  nozzle  configurations  and 
flow  conditions  for  which  power  spectra  have  been  determined  are  compared  in 
Figure  1.24,  where  the  differences  between  measured  and  predicted  levels  are 
plotted  against  flow  exit  velocity. 

Several  observations  can  be  made  concerning  the  results  of  this 
investigation.  The  shapes  of  the  predicted  sound  power  spectra  agree  very  well 
with  measured  values,  and  the  amplitudes  agree,  in  general,  within  the  limits 
of  experimental  accuracy.  The  spectrum  shape  results  imply  that  the  assumed 
distribution  of  sound  sources  within  the  jet  wake  is  adequate,  and  that  the 
evaluation  of  nozzle  equivalent  diameter  by  Equation  3  ia  vilid.  The  degree  of 
accuracy  achieved  in  the  prediction  of  absolute  noise  levels  (for  spectra  and 
total  acoustic  power)  verifies  the  applicability  of  the  fully  expanded  flow 
velocity  as  the  significant  nozzle  exit  velocity,  and  furthermore  indicates 
that  the  enpirical  constants  employed  in  the  correlation  are  applicable  over 
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the  range  of  flow  conditions  tested.  Although  somewhat  greater  amplitude 
deviation  than  can  he  attributed  to  experimental  error  is  noted  for  the  shrouded 
plug  nozzle  geometry,  it  is  believed  that  the  discrepancy  ir.  due  to  generation  of 
noise  (other  than  jet  noise)  by  the  shroud.  The  nanner  in  which  this  noise  is 
generated  has  not  been  investigated  here,  but  the  concision  is  in  agreement 
with  the  majority  of  past  experience  which  shows  that  addition  of  a  shroud  may, 
under  boom  circumstances,  result  in  production  cf  a  nusber  of  discrete  frequency 
noise  coaponents.  The  possibility  that  noise  generated  in  the  inner  regions  of 
a  complex  jet  is  dissipated  in  the  surrounding  turbulent  flow  is  more  fully 
discussed  in  Section  k.O  of  this  report.  However,  the  applicability  of  the 
aerodynamic -acoustic  correlation  to  the  variety  cf  suppressor  configurations 
evaluated  confirms  the  result  of  that  concurrent  study  that  no  significant 
absorption  tabes  place. 

1.5  Conclusions  and  Recoanendatlona 

The  results  of  this  investigation  further  confirm  the  hypothesis  presented 
in  Reference  1,  which  states  that  the  noise  generated  by  a  jet  is  uniquely 
related  to  the  mean-flow  velocity  distribution  within  the  wake.  The  technique 
for  calculating  sound  power  spectra  from  aerodynamic  properties  cf  the  jet  wake 
is  outlined,  and  the  generality  of  the  technique  is  demonstrated  for:  (l) 
nozzle  geometry  variations.  Including  circular  nozzles,  plug  nozzlea,  and 
suppressor  designs  whose  individual  elements  are  not  nil  equal  in  exit  dimen¬ 
sions;  (2)  flow  pressure  ratio  up  to  3*0;  and  (3)  flew  tccptrature  up  to  l800*R. 
Jet  noise  spectrum  shapes  and  levels  can  be  adequately  predicted  within  the 
limits  stated. 


9- 


Extension  of  the  applicability  of  the  jet  noise  suppression  theory  to 
actual  afterburner  teaperatures  (beyond  the  scope  of  the  present  program)  is 
discussed  in  Section  3.0  of  this  report,  and  it  is  probable  that  the  aerodynamic 
prediction  technique  can  be  applied  to  this  condition,  though  possibly  with 
8 one  modification.  Toward  that  goal,  the  analytical  techniques  for  flew  profile 

a 

prediction,  presented  in  Section  2.0,  could  be  employed  within  the- limits  of  * 

applicability  to  calculate  flow  data  for  correlation  with  noise  data  at  after- 

'  • 

burning  teaperatures. 

Finally,  the  equivalent  mean  flow  velocity  profiles  presented  far  the 
various  noise  suppressor  nozzles  may  suggest  other  geometric  designs,  more 
applicable  to  specific  problems,  which  would  generate  desired  profiles,  and 
thus,  low  noise  levels. 
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Figure  1.17*-  Maxima*  Velocity  Decay  in  Axial  Direction 
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Figure  1.17c-  Maximum  Velocity  Decay  in  Axial  Direction. 
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Figure  l.i7e*  Maxinmn  Velocity  Decay  in  Axiel  --rectloru 
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Figure  1.17f  -  Maxinun  Velocity  Decay  in  Axial  Direction. 
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Figure  l.lga  -  Comparison  of  P reclamed.  and  Meas'ured  Sound  Power  Specti-a 
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Figure  1.18b  -  Comparison  of  Predicted  and  "ensured  Sound  Power  Spectra 
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Figure  1.19  ~  Compari  on  of  Predicted  and  Measured  Sound  Power  Spectra. 
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p^r-^re  1.20“  Comparison  of  Predicted  ar.d  Measured  Sound  Power  Spectra. 
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Figure  1.21  -  Comparison  of  Predicted  and  Measured  Sound  Power  Spectra. 
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Figure  1.22b  -  Comparison  c F  Predicted  and  Measured  Sound  Power  Spectra. 
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1  igure  1,23  “  Comparison  of  Predicted  and  Ileasured  Sound  Power  Spectra. 
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C,°  ANALYTICAL  STUDY  ON  THE  HtEDICTIQN  OF  AERDDYANMIC  FROPERTTRR  nv  free  JET 
FLOWS  ' — . . . . 


2.1  '  Introduction 

In  the  preceding  section  of  this  report,  as  well  as  in  Reference  1, 
the  unique  relationship  between  aerodynamic  noise  generation  and  Jet  flow  field 
properties  has  been  demonstrated.  Thus,  from  knowledge  of  the  flow  field  the 
Jet  noise  sound  power,  spectrum  can  be  determined  by  analytical  means.  In  order 
to  relate  jet  noise  to  nozzle  geometric  parameters  by  the  previously  developed 
acoustical  theory,  a  means  is  first  required  for  prediction  of  the  flow  field. 

The  subject  of  free  jet  mixing  and  resultant  noise  generation  has 
been  the  object  of  extensive  experimental  and  analytical  research  for  a  number 
of  years.  In  the  field  of  free  jets,  the  most  complicated  practical  problem  is 
found  in  the  f  islysis  of  the  free  Jet  flow  emanating  from  sound  suppressor 
configurations.  It  has  been  found  through  research  that  good  sound  suppression 
performance  requires  a  significant  amount  of  geometric  complexity.  The 
analysis  of  the  complicated  jet  flows  produced  by  such  configurations  must 
depend  heavily  upon  a  very  simple  theory  to  be  at  all  manageable.  Such  a  theory 
is  the  approach  suggested,  and  used  by  Reichardt,  and  is  the  one  the  authors 
have  demonstrated  to  be  of  practical  usefulness  to  the  analysis  of  the  flow 
fields  from  jet  suppressors. 

In  this  study  the  development  of  analytical  techniques  for  jet  flow 
field  prediction  is  undertaken  in  order  to  provide  a  working  tool  for  suppressor 
research.  The  theory  is  applied  to  the  problem  of  predicting  aerodynamic 
properties  of  Jet  flows  from  suppressor  nozzles  through  the  use  of  a  computer 
program.  This  section  presents  the  analysis  and  computer  programs  that  have 
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been  developed  around  these  theories  toward  the  goal  of  analytically  predicting 
the  jet  flow  fields  for  complex  nozzle  geometries,  from  which  the  jet  acoustical 
properties  can  be  determined  using  existing  techniques . 

2.2  Theoretical  Analysis  and  Computer  Program  Development 

The  analytical  theory  and  method  presented  in  this  report  are  exten¬ 
sions  of  the  theory  end  method  originally  developed  in  Reference  1.  The 
original  work  produced  satisfactory  results  in  the  prediction  of  the  flow  field 
and  noise  produced  by  elemental  jet  types  (i.fi.,  jets  from  circular  and  rect¬ 
angular  exits).  That  work  was  culminated  by  a  computer  program  which  provided 
the  means  for  evaluating  the  theory  for  complicated  jet  flows.  It  is  in  the 
prediction  of  complicated  jet  flows  (such  as  from  the  nozzle  geometry  shown  in 
Figure  2.1  the  ;  uhe  orig  oil  worlt  has  been  found  to  be  marginally  satisfactory. 
Predicted  and  experimentally  determined  profiles  for  the  nozzle  shown  in  Figure 

2.1  are  compered  at  a  number  of  flow  cross  sections  in  Figure  2.2.  It  is  the 
task  of  the  study  reported  here  to  improve  the  original  work  in  order  that 
deficiencies  can  be  removed. 

2.2.1  Theoretical  Analysis  of  the  Jet  Flow  Field 

This  section  of  the  rt,  ■  ;sents  the  analysis  of  an  arbi¬ 
trary  jet  flow  field  in  which  the  solutions  for  velocity,  shear  stress,  and 
other  properties  of  the  jet  flow  are  considered  for  the  most  general  case.  The 
formal  and  detailed  derivation  of  the  equations,  which  supports  the  broad 
treatment  given  the  analysis  in  this  section,  are  presented  in  Section  2.6. 

The  principle  problem  area  of  this  study  has  been  the  descrip¬ 
tion  of  the  manner  in  which  jets  that  have  different  properties  (i.e.,  velocity. 
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density,  etc.)  intermix.  More  specifically,  the  problem  has  been  the  proper 
analytical  formulation  for  the  rate  of  Jet  spreading  (or  mixing)  under  condi¬ 
tions  of  relative  flow.  Once  the  mixing  process  is  accurately  described,  the 
determination  of  such  properties  as  velocity,  shear  stress,  and  density  is 
easily  accomplished. 

The  original  treatment  of  the  Jet  flow  problem  (Kef.  l)  used 
the  Reichardt  hypothesis  which  states  that  the  spreading  property  of  free  Jet- 
mixing  is  dependent  only  on  the  downstream  distance  (i.e.,  the  x-coordinate). 
This  hypothesis  was  used  in  the  Reference  1  analysis  even  for  complicated  flam 
such  as  those  emanating  from  sound  suppressors.  There  is  an  extensive  body  of 
literature  in  existence,  some  of  which  is  referenced  in  this  report,  which 
supports  the  theory  that  relative  velocities  determine  the  rate  at  which  mixing 
regions  spread.  Physically,  one  may  think  of  the  mixing  region  as  one  composed 
of  large  eddies  (of  size  comparable  to  the  mixing  region  itself)  which  cause 
the  transport  of  momentum,  energy,  and  concentration.  The  size  of  these  eddies 
is  a  function  of  the  shear  forces  forming  it  and  the  time  per  unit  distance 
available  for  their  formation  (convective  effect).  The  first  effect  is  often 
described  by  an  expression  involving  the  Jet  Mach  number  such  as  the  one  below: 

g  .  Cu  (X  .  Cb2  Mc)  <*•> 

do 

m 

This  expression  shows  that,  the  higher  the  Jet  Mach  number,  the  smaller  the 
eddy  size  at  a  given  station  ('lefs.  57  and  58)* 

Research  at  Stanford  has  shown  that  the  velocity  profiles  for 
simple  relative  jet  flows  can  be  predicted  by  using  a  simple  empirical 
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differential  equation  relating  the  grovth  rate  to  the  relative  velocity 
difference  (Ref.  2).  In  this  research  activity,  it  has  been  found  that  the 
relative  or  convective  velocity  effect  on  mixing  is  veil  approximated  by  the 
equation  given  below: 


Pl  +  °2 
Ul-P2 


(lb) 


During  the  time  when  the  Work  of  this  report  was  in  progress, 
no  information  was  available  as  to  the  manner  in  which  these  two  effects  were 
interrelated.  As  a  consequence,  for  the  purposes  of  this  study,  the  two  effect* 
were  geometrically  superimposed  to  give  the  following  equation: 


dx 

db 


Si*1 


Cb21'c> 


<V 


V 


(2) 


Previousc  researchers  have  also  found  this  relative  velocity  effect.  In  partlc- 

*  * 

ular,  Willis  and  Glassman  (Ref.  3)  found  a  strong  relative  velocity  effect; 
however,  they  were  not  able  to  derive  an  empirical  equation  describing  the 
effects  they  found.  Similarly,  Marshall  and  Bfeiley,  and  Weinstein  (Refs.  4  and 
5)  report  a  relative  velocity  effect  on  the  spread  of  the  mixing  region  hut 
again  do  not  develop  equations  describing  this  effect.  Kendall  (Ref.  l)  found 
that  the  centerline  velocity  distribution,  using  the  Eeichardt  method,  was  well 
described  using  the  rate  of  spread  given  by  Equation  (l)  for  the  case  of  zero 
external  flow  velocity. 

Section  2.2.3  gives  the  results  of  this  analysis  vising  the 
linearly  superimposed  solution  (Fq.  (2))  in  a  comparison  of  theory  am1  data. 
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For  simple  Jets,  the  relative  flow  velocity  (or  density) 


effect  produces  no  difficulty  so  far  as  the  theory  of  Reichardt  is  concerned. 

It  is  only  when  store  than  two  simultaneous  flows  are  concerned  (such  as  might 
occur  In  a  triple  annulus  of  three  different  flows  or  in  multi-element  Jet  exit 
shapes)  that  difficulty  with  the  Reichardt  method  is  found.  Before  delving 
Into  a  detailed  discussion  of  this  difficulty,  it  Is  perhaps  well  to  give  an 
example  showing  the  practical  significance  of  the  problem.  The  most  obvious 
example  occurs  in  the  superlmposltlon  technique  for  nozzle  suppressors  that  ia 
used  in  the  Reichardt  method  Itself.  Consider  the  suppressor  element  shown  in 
the  sketch  below  surrounded  in  part  by  Induced  flow  resulting  from  the  mixing 
action  of  the  Jet: 


Sketch  (a).  -  Suppressor  Flow  Approximation. 

In  the  approximation  technique  used  with  the  present  theory  (to  be  explained  in 
subsequent  sections),  the  induced  flow  is  conveniently  broken  up  into  arbitrary 
fictitious  Jet  flows,  each  having  a  unique  velocity  and  other  properties.  The 
relative  velocity  results  suggest  that  all  three  flows  (Vj,  V^,  V2)  spread  at 
different  rates  as  they  contact  one  another. 


In  the  multitube  type  of  suppressor,  a  similar  situation  arises 
as  the  mixing  regions  of  each  of  the  individual  tubes  interact.  Since  the 
interaction  is  not  uniform,  the  convective  velocity  and  hence  the  rate  of  mixing 
vary  in  all  dimensions. 

Another  example  is  a  multi -element  suppressor  having  both 
primary  Jet  flow  and  a  sound  suppressing  lcw-energy  secondary  flow,  all  embedded 
in  a  finite  external  flow  field.  Such  cases  as  these  cause  difficulty  in  the 
Reichardt  method  which  in  all  other  respects  has  proven  to  he  extremely  succes¬ 
sful  in  dealing  with  the  complicated  flow  fields  fround  in  suppressor  nozzle 
configurations. 


In  his  study  of  Jet  flows,  Reichardt  found  that  the  simplified 
constant  pressure  axial  momentum  equation  for  turbulent  flow 
2 

7)  T^r  +  1  2  (r  T^w)  .  0  (3a) 

r  ?r 

could  be  satisfied  by  an  elemental  exponential  equation  (as  ascertained  from 
analyzing  experimental  data  of  conical  free  Jet  mixing). 

2  2 

_P  lp«"rA»  (3b) 

if  certain  arbitrary  conditions  were  satisfied,  specifically, 

_ =  X  dTT  (3c) 

puv  2  dx  $  r 

In  this  elemental  relation,  there  appears  the  variable  b^  upon  which  the  rate 
of  spread  of  the  mixing  region  depends.  It  can  be  shown  that  the  differential 
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equation  for  axial  momentum  (Eq.  (3a))  is  satisfied  by  Equation  (3b)  if,  and 
only  if,  the  parameter  b^is  solely  dependent  upon  the  x- coordinate-  Here  the 
variable  b^  is  the  same  as  that  appearing  in  Equations  (la),  (lb),  oad  (2). 
That  the  variation  of  b^  in  the  other  coordinate  planes  does  lead  to  fallacious 
results  can  be  readily  shown  as  illustrated  below:  '* 


Sketch  (b).-  Y  dependent  b^  conflict. 

Consider  the  intermixing  of  two,  two-dimensional  jets,  each  spreading  at  a 
different  rate  (Sketch  (b)).  According  to  the  superposition  theory,  which 
results  from  the  linearization  of  the  differential  Equation  (3&),  this  momentum 
in  the  mixing  region,  shown  cross-hatched,  is  given  by  the  si®  of  the  momentum 
of  tne  individual  jets  (which  spread  according  to  Equation  (3b).  In  the  poten- 
tir.l  core  of  the  jet  undergoing  the  least  spread,  there  is  found  momentum  flux 
from  the  other  jet  which  is  spreading  faster.  This  region  is  shewn  by  the 
double  cross-hatching  in  Sketch  (b).  Thus  it  is  found  in  this  particular  region 
that  the  momentum  flux  is  the  sum  of  the  momentum  flux  of  the  potential  core  of 
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the  first  je t  plus  the  momentum  flux  spread  from  the  second  jet  -  and  this  sum 
is  greater  than  the  momentum  flux  from  the  first  jet  alone.  Such  a  result 
yields  velocities  in  the  potential  core  of  the  first  jet  (i.e.,  the  jet  that  is 
spreading  the  least)  which  are  greater  than  those  found  at  the  exit  plane,  a 
result  contrary  to  physical  reality.  This  discrepancy  represents  a  limiting 
condition  of  prediction  for  Jets  having  significant  core  interference. 

The  method  used  is  an  implicit  solution  for  the  x -coordinate 
using  the  spreading  parameter  bffl  as  a  dummy-dependent  variable  as  suggested  in 
the  original  work  by  Kendall  in  Reference  1.  In  this  technique,  the  properties 
of  the  jet  are  determined  for  various  b^  planes,  and  the  x- coordinate,  which  is 
now  specially  dependent  in  the  b  planes,  is  determined  using  Equation  (2) 
above.  In  the  bQ  planes,  the  momentum  and  energy  equations  are  uniquely  satis¬ 
fied,  and  the  contradiction  demonstrated  in  Sketch  (b)  (which  cannot  occur  in 
the  b  planes)  is  necessarily  avoided  but  the  continuity  of  the  x -plane  is  not 

IB 

guaranteed.  Even  if  the  x -plane  is  continuous,  the  satisfaction  of  momentum 
and  energy  conservation  principles  in  a  plane  of  constant  x  is  not  necessarily 
guaranteed. 


This  particular  problem  area  is  a  subtle  one  and  perhaps  not 
too  important  to  the  general,  class  of  free  jet  flows;  and  yet  it  is  of  utmost 
importance  to  the  successful  analysis  of  the  acoustic  suppression  phenomena 
found  in  suppressor  nozzles. 

A  closely  related  subject  to  the  foregoing  is  the  relative 
rate  of  spread  of  momentum  and  energy.  Experiments  of  previous  years  led  to  the 
well-known  result  that  energy  spreads  faster  than  momentum  in  free  jet  turbulent 
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mixing.  The  turbulent  Prandtl  number  which  relates  the  momentum  transport  to 
energy  transport  is  commonly  given  the  value  of  0.7  to  describe  this  phenomenon. 
The  rate  of  spread  of  energy  can  be  analytically  specified  in  a  manner  quite 
analogous  to  the  spread  of  momentum.  Reference  15  shews  that  the  Relchardt 
technique  linearizes  the  energy  and  concentration  equations  in  the  same  manner 
as  the  momentum  equation.  The  analysis  of  concentration  spreading,  in  this 
study  is  limited  to  concentration  spreading  at  the  same  rate  as  energy  (i.e, , 
turbulent  Lewis  number  of  unity).  Thera  is  a  spreading  parameter  b^  associated 
with  the  energy  spreading  as  there  is  a  b^  in  momentum  spreading.  In  this 
study,  the  analysis  was  extended  to  relate  b^  to  bffi  through  Prandtl  number. 

This  was  accomplished  using  the  definition  of  Prandtl  number  and  the  equations 
for  shear  stress,  etc.,  that  are  derived  from  the  Relchardt  method.  Several 
approximations  regarding  shear  stress  derivatives  must  be  made  to  achieve  the 

stated  result  which  has  a  simplicity  in  keeping  with  the  state-of-the-art  in 

♦ 

spreading  property  kifcwiedge.  The  approximate  derivation  is  presented  in  Section 
2.6.1  of  this  report,  and  the  result  is  shown  to  be  that  the  spreading  para¬ 
meters  are  related  by  the'  square  root  of  Prandtl  number: 

jm  a*  YPt  (4) 

This  result  has  been  programmed,  into  the  computer  solutions  contained,  in  the 
report. 


The  solution  for  the  properties  of  the  jet  flow  downstream  of 
the  exit  plane  has  been  extended  in  gene:  ility  from  what  was  developed  in 
Reference  1.  The  basic  technique  involved  is  the  summation  of  momentum,  energy. 
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and  concentration  flu*  using  the  sere  Integra?,  fora  as  that  developed  in 
Reference  1.  Only  the  general  method  cf  deveiqpaant  will  he  presented  here  as 
the  details  of  the  derivation  are  presented  in  Section  2.6.2. 


Using  the  superposition  technique,  the  momentum,  energy,  and 
concentration  flux  at  any  point  in  spec*  downstream  of  a  jet  exit  eontour(s)  can 
he  written  as  integrals  having  the  fera  shown  below  in  Equations  (f<),  (6),  and 

(7): 


2  *.2 

(1  -  e'r  K  )  a© 


^  ^ • f 

contours  _-/0 

2  2 

Em  hr  SL  (/^V)  •  f  (1  -  e‘r  ^0 

contours  0 

Ki  *  bf?  ^  e  (l  -  e"r  /^h) 

contours  " 


d9 


de 


(5) 


(6) 


(7) 


By  means  of  the  energy  equation,'  the  velocity  can  he  expressed  in  terms  of 
these  quantities  as  shown  in  Section  2.6.2,  The  resulting  expression  for 
velocity  is  a  quadratic  in  which  csly  the  positive  root  has  significance.  The 
equation  is  given  below. 


p  M  c—  T  ^ 

*  h_  +  (jr_  +  f  -*0  o  ) 

U  2 1  442  £ 


(8a) 


where 


£=  Mf  +  %  P 


(8b) 


'-SO'- 


Once  the  velocity  is  known*  the  other  properties  of  the  jet  such  as  temperature* 
density,  enthalpy,  etc.,  can  he  immediately  obtained.  The  solution  for  the 
velocity  and  other  jet  properties  has  been  generalized  to  the  point  of  handling 
n  number  of  different  gases,  each  having  its  owr  molecular  weight  and  specific 
heat  relationship.  This  means  that  the  solution  is  applicable  to  the  study  of 
hot  turbojet  exhausts  in  the  analysis  of  sound  suppressors.  The  only  restric¬ 
tion  that  is  made  on  the  thermodynamics  of  the  solution  is  that  the  Jet  flows 
at  the  exit  be  chemically  frozen. 


In  Reference  1,  it  was  found  necessary  to  include  two-dimen¬ 
sional  effects  in  the  equation  for  shear  stress  in  order  to  obtain  analytical 
results  comparable  with  experimentation.  At  that  time,  the  amended  shear-stress 
equations  were  not  included  in  the  computer  solution.  This  condition  has  been 
rectified  and  the  equation  for  shear  stress,  in  the  course  of  this  inclusion, 
has  been  refined.  The  exact  derivation  of  the  shear-stress  equations  stated 


below  is  given  in  Section  2.6.3.  For  the  radial  shear: 


contours 


X  (/>"g)  «  (0 »  e  -jfcrt  (r/bB) 

ntnurs  T?  J  ___ 


-TIT  erf  (  r  )  -  r  ~(r/V 
2  b  'Y2’  b  6 


and  a  similar  relation  for  the  circumferential  shear  stress  with  sin  e  replacing 
cos  0.  The  “total  shear  at  any  flow -field  point  is  given  simply  by  the  square 
root  of  the  sum  of  the  two  component  shears  squared.  A  brief  discussion  of  the 
results  and  comparison  with  the  test  data  is  given  in  Section  2.3.  It  is 
pertinent  to  note  here  that  the  computer  programs  can  be  used  as  research  tools 
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in  this  area  as  well  as  in  the  spreading  coefficient  problem  previously  mentioned. 
Here  the  effect  of  density  on  shear  stress  can  be  analyzed  and  the  correctness 
of  the  present  formulation  for  shear  stress  determined. 


2.2.2  Analysis  of  Jet  Flow  Acoustical  Properties 

In  this  section,  extensions  to  the  acoustic  theory  presented  in 
Reference  1  are  considered.  The  refinement  of  the  formulation  for  the  genera¬ 
tion  of  aerodynamic  sound  power  proceeds  from  the  point  ip  Reference  1  where 
equation  (A. 36)  is  presented.  The  more  involved  derivations  are  presented  in 
Section  2.6  of  this  report. 


In  the  present  work,  more  dependence  is  placed  on  the  Reynolds 
shear  stress  parameter.  This  is  a  favorable  development  in  that  a  suppressor 
is  likely  to  have  a  bigger  effect  on  the  variation  of  shear  stress  than  cm  any 


other  aerodynamic  parameter  in  the  jet  flow.  The  development  is  given  in 
Section  2.6.4  where  the  result  is  shown  to  be 
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The  analyses  of  this  report  show  that  the  maximum  power  is  generated  at  the 
terminus  of  the  potential  core  in  agreement  with  results  of  the  majority  of 
experimental  researchers  in  this  field. 


The  largest  deficiency  of  analytical  techniques  for  the  predic¬ 
tion  of  jet  acoustic  power  is  the  inability  to  accurately  predict  the  reduction 
in  acoustic  power  which  is  experimentally  found  with  the  multiple-exit  or 
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contoured  exit  nozzles.  This  acoustic  attenuation  has  been  found  to  depend  on 
both  geometric  and  aerodynamic  properties  of  the  jet  flow.  Tests  havs  shown 
that  sound  reduction  can  he  produced  even  with  simple  jets  by  placing  two 
simple  exits,  such  as  rectangles,  in  proximity  to  each  other. 

As  a  backup  method,  in  case  the  refined  aerodynamic  techniques 
described  above  do  not  suffice  to  give  the  proper  attenuation,  an  attenuation 
mechanism  and  an  analytical  method  have  been  established  in  this  study.  The 
most  obvious  mechanism  is  the  well-established  attenuation  of  sound  energy 
through  viscous  dissipation  and  is  the  one  used  here.  In  the  case  of  a  jet 
flow,  the  attenuation  is  presumed  to  occur  in  the  mixing  region  where  the 
apparent  turbulent  viscosity  is  the  greatest.  Such  a  mechanism  has  intuitively 
satisfying  characteristics  which  explain  several  of  the  experimentally  observed 
effects.  These  are: 


(a)  Sound  suppressors  which  seem  to  work  best  are  those  having 

inner  jet  flows  (e.g.,  the  multi-tube  supgpressar) .  % 

the  above  notion,  the  sound  energy  produced  by  the  inner 
jets  is  absorbed  by  the  mixing  regions  surrounding  these 
jets. 

(b)  As  the  jet  velocity  increases,  more  of  the  acoustic  energy 
is  passed  through  the  mixing  region  by  refraction.  This 
plus  the  greater  turbulent  viscosity  produces  increased 
attentuation  with  increasing  Mach  number. 

(c)  Attentuation  of  sound  energy  produced  by  suppressors  is 
known  to  be  frequency  dependent.  The  formulation  that  is 
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developed  in  this  study  is  found  to  have  frequency 
entering  to  the  second  power. 

However,  the  theory  of  the  viscous  damping  Mechanism  suffers 
in  two  areas.  The  first  and  most  important  is  that  the  phenomenon  is  not  in 
agreement  with  experimental  evidence.  Secondly,  the  mechanism  does  not  lend 
itself  to  mathematical  formulation  since  an  exact  analysis  would  detail  the 
history  of  the  wave  front  from  each  equivalent  source  in  three-dimensional 
space.  Because  of  the  first  deficiency,  such  an  exact  solution  is  not  Justi¬ 
fied. 

The  development  of  the  attenuation  phenomenon  presented  herein 
is  based  on  the  Kirchhoff  equation  for  plane  waves  (Bef.  6).  In  the  Klrchhoff 
development  used  in  this  report,  the  viscous  coefficient  is  replaced  by  the 

eddy  viscosity,  and  the  eddy  Prandtl  number  Is  substituted  into  the  equation. 

# 

Further,  a  measure  of  the  length  uve#  fSiich  the  attenuation  takes  place  is 
estimated  using  the  mixing  region  width  parameter,  bffl.  To  account  for  the  fact 
that  more  of  the  attenuation  would  seem  to  occur  in  the  innermost  portions  of 
the  jet  (a  logical  consequence  of  the  theory),  a  porosity  factor  is  applied  to 
the  attenuation;  this  porosity  factor  relates  the  attenuation  of  each  source  in 
the  acoustic  integration  tc  its  location  within  the  jet  flew  purely  on  a 
geometric  basis.  The  developments  of  this  theory  arid  the  equations  below  are 
given  in  Section  2,6,4.  These  equations  are 

(a)  The  attenuation  formula  is: 


(11a) 
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(b)  The  mean  length: 
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(c)  The  porosity: 
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The  above  set  of  Equations,  (10)  to  (13),  pl«8  Equation  (2),  contain  five 
arbitrary  constants,  the  determination  of  which  is  discussed  in  Section  2-3- 
Some  preliminary  results  obtained  vising  these  formulas  in  the  circular  exit 
solution  program  are  given  in  Section  2.3  of  this  report. 


2.2.3  Computer  Program 

The  equations  of  the  foregoing  section,  along  with  many 
subordinate  equations,  have  been  written  into  two  separate  programs  in  the 
process  of  the  development  of  the  solution.  The  need  for  two  programs  has  a 
historic  background  which  can  be  briefly  outlined  at  this  point.  The  formula¬ 
tion  of  the  equations  was  performed  early  in  the  study  contract  and  it  was  at 
that  time  obvious  that  the  determination  of  the  five  arbitrary  coefficients 
involved  in  the  solution  could  be  obtained  only  through  a  trial -and -erxoar 
procedure  using  a  computer  solution.  Rather  than  modify  the  original  program 
at  this  point,  it  was  decided  to  create  a  small  program  which  would  fulfill 
the  following  three  goals:  + 
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(a)  Provide  an  economical  means  for  obtaining  the  values  for 
the  five  arbitrary  coefficients. 

(b)  Provide  means  for  ascertaining  the  correctness  of  the 
relative  velocity  effect  on  sound  attenuation. 

(e)  Provide  program  logic  for  insertion  into  the  main  program 
during  its  revision. 

It  became  obvious,  however,  as  the  programming  of  the  "simple 
program*'  commenced  that  the  problem  had  become  complicated  to  the  point  where 
even  a  simple  circular  exit  solution  required  sophisticated  and  extensive 
programming,  especially  if  it  was  to  be  used  in  the  revised  main  program.  When 
it  became  obvious,  moreover,  that  the  simple  solution  was  too  large  for  the 
1620  computer,  several  of  the  first  goals  were  compromised  and  reliance, 
instead,  was  placed  upon  the  use  of  the  final  major  program.  The  circular  exit 
program,  then,  as  it  has  turned  out,  has  had  the  function  of  performing  (a)  and 
part  of  (b)  noted  above.  It  has  yielded  preliminary  values  for  the  five  arbi¬ 
trary  coefficients,  and  major  portions  of  the  program  are  being  used  in  the 
main  program  itself. 

The  simple  circular  exit  program  cannot  be  used  for  multiple- 
exit  configurations,  that  is,  double  annuli,  but  it  can  be  used  to  evaluate  the 
acoustic  attenuation  observed  when  a  secondary  flow  of  large  expanse  (compared 
to  that  of  the  jet  itself)  surrounds  the  circular  jet.  This  encompasses  by  far 
the  majority  of  the  experimental  data  available  on  jet  flows  including  those 
data  available  on  mixing  and  the  spread  of  mixing  regions  in  jets  surrounded  by 

finite  external  flow  fields. 
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In  the  following  sections,  the  philosophy  of  the  two  programs, 
some  of  the  programming  methods  used,  and  the  broad  physical  features  of  the 
pregrams  will  be  discussed.  In  Section  2.6.5,  a  more  detailed  description  of 
the  programs  is  given. 

Preliminary  Computer  Program 

The  solution  technique  of  the  simple  circular-exit  program  is, 
in  general,  parallel  to  that  of  the  main  program,  and  for  this  reason  it  will 
be  discussed  only  briefly.  In  the  simple  circular -exit  program^  the  inverse 
procedure  is  used  wherein  the  solution  is  first  performed  for  various  b  planes' 
and  then  x  is  explicitly  calculated  using  Equation  (2).  This  technique  offers 
the  advantage  that  the  flux  integrations  (which  are  b  dependent)  need  be 
performed  but  once  for  any  exit  configuration  (in  this  case  the  exit  radius); 
then  the  results  of  these  integrations  are  used  for  all  subsequent  running  of 
that  exit  configuration.  Such  a  technique  offers  a  factor  of  four  saving  in 
computer  time  compared  to  doing  the  integrations  each  time  the  program  is  run. 

In  the  main  program  a  similar  technique  is  used.  The  velocity  is  obtained 
through  an  iteration  procedur  it  is  in  the  main  program,  and  the  sound  power 
calculation  is  performed  explicitly  after  the  x's  (and  therefore  frequency) 
have  been  determined. 

Main  Computer  Program 

In  this  section  is  presented  the  philosophy  underlying  the 
computer  programming  of  this  project  and  an  outline  description  of  the  program. 
This  program  combines  the  geometric  generality  of  the  previous  computer  program 
reported  in  Reference  1  with  aerodyuamic  generality  and  sophistication  estab¬ 
lished  in  this  study.  The  aerodynamic  generality  stems  from  the  fact  that  the 
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program  can  handle  any  number  of  gases  haring  different  expositions  and  heat 
capacities.  At  present,  its  only  thermodynamic  limitation  is  that  the  gas  flows 
at  the  nozzle  exit  must  be  chemically  frozen.  Except  for  the  "moderate  compres¬ 
sibility  limitation"  on  the  development  of  the  shear-stress  equation,  the 
inherent  thermodynamic  limitations  placed  on  the  differential  equation  solution 
(i.e.,  the  original  simplified  differential  equation  due  to  Reichardt),  the 
program  makes  no  restrictions  as  to  the  density,  temperature,  and  the  Mach 
number  of  the  cases  it  can  handle-  Specifically,  the  program  as  it  cow  stands 
should  be  capable  of  providing  reasonably  accurate  solutions  to  actual  turbojet 
suppressor  problems.  The  only  aerodynamic  restriction  that  is  placed  on  the 
solution  is  that  the  overpressure  ratio  of  the  jet  at  the  exit  plane  be  nearly 
unity.  Actually,  the  program  can  be  used  in  certain  circumstances  far  the 
overpressure  case  and  for  providing  approximate  solutions  to  the  ejector  flow 
problem  through  clever  input  to  the  program.  These  techniques  are  discussed  in 
Section  2.4. 

The  biggest  stumbling  block  to  the  program  development  has 
been,  as  mentioned  pr«- /iousiy,  the  technique  for  the  relative  velocity  effect. 
Another  problem  has  been  the  development  of  a  method  for  dividing  the  flow  field 
so  that  the  integrations  in  the  sound-power  calculation  will  have  good  accuracy. 
The  program  now  integrates  200  points  (approximately)  in  a  typical  quadrant  of 
symmetry.  For  most  suppressor  configurations,  this  fineness  is  not  sufficient 
to  provide  accurate  sound-power  integrations  for  freq_enci.es  higher  than  those 
found  at  the  terminus  of  the  elemental  potential  cores.  For  example,  a  typical 
multi-tube  suppressor  would  require  about  an  order-cf-uagni tude  increase  in  the 
number  of  mesh  points  to  obtain  accurate  acoustic  integrations  at  the  higher 
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frequencies-  Such  an  increase  is  impossible  to  attain  in  the  present  program 
because  of  the  memory  storage  limitations.  It  is  thought*  however,  that  the 
200-point  solution  will  provide  sufficient  accuracy  to  at  least  determine  the 
maximum  power  output  of  these  sound  suppressors. 

The  program  calculate ,  flow  properties  and  acoustic  properties 
for  a  number  of  b  planes,  the  number  depending  upon  the  length  ratio  and  the 

SI 

acoustic  bandwidth  specified  in  the  input.  For  example,  if  data  are  desired 
from  the  exit  nlane  to  a  point  downstream  corresponding  to  x/R  of  about  30 

-  CjC 

and  the  acoustic  spectrum  is  to  be  in  one-third  octave  bandwidths,  then  about 

19  b  planes  (exponentially  spaced)  will  be  used.  This  technique  will  assure  a 
m 

high  probability  of  calculating  an  acoustic  point  for  each  acoustic  band. 

The  program  makes  a  "pass"  through  the  200  points  of  each  b^ 

plane  three  times.  On  the  first  pass,  it  calculates  the  flow-field  velocity  and 

a 

other  properties  such  as  density,  temperature,  ei£i  places  the  velocity  values 
into  memory  storage,  and  then  puts  it  and  other  property  values  on  the  output 
t-.ape-  On  the  second  pass  (currently  bypassed)  the  program  interrogates  the 
magnitude  of  the  velocity  of  neighboring  mesh  points  (surrounding  the  point  in 
question)  to  determine  whether  or  not  to  subdivide  the  existing  mesh  into  much 
finer  increments  to  provide  a  more  accurate  solution  in  a  way  which  avoids  the 
machine  memory-capacity  problem.  The  main  mesh  can  be  divided  in  care  of  three 
ways  1  a  horizontal  division  by  eight,  a  vertical  division  by  eight,  or  a 
combined  horizontal  and  vertical  division  by  eight  to  provide  6k  points  within 
the  existing  mesh.  The  size  of  this  inner  mesh  is  increased  as  the  program 
proceeds  to  calculate  b^  planes  further  downstream.  For  example,  at  an  x/R^ 
of  approximately  eight,  the  program  divides  the  major  meshes  into  16  points 
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instead  of  64,  the  general  idea  being,  of  course,  that  the  mixing  regions  at 
this  point  are  so  large  that  16  points  provide  as  much  accuracy  as  64  points 
provide  near  the  exit  plane-  There  is  also  a  practical  reason  for  this  decrease 
in  mesh  size.  As  the  program  proceeds  to  hffl  planes  located  further  downstream, 
the  number  of  meshes  increases  and  would  soon  exceed  the  capacity  of  the  7094 
memory.  %  increasing  the  mesh  size,  the  mesh  solutions  can  actually  he  nested 

in  the  machine  logic.  Figure  2.3  shews  a  b  plane  and  some  typical  major  mesh 

m 

points  having  inner  mesh  divisions  at  several  physical  locations  in  the  b 

m 

plane,  end  shows  how  these  inner  meshes  are  nested  to  conserve  memory  storage. 

On  the  thin1  pass,  the  program  interrogates  the  velocity 
points,  which  are  stored  in  memory,  within  "the  cone  of  influence"  from  the 
point  in  question.  This  cone  of  influence  has  a  physical  size  approximately 
equal  to  twice  the  size  of  the  typical  mixing  region  at  that  point  in  space. 
Figure  2.4  shows  the  cone  of  influence  interrogation  on  a  figurative  basis.  The 
program  determines  the  relative  maximum  and  minimum  of  the  velocity  points  in 
the  plane  within  the  cone  base  circle  and  uses  these  relative  maximum  and 
minimum  in  a  calculation  of  spreading  properties  of  the  jet  flow  and  the 
acoustic  properties.  From  these  properties,  a  unique  value  of  x  and  frequency 
is  obtained  for  each  of  the  major  mesh  points. 

At  each  point,  the  sound  power  which  is  calculated  is  assigned 

a  frequency  index  which  is  related  to  the  frequency  calculated  at  that  point. 

This  s<- -l 'd-power  calculation  (which  at  this  time  is  in  terms  of  energy)  is 

added  to  all  other  energy  calculations  in  the  plane  in  question  having  like 

frequency  indices.  These  frequency  indices  have  a  bandwidth  tolerance  of  one- 

tenth  octave.  At  the  termination  of  a  b  plane,  the  sound -power  calculations 

m 
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are  ordered  in  decreasing  levels  of  sound  power  and  those  sound-power  sums 
having  frequency  indices  within  one-third  of  an  octave  of  each  other  are  added 
and  the  frequency  of  this  sum  is  normalized.  Figure  2.5  presents  a  basic  flow 
chart  for  the  program,  outlining  the  discussion  just  presented.  It  in  no  way 
represents  the  entire  logic  of  the  program  which  is  far  more  extensive.  This 
flow  chart  is  keyed  to  the  program  listing  in  the  more  complete  discussion  of 
the  program  in  Section  2.6.5. 

2 , 3  Discussion  of  Res  -'Its 

This  section  will  discuss  the  results  that  have  been  achieved  in  this 
study  program.  Included  in  this  discussion  are  the  evaluation  of  the  coef  • 
f icients,  comparison  of  theoretical  velocity  profiles  with  the  data,  comparison 
of  theoretical  shear  stresses  with  data,  and  comparison  of  the  theoretical 
sound-power  production  with  data. 

Evaluation  of  Constants 

The  analytical  theory  developed  in  this  st-2y  involves  five 
arbitrary  constants  which  must  be  determined  i  rom  the  data  existing  in  the 
literature.  Two  of  these  coefficients  are  in  the  equation  which  relates  the 
rate  of  spread  of  the  mixing  region  to  the  axial  coordinate  and  the  convective 
Mach  number  (Eq.  (3))  and  one  of  each  are  in  the  equations  for  sound -power 
generation,  attenuation,  and  frequency.  In  the  evaluation  that  lias  been  conducted 
thus  far,  the  constants  have  been  found  to  be  coupled  so  they  cannot  be  deter¬ 
mined  independently,  thus  increasing  the  difficulty  of  obtaining  values  for 
the  constants. 
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The  constants  for  the  rate  of  spread  equation  have  been  ■ 
determined  largely  from  the  work  of  Lawrence  (Ref.  ?)  following  the  work  orig¬ 
inally  done  in  this  area  (Ref.  l).  Kendall  based  the  value  of  the  coefficients 
in  Equation  (l)  upon  the  spreading  of  the  flux  on  the  jet  centerline  where  the 
cylindrical  solution  is  completely  analytic.  In  the  present  study,  the  velo¬ 
city  profiles  for  all  space  have  been  compared  with  the  data  of  Lawrence  and 
the  best  fit  of  the  velocity  profiles  has  yielded  the  coefficient  value.  It 
has  been  found  that  the  coefficients  are  only  slightly  different  from  the  values 
presented  in  Ref.  (l).  These  values  are: 

cbl  =  13-5  (1%) 

% 2  *  1-5  (14b) 

The  above  constant  values  produced  reasonably  good  predictions  of  the  velocity 
profile  data  at  all  points  in  the  flow  field  for  the  two  cases  that  were  run 
(refer  to  Section  2.3).  The  acoustic  prediction  at  Mach  0.7  is  extremely  good 
concurrently.  However,  the  acoustic  prediction  at  Mach  1.05  was  found  to  he  2db 
lower  than  the  data  (Fig.  2.6).  This  result  is  directly  attributable  to  the 
reduction  in  mixing  predicted  by  Equation  (3)  for  the  higher  jet  Mach  number. 
With  C^2  set  equal  to  zero,  the  prediction  of  jet  noise  is  seen  to  be  much 
better. 

The  work  done  by  researchers  in  the  field  of  acoustics  has 
generally  supported  the  contention  that  the  frequency  of  the  sound  produced  at 
a  given  x  location  is  proportional  to  the  eddy  velocity.  On  this  basis,  the 
value  of  the  coefficient  would  be  expected  to  have  a  value  of  2.0.  Using  the 
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sound-power  profiles  such  as  Figure  A-U  in  Reference  1  and  the  results  of  the 
circular  solution,  the  value  of  the  coefficient  has  been  found  to  be 

CfSk.5  (15) 

These  same  data  have  yielded  values  for  the  sound-power  coefficient  and  the 
attenuation  coefficient  of 
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(16) 

250 

(17) 

If  the  analysis  were  precisely  correct,  these  values  for  the  coefficients  would 
apply  for  all  manner  of  configurations  and  conditions.  However,  since  the 
analysis  is  at  best  approximate,  it  may  be  found  j.n  the  course  of  running  other 
configurations  that  minor  changes  in  the  value  of  these  coefficients  will  yield 
better  overall  results . 


Velocity  Profiles 

Figures  2.7(a)  through  2.7(e)  present  the  results  of  the 
theoretically  predicted  velocity  profiles  determined  with  the  theory  presented 
herein  as  compared  with  the  data  of  Reference  7*  It  is  seen  that  the  theory 
indicates  less  spread  of  momentum  flux  at  both  high  and  low  values  of  the  axial 
distance.  Since  the  theory  conserves  momentum,  the  shape  of  che  velocity  profile 
suggests  that  the  data  are  in  error  since  the  data  definitely  do  not  satisfy 
the  conservation  of  momentum  condition  at  axial  distances  near  the  exit. 


The  comparison  of  theoretical  and  experimental  determination 
of  the  shear  stress  found  in  simple  flows  is  made  in  Figure  2.8.  It  is  seen 
that  the  theoretical  prediction  is  in  reasonable  agreement  with  the  data  (in 
keeping  with  the  prediction  of  Ref.  l). 

Sound-Power  Prediction 

Figure  2.6  presents  the  acoustic  predictions  obtained  with  the 
theory  and  the  values  of  the  constants  given  previously.  At  Mach  0-7*  the 
prediction  is  by  all  standards  excellent.  The  premature  dropoff  in  sound  power 
at  the  low  frequencies  found  previously  has  been  eliminated.  At  Mach  1.0,  the 
results  have  not  been  quite  as  satisfactory  for  reasons  suggested  before.  Bote 
that  the  difference  between  theory  and  data  is  quite  uniform  and  that  here  again 
the  theory  seems  to  he  predicting  the  correct  profile  and  frequency  shift. 

2 . 4  Additional  Approximations  for  :~joctor,  Overpressure,  acyWTnduccri  Flow 

Ejector  Approximation 

Because  of  the  problem  encruntered  in  the  mixing  rate  diffi¬ 
culty  discussed  in  Section  2.1,  the  extension  of  the  analysis  to  encompass 
ejector  (i.e.,  shrouded  nozzle  configurations)  was  not  performed  in  this  study. 
However,  sufficient  work  was  done  to  establish  an  approximate  procedure  that 
be  used  with  the  existing  program  to  determine  the  flow-handling  character¬ 
istics  of  ejector-type  sound  suppressors.  The  method  is  approximate  because  the 
theory  behind  it  assumes  constant  pressure  mixing  and  because  it  uses  a  ficti¬ 
tious  geometry  surrounding  the  primary  flow.  This  pseudo-geometry  is  shown  in 
Figure  2.9,  wherein  a  hexagonal  outer  nozzle  surrounds  the  primary  nozzle,  the 
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hexagonal  configuration  replacing  the  actual  circular  ejector,  (such  a  configu¬ 
ration.  is  of  little  use  in  those  sound  suppressors  utilizing  an  ejector  sur¬ 
rounding  a  complicated  exit  configuration  such  as  the  petal  suppressor  shown  in 
Sketch  (a) .  The  hexagonal  configuration  is  necessary  to  deceive  the  program 
into  generating  reflected  velocity  profiles  at  the  hexagonal  boundaries  which 
then  approximate  a  shearless  wall  velocity  profile.  Such  a  solution  for  simple 
ejector  geometries  should  produce  qualitative  flow-rate  ratios  for  those  ejector 
configurations  in  which  the  secondary  shroud  length  is  quite  short. 


Overpressure  Approximation 

The  overpressure  approximation  that  can  be  used  is  again 
restricted  to  simple  configurations  such  as  a  multitube  suppressor  and/or  to 
low  overpressure  ratios .  The  criterion  is  the  degree  of  flew  interaction  and 
the  attendant  shock  structure  strength  in  the  downstream  flow  field.  Such 
shock  structures  obviously  cannot  be  handled  by  the  existing  program.  The 
approximation  that  can  be  made  for  these  simplified  cases  is  to  increase  the 
exit "area  to  one  corresponding  to  the  exit  Mach  number  based  upon  one-dimen¬ 
sional  expansion  of  the  exhaust  gases  to  the  ambient  static  pressure.  The 
modified  area  coordinates  are  then  substituted  for  the  physical  exit  coordinates 
in  the  input  to  the  program. 

Induced -Flow  Approximation 

A  possible  cause  of  the  suppression  of  jet  noise  found  with 
complicated  jet  exit  configurations  is  the  inducement  of  a  fir-  in  the  mixing 
regions.  Sacn  flow  inducement  theoretically  reduces  the  shear  stress  and  there¬ 
by  the  noise  level.  To  auA  Licate  this  effect  analytically,  it  is  necessary  to 
input  to  the  computer  program  a  fictitious  exit  shape,  modifying  the  actual 
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suppressor  exit  shape.  In  the  modified  shape,  estimated  velocity  values  are 
used  to  produce  the  induced  flow.  Referring  to  Sketch  (a),  which  is  a  partial 
view  of  a  petal-type  suppressor  exit,  the  fictitious  exit  shape  that  is  added 
is  shown  by  the  dashed  lines  and  the  approximate  flow  in  this  portion  is  shown 
by  the  cross  hatching.  Such  a  technique  allows  the  computer  to  consider  the 
added  jet  flow  as  the  entrained  flow.  This  technique  remains  to  be  verified 
using  the  new  computer  solution  of  this  study. 


2.5  Detailed  Analyses  ar.d  Computer  Programs 


2.5*1  Derivation  of  Prandtl  Relation 

The  Prandtl  number  relates  the  transport  of  energy  to  the 
transport  of  momentum.  By  definition,  the  turbulent  Prandtl  number  is  given  by 
the  following  equation: 

Pr  -  cp?e  (13) 

ke ' 

The  definition  of  the  eddy  viscosity  and  heat-conduction  coefficients  in  this 
equation  are  as  follows: 


-  r_ 

d  u/<3r 


,  Pcp  v’h°* 

K  “  - 

e  J7/fr 

How  as  shown  in  Section  2.5«3»  the  shear  stress  is  given  by 

7-  K\i*K\  (^pS' 


dx  |  1 

Substitution  into  Equation  (19)  yields 
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dx 


(15) 

(20) 
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To  put  the  expression  for  kft  in  a  similar  form,  it  is  necessary  to  make  the 


following  approximations 


f>v» h°V  -  fVh°  (23 

which  will  hold  reasonably  well  for  moderately  incompressible  mixing.  By 
analogy  to  the  linearisation  of  momentum  flux  (which  leads  to  (21>,  one  can 


show  that 


Pv«h°  =  .  e)  (ftth5) 

2  1  dx  ]  ^r 


where  the  use  of  the  enthalpy  spreading  parameter  ^  can  be  noted.  Substitu¬ 
tion  into  (20)  produces  the  following  equation: 


_  hi  (dhh\ 

£em  p  2  \  dx  /  L 


3  tft*h°) 
3  x  _ 


Substituting  the  relations  for  the  transport  of  momentum  and  ir.ergy  (?2)  and 


(2$)  into  (ifi)  yields  the  following  equation: 


Vd  v 


Llr 


\  dbh 

dx  3 7 


Equation  (2)  shows  that  db^/dx  is  approximately  constant,  at  least  up  to  the 
end  of  the  potential  core.  The  presumption  made  here  is  that  the  derivative 


will  behave  similarly.  If  such  is  the  case,  the  derivative  ratio  is  propor¬ 


tional  to  the  r*Uo  of  the  spreading  parameters.  With  this  approximation  and 
further  expanding  the  derivatives  in  (26)  (neglecting  radial  gradients  of 
density),  the  following  result  obtains: 


m  A  *  *  .  %  *  . 


This  expression  relates  the  square  of  the  spreading  parameter  ratio  to  the 
Prandtl  number  and  another  term  containing  derivatives  of  velocity  and  enthalpy. 
This  latter  term  can  be  written  in  terms  of  logarithmic  derivatives 


2 _ _ 

1  +  3i  3-n  5L 

d(  In  h*)  _ 


(28) 


The  approximation  which  is  made  at  this  Juncture  is  that  the  logarithmic 
derivative  ratio  is  approximately  unity.  This  approximation  leads  to  the  final 
result  that  the  spreading  coefficients  are  related  to  the  square  root  of  the 
Prandtl  number. 


ja  •  (29) 

The  validity  of  the  above  approximation  would  seem  t®  he  most 
open  to  question  under  isoenergetic  mixing  conditions  where  theoretically  there 
is  no  gradient  in  stagnation  enthalpy. 

2.5.2  Derivation  of  Velocity  and  Related  Expressions 

At  any  point  in  the  jet  flew  field  the  flux  of  momentum, 
energy,  and  gas  concentration  ere  known  as  a  consequence  of  the  theory  used  in 
this  study  (see  Eqs.  (5),  (6),  and  (7)).  These  three  flux  quantities  will,  be 
termed  M^,  H,  and  K  in  the  subsequent  development.  It  can  be  shown  that  the 
mean  of  the  product  of  the  unsteady  quantities  in  these  flux  functions,  namely, 


A2 


H 


/°^r 


(30) 

(31) 

(32) 


can  be  well  approximated  by  the  product  of  the  means .  The  error  is  generally 
less  than  10  percent.  Reference  15  presents  an  exact  treatment  of  this  problem 
and  it  will  not  be  considered  further  here.  Also,  the  concentration  expres¬ 
sion  commonly  vised  for  the  mixing  of  two  dissimilar  gases  is  extended  in  this 
analysis  to  an  indefinite  number  of  gases  without  the  benefit  of  a  formal 
proof.  Thus 


ko-X 


(33) 


# 


where  kQ  is  the  total  concentration  of  ambient  gases  at  a  point  and  the  sum 
represents  the  calculated  total  concentration  of  jet  gases. 


From  the  energy  equation,  using  the  Crocco  relative  enthalpy 
for  boundary  condition  satisfaction,  one  obtains 


h  0  =  h  -  h  °  ♦  u2 
r  T 


(3*0 


Now  using  the  concentration  expression  given  above,  the  local  specific  heat 
for  an  arbitrary  gas  mixture  is  given  by 


c  -  c  k.  +  /l-  \  k.\ 

P  i  (  L>  iJ 


(35) 
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Taking  the  base  enthalpy  point  at  absolute  zero  reduces  the  above  equation  to 
the  one  used  in  this  study 
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Substitution  of  this  mean  specific  heat  into  the  energy  equation  produces 


r u  v*  r 


e  T  -  c 

p»  I 


T°  +  u2 


H 


(4l) 


The  following  algebraic  manipulations  and  substitutions  can  be  performed  to 
produce  a  quadratic  expression  in  the  desired  velocity  variable: 


p- 


H  -  u  c  P/Rg(k) 


u^/2  -  u  c  T 
pQ  o 


(42) 


(43) 


) 


(44) 

(45) 


In  this  quadratic,  only  positive  roots  have  physical  significance  nd  the  final 
equation  for  the  velocity  is 


u  -  H  + 

2T 


'h2 

IF 


M.  c  T 
+  f  P0  0 
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1/2 


(46) 
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where: 


M.  cp  P 
=  f  +  ym 
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(47) 
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For  the  general  esse  involving  gas  mixtures  pf  different  molecular  weights  and 
specific  heat  coefficients,  the  solution  is  iterative  around  the  variables, 
c  and  R  .  This  is  done  through  the  solution  for  concentration 

Ti 


the  gas  constant 


v  Vi  +  (1 
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and  Equation  (to)  far  the  mean  specific  heat,  c  .  When  the  molecular  weight 

pm 

and  specific  heat  are  constant  and  equal  for  all  the  gaBes  entering  into  the 
problem,  then  the  solution  does  not  involve  iteration.  In  running  the  programs, 
the  convergence  of  the  velocity  solution  has  been  rapid. 
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Shear-Stress  Derivation 

The  derivation  he a  as  its  basis,  the  representation  of  shear 


stress  by  the  well-known  equation  for  Reynolds  stress 


f  u'v* 


Now  it  can  readily  be  shown  that  the  mean  of  the  fluctuating  quantities  in  (l) 
can  he  replaced  with  differences  of  mean  quantities.  Thus, 


If  the  presumption  of  nods’-ete  incompressibility  made  in  Reference  1  is  made 
again  here,  then  the  results  of  (51)  substituted  into  (50)  yields  the  fundamental 
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relation  used  for  the  shear  stress  in  the  derivation 
t=  *  pxv  -  f>  u  v 


(52) 


Reichardt  found  a  relationship  between  the  first  term  of  (52)  and  the  radial 
gradient  of  momentum  flux  in  his  mixing  solution.  For  an  elemental  Jet,  this 
relation  is 


(53) 


The  momentum  flux  at  any  point  in  space  for  an  elemental  Jet  from  the  Reichardt 
solution  la  given  by 


P»2'  ( P»2)eAe« 


<50 


and  the  derivaiive  of  (5*0  with  respect  to  r  yields  the  following  relation: 

alfcLl  =  -2  (t Aere*r^  /v'b  '  (55) 

dr  /  * 

Substitution  of  this  result  into  (3)  gives  a  relation  far  the  first  term  in 
quantities  which  can  be  evaluated  through  a  ccoputer  integration. 

The  second  term  is  reduced  by  applying  the  continuity  equation 
to  the  flow  field  of  the  elemental  jet.  In  polar  coordinates,  this  expression 
is 

P  v  +  d  (P  v)  +  d  Pu  =  0  (56) 

r  r  dx 
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At  a  given  x-plpne,  the  gradient  of  mass  flux  with  respect  to  x  can  be 
considered  to  be  a  function  of  r  only.  Thus,  the  above  differential  equation 
becomes  a  simple  first-order  differential  equation  in  one  unknown,  having  as 
its  solution 


1  V  as  i 


0131 


where  the  constant  of  integration  is  zero.  Multiplication  through  by  u  leads 


to  the  result 


_ u 

^v  u  «  r 


3(Pu)  r’dr* 

?)* 


Using  again  the  elemental  Jet  momentum  flux  relation  (5*»)  and  taking  the  differ¬ 
entiation  outside  of  the  integration,  the  following  sequence  of  equations  can 
be  written: 


v  -  i  2~ 

r  *x 


r*  dr* 


- —  (  P  u  )  A 

It  V  .  11  '  1  't*  i 


Pm  v  =  u 

/  r. 


L  I  r*.  i 
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,2.  2  v 

r  r’dr’  \  (6o) 


Introducing  (5*0  (with  the  mean  of  the  product  assumed  equal  to  the  product  of 
the  means)  u  can  be  eliminated,  namely. 
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1r\r*  ix 


-r,2/b2 


i  '  r1  dr 


V/r 


Presuming  that  the  x-directed  gradient  of  density  is  nearly  zero  Justifies  the 


insertion  of  the  density  square  root  into  the  derivative.  The  density  expres¬ 


sion  may  b«  removed  altogether  if  it  is  presumed  that  the  local  density  is 


approximately  equal  to  the  mean  density  over  the  integration.  Such  simplifi¬ 


cation  reduces  Equation  (6l)  to  the  x'ollow'ing  (after  integration): 


I  FV - (p  L"  f 


-r2/2b  2 
e  '  m 


Performing  the  indicated  partial  differentiation  and  combining  terms  leads  to 


the  following  final  expression  for  the  second  term  of  the  shear-stress  equation 


for  an  elemental  jet  in  terms  of  quantities  which  can  be  evaluated: 


pr  uv=<^>eAe  fA\ 
ir\r  [  & 


2  2  /  2  .  ,  \  -r2/b  2 
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Addition  of  the  two  terms  produces  this  result  for  the  elemental  jet  shear 


stress: 


T-  (^u  }eAe 
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The  shear  stress  at-  a  point  in  the  flow  field  of  a  finite  jet  is  obtained  by 
the  linear  super  position  of  the  solutions  obtained  for  the  element  a!  jet. 

(The  exit  is  decomposed  into  a  matrix  of  elemental  jets.)  In  the  computer 
program,  the  super  position  is  replaced  by  an  integration.  The  radial  and 
tangential  components  of  the  shear  stress  are  obtained  by  multiplying  the 
integral  relation  by  the  cosine  or  sine  of  the  angle  between  lines  connecting 
the  point  in  question  to  the  origin  and  to  the  elemental  jets  being  integrated. 
(See  Fig.  2.10). 


Again  utilizing  the  super imposeble  nature  of  the  equations, 
integration  over  the  exit  plane  produces  the  total  shear  in  the  two  coordinate 
directions  (s  and  0)',. 


f  TP?)  it  -r2AZ 

/  e  _ El  e  '  m  cos  0 

Tfb  r’  dx 
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(65) 


Integrating  in  polar  coordinates  (dA  =  r'dr'dQ)  with  respect  to  r*  yields 

's 

contours 


Tie  =  /  ■[ ^  U^e  d  bm  f  cos  0  J  (r/b^) 
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and  a  similar  expression  for  0  with  sin  0  replacing  cos  6. 
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Sound-Power  Derivation 


2.5.4 


In  Reference  1,  the  following  equation  was  derived  from  the 
theory  of  Lighthill: 


Pa(x)~(AV  fe 


*2*) 


(67) 


This  equation  relates  the  sound  power  produced  within  a  volume,  a7,  about  a 
point  in  the  flow  field  to  the  eddy  oscillation  frequency  and  the  shear  stress 
generating  the  eddy.  The  volume  at  the  point  is,  as  in  Reference  1,  assumed  to 
have  as  its  characteristic  dimension,  the  mean  size  of  the  eddy,  CT .  Utilizing 
the  relation  below  (again  from  Ref.  l)  yields  the  following  equation  for  the 
volume  in  question.  The  frequency  equation  derived  in  Reference  1  is  also 
used  here. 


A  v~a3 


f  1  /tt\ 
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1/2 


The  product  of  these  two  quantities  is  then: 


AV  f  i  /  % 

x  *  »-' 
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(T  (  e'-3)  dV 


(68) 

(69) 


(70) 


By  placing  the  shear  stress  and  density  variables  into  the  integrals  of  (67) 
and  simplifying,  there  results 


(71) 
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This  volume  integral  can  be  replaced  -with  an  area.  integral  by  considering  the 


sound  power  to  be  produced  in  unit  length. 


dP 


TJ 


d  In  x 


Pa' 
I  o  o 


dA 


(72) 


From  the  frequency  relation 
t  ^  Uc/x 


(73) 


the  following  useful  equation  can  be  derived: 

d  In  f  rJ  d  In  Pc  -  1 
4  In  x  d  la  x  ' 


(7*0 


Substitution  of  the  last  relation  into  Equation  ^72)  yields  the  sound  power 
per  cycle  of  frequency  (or  bandwidth).  Introducing  the  attenuation  coefficient 
and  the  arbitrary  constant  into  the  equation  produces  the  equation  for  sound- 
power  production  used  in  this  study 


(75) 


The  acoustic  power  developed  by  free  jets  (sound  suppressors 
in  particular)  has  two  characteristics:  (l)  a  Kach  number  effect  and  (2)  a 
shielding  effect  wherein  the  sound  produced  by  eddies  internal  to  the  overall 
jet  structure  are  apparently  masked,  shielded,  or  attenuated  in  some  manner  by 
the  surrounding  flow  field.  An  exact  solution  would  describe  the  attenuation, 
refraction,  and  convection  effects  of  the  sound  produced  at  each  point  in  the 
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flew  field.  To  avoid  the  difficulties  of  performing  such  a  solution,  two 
parameters  will  he  artificially  introduced  into  the  standard  viscous  dissipa¬ 
tion  formula,  which  (it  is  hoped)  will  provide  the  desired  effects.  These  are: 
(l)  the  jet  porosity  which  simply  relates  the  point  to  the  jet  flow  field  in  a 
geometric  fashion  and  (2)  an  extent  effect  which  should  provide  a  measure  of 
the  time  the  acoustic  energy  is  under  the  influence  of  the  viscous  dissipation. 

Forming  an  attenuation  coefficient 

at*  e  "^kXe  (76) 


where  is  given  1/ 


Xe^  bm  Ef  Mc 


and  where  is  obtained  from  the  Kirchhoff  viscous  dissipation  equation 


P  oao 


A+  r  - 1\ 

v3 


This  equation  is  slightly  different  from  the  form  usually  found.  Here  the 
eddy  viscosity  replaces  the  kinematic  viscosity  and  has  been  factored  out  of 
the  expression  through  the  use  of  the  Prandtl  relation.  The  eddy  viscosity  in 
Equation  (78)  is  taken  to  he  some  mean  viscosity  over  the  mixing  region.  To 
obtain  this  effect,  the  following  arbitrary  equation  has  been  formed. 


(V  4  +  7. 


The  local  eddy  viscosity  is  given  by  the  standard  definition 
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and  using  the  approximate  relation  for  the  shear  stress  neglecting  density 
effects  produces 

_  pU  h  d  b  v 

•b.  =  r  c  a  _ a  (81) 

/e  2  dx. 


It  should  be  noted  that  the  above  relation  for  the  eddy  viscosity  is  an  approx¬ 
imation  for  the  elemental  jet  and  not  for  the  general  flow  field  which  would 
involve  integration  of  such  quantities  (as  in  Section  2.6.2). 

In  light  of  the  preceding  approximations  in  the  attenuation 
development,  such  an  integration  is  not  warranted. 

0 

«•* 


In  effect,  this  equation  arbitrarily  de-emphasizes  the  sound  energy  produced  by 
mixing  regions  close  to  the  jet  axis. 

Clearly,  the  use  of  such  arbitrary  expressions  reflects  the 
complete  lack  of  experimental  knc/wlsdge  in  this  area  which  could  direct  such 
analytical  efforts  and,  moreover,  uhe  complexity  of  the  situation,  which 
discourages  a  more  detailed  and  fundamental  treatment. 


Finally,  the  arbitrary  equation  for  the  porosity  effect  that 
has  been  tentatively  used  in  this  study  is  given  below: 


_  R’  A 
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2.5.5 


This  section  of  the  report  presents  the  description  of  the  two 
programs  that  were  developed  in  the  course  of  the  study  contract.  Of  the  two 
programs,  the  second  or  main  program  wiU.  receive  the  most  attention  since  it 
is  the  one  having  the  most  significance  to  the  analysis  of  suppressors  and, 
indeed,  is  the  desired  product  of  the  study.  These  programs  are  presented  in 
this  section  and  the  reader  is  asked  tr  'refer  to  these  sections  during  the 
descriptions  which  follow.  As  p.rt  of  the  program  description  typical  input 
and  output  for  the  two  programs  and  the  manner  in  which  the  programs  are  used 
will  be  discussed. 


Single  Circular  Exit  Program 


The  single  circular  exit  program  and  its  two  companion  sub¬ 
routines  are  presented  in  this  section.  The  first  of  these  sub-routines  (Sub¬ 
routine  Flux)  is  vised  once  for  a  given  exit  radius  and  the  second  subroutine  is 
used  thereafter.  The  exception  to  this  procedure  is  if  different  bandvidths  in 
the  acoustic  spectrum  are  desired  or  if  a  different  starting  point  in  the  bffl 
dummy  coordinate  is  desired.  These  cases  require  that  the  longer,  first -mentioned 
subroutine  be  used  (again  for  the  first  run  only).  The  running  time  for  the 
program  using  the  second  subroutine  is  a  factor  of  four  less,  so  that  there  i3 
a  considerable  advantage  to  using  it.  A  typical  running  time  for  the  program 
using  the  second  subroutine  is  about  1  hour  on  the  1620  computer  (evaluating  19 

b  planes).  The  second  subroutine  (Subroutine  Short  Run)  is  used  in  conjunction 
© 

with  the  output  from  the  first  subroutine  (the  output  used  is  a  part  of  the 
input). 
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The  following  brief  paragraphs  give  a  functional  resume  of  the 
major  portions  of  the  program. 

Preceding  the  exit-conditions  calculation,  input,  and  normal 
program  preliminaries,  is  the  sound-power  function  which  is  the  main  equation 
for  the  noise  prediction.  It  i3  basically  the  same  as  Equation  (lO).  In  the 
exit -conditions  calculation,  the  properties  at  the  exit  plane  such  as  velocity, 
temperature,  and  the  fluxes  of  mass  momentum  and  energy  are  determined.  The 
flux  values  are  necessary,  of  course,  in  determining  properties  in  the  jet  flour. 

The  determination  of  initial  conditions  is  primarily  for  the 
various  integrations  which  take  place  later  in  the  program.  Of  particular  note 
is  the  frequency  initial  ization -procedure.  Here  the  frequency  bands  are 
specified  and  the  sound  power  in  them  clear  _d  to  zero. 

Next  follows  the  determination  of  the  number  of  b  planes 
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which  will  he  calculated  and  the  specification  of  limits  for  the  various  loops 
in  program  logic. 

Entering  into  the  radius  loop  (wherein  the  points  along  a 
radial  ray  are  calculated)  the  flux  calculation  loop  is  first  found  and  it  is 
at  this  point  that  either  of  the  two  subroutines  to  the  program  are  called. 

From  the  flux  calculation,  the  velocity  is  determined. 

The  calculation  of  the  local  aerothermal  acoustic  properties 
is  followed  sequentially  by  the  calculation  of  sound  power.  The  sound-power 
calculation  procedure  is  more  complex  than  need  be  for  this  program  since,  in 
the  case  of  a  simple  circular  exit,  there  is  only  one  nixing  region  and  hence 
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only  one  characteristic  frequency  in  each  b  plane.  The  procedure  shown  is 

n 

used  in  the  main  program  where  more  than  one  characteristic  frequency  is  found. 

The  radios  loop  is  concluded  with  the  integration  of  the 
momentum  energy  and  mass  flux  in  the  h  plane.  These  integrations  are  used 
primarily  to  check  the  program's  accuracy  since  energy  and  momentum  are  theo¬ 
retically  conserved.  The  radius  loop  is  terminated  when  the  velocity  difference 
(jet  minus  free  stream)  approaches  zero  and  the  program  is  completed  with  the 
calculation  of  sound  power  in  decibels  and  the  final  output. 

Because  the  program  is  quite  simple,  the  only  procedure  that 
might  require  a  little  explanation  is  the  manner  in  which  the  two  subroutines 
are  used.  The  program  also  has  an  idiosyncrasy  in  the  subroutine  transfer 
which  also  bears  some  explanation. 

Presume  that  the  program  has  been  run  once  with  the  first 
subroutine  and  it  is  now  desired  to  run  it  with  the  short  flux  subroutine.  It 
is  necessary,  then,  to  take  the  output  cards  from  the  first  running  and  remove 
all  acoustic  output  data ,  This  is  easily  accomplished  since  these  data  axe 
readily  distinguished  blocks  of  output.  The  modified  output  should  then  consist 
of  a  sequence  of  two  output  data  cards  followed  by  a  blank  card  etc.,  from 
beginning  to  end.  In  this  form,  it  can  be  considered  input  for  the  short-flux 
subroutine  running  of  the  program.  Figure  2.11  shows  the  respective  positions 
of  the  input  as  it  would  be  loaded  in  the  •'ourse  of  running  the  program.  Notice 
that  the  modified  output-input  follows  the  regular  input  to  the  program. 

A  peculiarity  may  be  found  in  running  the  short  flux  program. 
Since  the  output  does  not  have  the  significant  figure  accuracy  carried  in  the 
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machine,  the  termination  of  the  radius  loop  previously  mentioned  does  not  occur 
for  precisely  the  same  radius  with  the  second  subroutine  because  of  a  minor 
difference  in  velocity  values.  This  naturally  destroys  the  correct  sequence 

of  operation  if  the  machine  overruns  on  a  b  plane  calculation.  The  machine 

n 

continually  checks  the  sequence  so  that  such  an  error  automatically  halts  the 
machine.  To  correct  this  situation  (which  is  found  to  occur  only  for  the  last 
few  1>  plane  calculations),  it  is  only  necessary  tc  set  (artificially)  the  flux 
value  in  the  output  to  zero  on  the  last  card  of  the  particular  troublesome  b^ 
plane  output,  (in  Figure  2.13  of  the  typical  output,  these  quantities  are 
circled). 


Figure  2.12  shows  some  typical  input  to  the  program,  which, 
for  the  most  part,  is  largely  self-explanatory.  The  units  are  in  ft-lbs-seeonds 
for  all  quantities  including  those  for  specific  heat  and  the  gas  constant.  As 
the  program  now  stands,  the  proper  value  for  the  characteristic  radius  is  the 
same  as  the  radius  itself.  The  output  is  presented  in  Table  I  which  again  is 
largely  self-explanatory.  The  program  has  been  modified  so  that  the  decibel 
calculation  is  performed  for  as  many  bandwidths  as  in  the  output  to  the  b^ 
planes  (19  in  the  case  shown  in  Figure  2.13). 
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Main  Program 

The  Compressible  Flow  Mixing  Profile  and  Acoustics  Program  and 
its  subroutines  (presented  in  this  Section)  is  a  greatly  expanded  version  of  the 
original  computer  program  developed  in  the  study  reported  in  Reference  1.  The 
major  contributor  to  this  expansion  has  been  the  logic  required  by  the  routine 
for  the  determination  of  the  spreading  characteristics  of  the  jet  flow  under 
relative  conditions,  and  the  routines  required  to  avoid  exceeding  the  xneaoxy- 
storage  capacity  cf  the  machine.  Much  of  the  flux  integration  logic  for  compli¬ 
cated  exit  shapes  that  was  developed  in  the  original  program  has  been  retained 
and  can  be  recognized  in  the  major  subroutine  to  the  program  -  Subroutine  FLUX. 
It  was  originally  thought  that  computational  time  could  be  significantly 
reduced  by  integrating  the  flux  from  simple  exits  with  simplified  methods  such 
as  used  in  the  Single  Circular  Exit  program.  In  particular,  the  intent  was  to 
calculate  circular  and  rectangular  exit  shapes  in  this  fashion.  It  has  turned 
out  though  that  only  the  circular  exit  can  be  so  treated,  due  to  the  complexity 
of  integrating  the  shear  equations  for  the  rectangular  exit  in  closed  form. 
Subroutine  CFLUX  is  this  routine.  It  is  somewhat  different  than  the  subroutine 
used  for  the  Single  Circular  Exit  program,  in  that  the  former  recognizes  the 
circular  exit  as  a  point  source  at  sufficient  removed  distances  from  the  exit 
(which  again  increases  the  speed  of  the  program).  For  this  reason,  subroutine 
CFLUX  will  provide  greatly  reduced  computational  times  for  suppressor  exit 
shapes  such  as  the  multiple  tube  nozzle. 

The  following  discussion  will  consider  first  the  main  program 
and  its  self-contained  subroutine  and  then  the  principle  subroutines  to  the 
program  itself.  As  in  the  preceding  program  discussion,  the  emphasis  will  be 
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placed  on  the  functional  aspects  of  each  of  the  major  portions  of  the  program 
and  little  or  no  mention  will  he  made  of  the  programming  logic  involved. 

Skipping  past  the  two  inner  subroutines  for  the  moment,  the 
first  section  of  technical  importance  are  the  exit  area  and  exit  conditions 
calculation  procedures.  In  the  first,  the  total  exit  area  of  the  nozzle 
configuration  under  consideration  is  calculated.  In  addition,  the  area  of  each 
of  the  elements  making  up  the  nozzle  (as  in  a  multiple  tube-type  suppressor 
nozzle)  is  stored  for  later  use  in  the  flux  integral  checks.  In  the  conditions 
calculation  procedure,  there  are  two  techniques  and  one  iteration  involved. 

The  iteration  is  necessary  because  in  this  program  only  the  exit  stagnation 
pressure  and  temperature  are  specified,  and  all  properties  are  obtained  from 
these  two  properties.  In  the  case  where  variable  gamma  is  being  considered,  an 
iteration  on  the  energy  equation  is  required  to  determine  the  exit  static 
conditions.  In  the  first  method  (constant  and  identical  gamma  and  molecular 
weight  for  all  gases  of  the  system),  the  iteration  procedure  is  not  required 
since  the  exit  static  properties  can  be  obtained  from  the  standard  one-dimen¬ 
sional  Mach  number  relations. 

As  in  the  previous  program,  certain  initial  calculations  and 
zero  settings  must  be  performed  and  these  are  found  next  in  this  program  es 
well.  The  frequency  initialization  is  much  the  same  as  before.  The  velocity 
is  set  equal  to  the  free -stream  value  to  assure  correct  interrogation  results. 

The  determination  oi  the  indices  on  the  major  loops  in  the 
program  logic  follow  the  initial  condition  determinations,  and  this  in  turn  is 
followed  by  the  major  loop  of  the  program  which  determines  the  set  of  b  planes 


in  which  the  values  for  the  acoustic  and  aerodynamic  properties  at  the  Jet  flow 
are  calculated. 

It  is  at  this  point,  that  the  control  for  the  b  plane  pass  is 

found.  The  program  makes  three  passes  in  each  b  plane  as  discussed  in  Section 
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2.2.3.  The  pa ss  is  controlled  by  the  value  of  a  special  program  signal  term  - 
in  this  case  the  coded  word  is  IBPAS.  In  every  pass,  the  program  logic  proceeds 
to  the  flow-field  matrix  which  determines  the  two-dimensional  point  matrix  in 
the  flow-field  b  planes.  These  points  are  determined  by  means  of  radius  and 
angle  loops  which  specify  points  approximately  at  the  corners  of  rectangles 
(refer  to  the  flow  chart.  Pig.  2.5). 

The  Section  entitled  "No  Mesh  Set  on  NMDOM"  is  artificially 
included  in  the  program  on  a  temporary  basis  to  insure  that  the  second  pass  of 
the  program  is  bypassed.  As  noted  previously,  this  section  must  be  bypassed 
because  the  companion  subroutine  to  the  matrix  subdivision  has  not  yet  been 
completed. 

The  flux  loop  procedure  follows  next,  and  it  is  at  this  point 
that  the  subroutines  are  called  into  the  main  program.  The  various  fluxes  are 
summed  over  the  number  of  individual  elements  that  exist  in  the  suppressor 
configuration.  Several  bypasses  are  included  in  the  flux  calculation  to  reduce 
the  program  running  time.  The  first  of  these,  entitled  "X  zero  shear  and 
iteration  bypass",  is  used  at  the  first  b  plane  calculation  which  is  coplanar 
with  the  ex.it  of  the  suppressor  nozzle  and  which  13  used  to  associate  the  exit 
velocity  with  corresponding  mesh  points  in  the  first  b  plane.  The  second 
bypass  is  on  the  shear  and  velocity  iteration  calculations  which  are  not  needed 
for  those  matrix  points  lying  within  the  potential  core  of  the  Jet  ?fflux. 


The  velocity  iteration  which  follows  is  the  sane  as  that  used 
in  the  preceding  program  with  the  addition  of  a  CP  and  PGC  (perfect  gas  constant) 
loop  which  provides  a  solution  for  gases  of  different  molecular  weight  and 
variable  specific  heat.  Two  coefficients  are  used  to  describe  the  functional 
relationship  between  specific  heat  and  temperature. 

The  next  major  section  of  importance  is  the  interrogation  loop 
wherein  those  points  lying  within  a  circle  whose  radii  are  twice  the  width  of 
the  mixing  region  are  compared  to  determine  the  relative  maximum  and  minimum  of 
this  set  of  points.  Since  the  major  mesh  points  are  not  subdivided  in  the 
present  program  the  inner  mesh  points  Cannot  he  interrogated,  and  for  this 
reason  this  portion  of  the  program  is  now  bypassed  as  indicated  in  the  listing. 
Once  the  relative  maximum  and  minimum  velocities  in  a  given  mixing  region  are 
known,  then  properties  such  as  the  convective  velocity,  the  local  value  of  the 
x-coordinate,  the  eddy  viscosity,  and  the  acoustic  properties  can  be  determined. 
The  attenuation  calculation  procedure  is  identical  to  that  used  in  the  preceding 
program,  and  only  slight  modifications  have  been  cade  to  the  x  determination 
and  the  eddy  viscosity  calculation  procedures. 

Tue  matrix  subdivision  section  of  the  program,  which  is  by¬ 
passed  for  reasons  previously  discussed,  contains  three  separate  routines 
depending  upon  the  manner  in  which  the  major  mesh  points  are  subdivided. 

The  first  two  sections  beginning  with  statements  numbered  167 
and  168,  respectively,  divide  the  major  matrix  either  circumferentially  or 
radially,  and  the  section  beginning  with  statement  169  performs  the  division 
both  ways.  In  all  three  cases  the  actual  subdivision  is  carried  out  by  one  of 
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the  self-contained  subroutines  (subroutine  mods).  This  subroutine  tells  the 
main  program  what  the  mesh  size  is  in  the  subdivision  -  for  hffl  planes  near  the 
exit  this  division  is  8  by  8:  fear  b  planes  near  the  end  of  the  potential  core 

n 

the  division  is  by  4;  and  downstream  of  this  point  the  division  is  by  2.  For 
b  planes  three  diameters  downstream  of  the  nozzle  exit  plane  there  are  no  inner 
mesh  sub-divisions. 

Following  this  calculation  procedure,  are  the  ratio  calcula¬ 
tions  for  the  output,  sound-power  calculations  (which  are  similar  to  those  in 
the  preceding  program),  the  mass  momentum  and  energy  integrations,  and  the 
termination  of  the  mesh  point  angle  and  radius  loops.  Next  is  found  the 
acoustic  search  and  addition  routine  which  is  another  major  portion  of  the 
program  logic  and  the  one  which  uses  the  other  internal  subroutine  for  the 
program  -  subroutine  orders.  It  is  here  that  the  program  searches  through  the 
frequency  indices  that  have  been  determined  at  a  given  plane  (one  frequency 
index  being  associated  with  each  mixing  region  having  different  relative  maximum 
and  minimum  velocities).  Those  sound-power  values  having  frequency  indices 
less  than  one-third  octave  apart  are  added  together  to  form  the  acoustic  output 
of  the  program.  The  band  assignment  and  shift  routine  on  the  following  page  is 
currently  not  used  in  this  program  due  to  memory  limitations.  In  it,  the 
machine  compensates  for  the  fact  that  the  sound  power  is  being  calculated  at 
discrete  points  and  not  integrated  throughout  the  whole  flow  field.  In  the 
band  assignment  and  shift  routine,  the  program  searches  through  the  bands  that 
have  been  specified  in  the  output  after  the  final  b^  plane  has  been  calculated 
to  make  sure  that  each  major  frequency  (i.e.,  the  major  mixing  regions)  contri¬ 
bute  to  each  bandwidth  once  and  only  once  after  they  appear  and  before  they 


disappear.  Finally,  all  the  sound  power  from  the  various  mixing  regions  a»-e 
added  at  each  bandwidth  to  give  the  total  sound  power  generated.  At  present, 
since  the  band  assignment  and  shift  is  not  operational,  this  function  will  have 
to  be  performed  by  hand,  using  the  output  from  the  acoustic  bands  ear  ch  and 
addition  routine. 

Following  the  main  program  listing,  is  the  listing  for  sub¬ 
routine  FLtiA  and  its  subordinate  subroutines.  This  subroutine  is  a  complicated 
bit  of  programming  logic  which  is  necessary  to  avoid  exceeding  7090  machine 
memory  storage.  It  also  serves  to  keep  the  cost  of  running  the  program  down  to 
a  reasonable  level.  In  this  subroutine,  the  geometric  integrations  that  are 
performed  for  each  matrix  point  are  stored  along  an  axial  my  so  that  for  each 
position  in  space  the  geometric  integrations  need  to  be  performed  but  once. 

Since  there  are  approximately  19  b  plane  calculations,  this  technique  alone 
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presents  a  19  to having  in  computational  time.  This  subroutine  also  stores 
the  value  of  the  flux  integrations  themselves  so  that  subsequent  runs  with  this 
particular  geometry  can  be  made  without  re-doing  the  time-consuming  integrations 
which  (as  in  the  preceding  program)  will  result  in  a  considerable  time  reduction. 

Despite  its  complexity,  the  main  program  is  easier  to  use  than 
the  first.  Any  number  of  cases  can  be  considered  at  the  same  time,  since  the 
program  uses  the  correct  subroutines  automatically.  Thus,  any  sequence  of 
nozzle  configurations  and  aerodynamic  cases  for  each  configuration  can  be 
handled  in  any  sequence  desired.  It  is  important  that  a  correct  case  number 
for  each  configuration  be  included  since  it  is  the  case  number  which  tells  the 
machine  which  subroutine  it  is  to  use.  For  instance,  case  Ho.  1  is  reserved 
for  the  first-run  integrations  which  are  stored  on  a  master  rape.  This  master 
tape  must  be  used  for  all  subsequent  running  of  this  configuration. 


2.6 


Conclusions 


In  the  course  of  conducting  this  study,  several  significant  areas 
worthy  of  further  analytical  and  experimental  research  were  uncovered.  These 
areas,  of  course,  are  refinements  to  the  existing  knowledge  and,  as  a  conse¬ 
quence,  will  he  less  rewarding  for  the  sane  amount  of  work  than  has  been  exper¬ 
ienced  in  the  past*  This  is  particularly  true  in  the  one  area  that  has  caused 
the  nost  concern  in  the  present  atudy  -  the  relative  velocity  effect  on  the 
rate  of  spread  of  jet  adzing. 

The  conclusions  of  this  study  are  as  follows: 

(1)  Reinhardt's  adzing  theory  is  not  adequate  in  fully  describing 
the  relative  rate  of  adzing  in  complicated  flow  situations. 

(2) .  The  jet  Mach  number  effect  observed  primarily  in  the  data  of 

Reference  10,  are  in  contradiction  to  the  observed  acoustic 
results  from  Reference  1.  (Presuming  the  present  aconstic 
theory  is  essentially  the  correct  one) .  As  a  consequence,  the 
value  for  the  second  rate  of  spread  coefficient  (constant  C^) 
is  found  to  be  zero  at  transonic  jet  Mach  numbers. 

(3)  Experimental  evidence  indicates  that  the  rate  of  spread  in  the 
mixing  region  is  to  some  degree  x  dependent.  This  effect  is 
not  included  in  the  present  analysis,  and  the  acoustic  predic¬ 
tions  for  circular  exits  indicate  that  such  a  dependence  is 
desirable  for  a  more  accurate  prediction. 


(ii)  The  modified  shear-stress  equations  successfully  predict 
experimentally  measured  shear  stresses. 

(5)  The  analytical  Telocity  profiles,  in  general,  satisfactorily 
agree  vith  experiment.  In  this  regard,  it  is  noted  that  the 
experimental  data  of  Reference  7  seem  to  have  a  deficiency  in 
that  the  data  are  at  odds  vith  the  conservation  of  momentum 
principle. 

(6)  In  general,  the  computer  programs  presented  herein  are  an 

* 

advancement  in  the  state-of-the-art  in  the  prediction  of  fluid 
dynamic  and  acoustic  properties  of  interacting  free  jets. 
Further  work  is  required,  however,  to  bring  these  programs  to 
their  full  state  of  usefulness. 
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2  READ  701.  PRAOL  .  CX  .  CA  «  CFREO  «  CB1  .  G32  *P«3XpR  «3M INT  .  (GAMAE  (L  )  ,  RGE  (L  J  « 
ICPcCLl  .TP£(L).TTE<L1.L«  1.2) 

FREO  *  G. 

PI  a  3.1M59 

TP  I  =  6 . 23318 

AREAC  =  PI *RADC*RADC  ' 

RAOL  a  RAOC 
SEMAS  *  0. 
sekom  "0. 

ScNLP  a  0. 

EANCZ  a  0. 

EXIT  CONDITIONS  CALCULATION 

N<P  =  NMNOZ  +  1 
DO  601  <-  1  .  NKP 

GAMEX  a  (GAMAE (K )  -  1 .) /GAMAE CK ) 

EMACH  =  2.* < ( (7PE <K >/P)**GAMEX)  -1  .)/< GAMAE CK >-l . ) 

TE  =  TTECK)  /  Cl.  -r  ( GAMAE  CK)  -  1.  3  *  EMACH/  2.  i 

AVELECK)  =  SORTF  C  GAMAE (K)  *  RGE  CK )  *  TE ) 

EMACH  =  SORTF  C EMACH) 

UECO  =  EMACH  *  AVELECK) 

RnOE(K)  =  P  /  (RGE  << )  *  TE  ) 

TNLPECO  =  CPE  CK )  *  TTECK)  -CPECNKP)  *  TTECNKPJ 
EM AS  CK)  =  RHOE(K)  »UECK> 

Ev.CMCK)  =  EMAS(K)  *U£(K>  ■* 

z XT  A  *  PI  * C  REKAX ( K )  +REMAXCK)  -  REMlNCK)  *R£MlN«>) 

SEKOM  a  EMOMCK)  *  EXTA  tSEMOM 
SEMAS  a  E.v.ASCK)*EXTA  +SEMAS 

SEnlP  =  EMASCK)  *TNLP£CK)  wEXTA  t  SENLP  +  .00001 

UD I  V  a  UE  c  n 
TDiV  »  TTEC 1 ) 
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RHOOV  *  RHOEd) 

SHOltf  s  RHOE(l)  *  WEtJJ  #UE(l) 

60 1  EANOZ  «  EANOZ  ♦  EXTA 

End  exit  calculation 

UMAX  *  UE ( 1 ) 

UM1N  »  UE(NKP) 

AVELC  «  AVELE(l) 

URAT  a  (UMAX  +UMIN>/2.0 

UCONV  *  URAT 

QM  a  BMINT 

SMOLD  «  0«0 

WADI  «  RAQC  /  20. 

INITIAL  CONDITIONS  LOOP 

RAO (11  -  0. 

DO  610  KJP*  I .99 
KJ  «  KJP  +1 

RAD  (KJ  )  *  RAD  (KJP)  +  'DR  AD  I 
610  EV I S (KJ)  *0. 

EVISU  )«0. 

XOLD  *  0. 

OUCNV  *  C.O 
NSC  «  0 

END  INITIAL  CONDITIONS  LOOP 

FREQUENCY  NEST  INITIALIZATION  LOOP 

DSPLd  )  *  0. 

Dflll)  «  0.  m 

FREQ I ( 1  I  »  125. 

NOCT  *  2 A 

FREQ J  *  1.259921 

DO  560  JFREQ  *  2. NOCT 

FREQ l <  JFREO )  •  FRE01 (JFREQ  -  1)  #  FREOJ 
DB (JFREQ I  a  0. 

560  DSPL  (JFREQ)  *  0* 

END  FREQUENCY  I Nl T 1 AL 1ZAT1 ON 


DSL  «  0. 

BMRAT  =  BMOR  *RADC/BM I NT 

NMBMa  3.  *  LOGF  (BMRAT  )  /  0.69315 

BM  LOOP 

DO  540  1=  l.NMSM 

CONVECTIVE  VELOCITY  CALCULATION 

UCONN  *  (UMAX  UMIN  )  /  2.0  +  DUCNV*  (BM-EMOLD) 
CMACH  a  UCONN  /  AVELC 
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SPREADING  COEFFICIENT  CALCULATION 

DXBM  =  (CB1  *  (1.0  +  C62* (  1.0  -  EXPF  {-  BM/ (0.2*RA0C ) ) )* 

1CMACH)*  (UMAX  +  UK!Nn/<UV.AX  -  UMTN) 

IF  (I-  1)  100.591.592 

591  DXBMC  =  DXBM 

592  X*  XOLD  + ( DXBM  DXSMO  )*(6K  -BMOLD)  /  2. 

BH  x  BM  /SORTF  (PRADD 

XRA  =  X  /  (XOLD  +  0.0001  ) 

UCONR  a  UCONN/UCONV 
DLNUC=LOGF (UCONR 1/LCGF ( XRA  » 

DUCNV  «  (UCONN  -  UCONV  )/  (X  -  XOLD) 

UCONV  «  UCONN 
I  SET  1  x  BM/(RADC*0. 15) 

NMDOJ  s  2  ** l SET  1 

IF  (NMDOJ  -  8)  593.594.594 
594  NMDOJ  »  4 

593  DOJ  «  NM.DOJ 

DERAD  «  0RAD1  *  DOJ 

NMJMX  x  (RADC  ♦  2.5»BM  )  /  DRADI  +3. 

TMIT  «  O. 

TMAST  *  0. 

TNTPT  «  C. 

FNSS  =  1. 

TYPE  705. X.AVELEd  )  .UE(1  ).BM,  TNLPE  ( 1  >  .  NMBM  .  NMJKX . NMDOJ 

RADIUS  LOOP 

DO  510  jxi,  NMJMX, NMDOJ 
DLI  *  CDIM 
RAD J*  QAD(J) 

END  COEFFICIENT  CALCULATION 

FLUX  CALCULATION  LOOP 

CL  =  RADJ/  REMAX ( 1  ) 

BM.T  IL  =  6H  /REMAX  (  1  ) 

SJ  «  -1. 

CALL  FLUX  (SHK.SJ) 
cMTIL  *  BM  /  REMAX  (1) 

SJ  *  0. 

call  flux  (SF.SJ) 

IF  (SF-0.0C1)  522.522,523 
523  TM  =  EMOM ( 1 )*SF 

M  s  EM AS ( l  >* TNLPE  t 1  >  *SHK 
SHEAR  =  EM0K<1>«SJ  /DXBM 

velocity  iteration 

RG  3  1700. 

U3  =  500. 

■  U2  =  0* 

UI  =  1 00000. 

GO  TO  97 

9 1  AP  3  TM  /2.0  +  CP  4  P/RG 

UN  s  H/(  2.0  *  AP)  +  SGRTF  <  H  +H/<4,0  *  AP*  AP  >  +  (  T«  #  CPE(N 
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1<P)*TT£(N<P>/AP  >) 

UD  «  UN  -  U3 

IF  (  ABSF  (UO)  -  1 •  )  96*  96.  92 

92  U2  «  <U2*U3  -UI+UN  >/(U2  +U3  -Ul-UN) 

95  U1  «  U3 
U3  *  U2 
U2  *  UN 
GOTO  98 

97  RHO  «  TM  / CU3*U3  > 

GO  TO  99 

98  RHO  *  TM  /(UN*UN> 

99  T«  P/( SG  *  RHO  ) 

82  PG  *  RGE ( I ) 

CP  ■  CPE  ( 1 ) 

GO  TO  91 

local  properties 

522  TMASI  »  0* 

TMI  ■  O. 

TNTPI  *  0. 

RAOJR  *  RADJ/RAOC 
SHERR  »  0. 

URAT  •  0. 

VMACM  *  0. 

GO  TO  505 

96  GAMMA«  CP  /  (CP  -  RG> 

TNTLP  *  H  /(RHO  *  UN) 

.TT  ■  (TNTLP  +  CPE(NKP)  *  TTE(NKP)  >/CP 
AVEl  *  SQRTF  (  GAMMA  *  RG  *  T) 

VMACH*  UN  /AVEL 

EDOY  VISCOSITY  AND  DERIVATIVE  CALCULATION 

EVISN  *  RHO  *  UCONN  *  8M  /(2*0  #D*BM  ) 

OEVIS  *  (EVISN  -  EV I S( J  > > / l X-XOLD )  . 

EVIS(J)  =  EVISN 

RATIO  CALCULATIONS 

RAOJR*RACJ/RAOC 

XRsX/RADC 

SHFRR  *  SHEAR/SHDIV 
RHOR  «  RHO/  RHODV 
TTR  «  TT/TOIV 
URAT  *UN  •  /  UD1V 

BMR*8M  /X 

UCONR  «  UCONN/UD I V 

ATTENUATION  CALCULATIONS 

IF  (RADL  -  RADJ )  I*  1.3 
)  RADP  *  RAOL 
GO  TO  «. 

3  RADP  «  RADJ 

2  EFUNC  *  l.C  -  RADP  *  EANOZ  /  (RADL  *  AREAC ) 
XEFEC  »  CX  *  8M  *  EFUNC  *  CMACH 
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EV1SM  *  evist JH-DEVIS*XEFEC/2.C 
FREOO=FR£0 

FREO  =  CFREQ  *  uCONN  /  X 
FREOUEN  *-Y  NEST  INDEX 

IFREO  3.  *  LOGF  IFREQ  /  125.  I  /. 69315  ♦  1. 

SOUND  *■  Ov#ER  CALCULATIOW 

ATKRK  =FRECl*-FREC*£ V  1  SM*  (4. /3.+  ( GAMMA-1  .  l/PRADL  1  ✓  (RHOE  (NKP  )*AV£LE  (N 
1KP)#*3. ) 

ATTN  a  EXPF  (-  ATKR<  *XEFSC ) 

IF(l-l)  100,467,49* 

487  IF  IJ-l)  100,493,494 

493  IFREO  «  IFREQ 
GO  TO  5 

494  FREQ  a  IFREO  -  IFREQ 

IF  (  ABSF(FRED)  -  0.05)  6,6,5 

6  GQ  TO  (7,7.7, 7,6 ) .NSC 

7  SC  «  1. 

NSC  a  NSC  1 

DOSPL  <=  OSPFCSC,  RHOE  (NKP > , AVELE(NKP)  ) 

DSL  =  DSL  +  DDSPL 
GO  TO  521 

8  SC  a  0,5 

DDSPL  *  DSPFlSC,  RHOE ( NKP ) . AVELE ( NKP) I 
DSL  a  OSL  -  0.5  *DSPLO 
DSPLO=  -  0.5  *DDSPL 
DSL  «  OSL  ♦  OOSPL 
GO  TO  521 

5  DSPL  (IFREO)  a  OSPL  (IFREO)  +  OSL 

PUNCH  710,  X«.  RAD JR,  1 FREO. NSC »D5PL (IFREO) .FREQO 
SC  ■  1. 

DDSPL  =  DSPF (  SC« RHOE (NKP),  AVELE (NKP 5) 

OSPLO  *  DDSPL 
DSL  «  DDSPL 
NSC  «  2 
IFREO  =  IFREQ 

MASS  MOMENTUM  AND  ENERGY  INTE  GRAT IONS 

521  TMASI  a  TPI  *  RAD J*FNSS  *  RHO  *  UN 
TMI  *  TPI  #  RAOJ*FNSS*TM 

TNTP I  s  TP!  «  RADJ*FNSS  »  RHO  *  UN  *  TNTLP 

TMAST  a  TMAST  +  TMASI 

TM1T  «  TMIT  +  TMI 

TNTPT  a  TNTPT  +  TNTPI 

IF  (J-l ) 100,483,464 

483  UMAX  «  UN 
AVELC  *  AVEL 
FNSS  »  2. 

484  FNSS  a  6,-  FNSS 
GO  TO  510 

TERMINATION  OF  RADIUS  LOOP 

>  505  RADL  «  RADJ 
XOLO  ■  X 
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I 

II 
■31 
3 
3 


TMIT  a  ITM1T  -  TM I /2»  »*t>£RAD  /3« 

TMAST  a  tTMAST  -  TMASI/2*  >*DERAD/3« 

TNTPT  *  (TNTPT  -  TNTPI/2.  )-*DERAD  /  3. 

.C 

C  INTEGRAL  CHECKS 

c 

FEMOM  *  (SEMOM  -  TWIT  J  /SEMOM 
FEMAS  *  CSEMAS  -  TMAST I  /  SEMAS 
ENRG  *  CSENLP  -  TNTPT  J  /SENLP 

Punch  71  1  .XR.RAOJR.URAT.ShEHR.SNK.  VMACH.SF.SJ.DLi.UCONR.ttMR.DXaM.T 
1TR.NS.N* 

PUNCH  7 1 2 » XR « FEMCM * ENRG  * S J  .  (KJP.rRECl  t*jP).DSPL(«JPl.  <UP  al.NO 
1CT1 

GO  TO  520 

510  PUNCH  7ll.-XR.RADJR.URAT.SHERR.5HK»VMACH.SF‘.SJ.DLl»UCONO.oHR.DXdl'UT 
1TR.NS.NK 
C 

C  END  RADIUS  LOOP 

C 

520  SMOLO  •  8M 

QX8MO  a  DXBM  J 

540  SM  =  1.2S9921*6M 

DSPL t IEREO)  ■  DSPL<lFRtO>  ♦  DSL 
DO  530  KJ  -  l.NOCT 

ORAT  a  DSPL<N-»J  /O.CGOCCOOOOOOOl  -+l  • 

530  D8  <KJ  I  =  10.  *LOGF (DRAT )  *  0.434 

PUNCH  713.  (K  .-PREOI  <<  >.DSPL<  «>»DBC  K>.K  *1.NCCT» 

c 

C  ENO  X  LOOP 

C 

803  IF  (SENSE  SWITCH  21  801.602 

801  GO  TO  571 

602  GO  TO  572 

C 

100  PAUSE 

GO  TO  803 
END 


n 

J 


li 
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subroutine  flux  calculation  procedure 


suanouTiNE  flux  t  psi«  tauj 

COMMON  NX. NS. PI «DL  I « CL «SMT IL»  PNS. XSO.EXPOC.TPSI 

IF  (TAU)  205.  201.201 
205  NSGO  *  2 
GO  TO  202 

201  NSGO  *  I 
TAU  *  0. 

202  XSQ  *  1.0/(8MTIL  *8MTIL1 
EXPOC  =  < 1 .-Cti*<l.-CL>*XSO 
PNS  *  4-. 

PSIP  *  0. 

PS!  *  0. 

TPSI  *  0. 

TTAu  *  0. 

DL  *  0  • 

IF  (CL-1. 1203.204,204 
204  VL*CL-1. 

NS  ■  1.  ✓  OLI  4  0.99999 
VN*NS 


OL  *  1 .✓(  VN  ) 

GO  TO  206 

203  IF (EXPOC  —20* 1222.221 .221 

221  PSIP  *  1. 

GO  TO  208 

222  PSIP*  1.  -  EXPF (-EXPOC) 
208  NS =CL/OLl+. 99999 

IF (NS)  231.231.224 


224  VN*NS 

DL  *  CL  /(  VN) 

VL  =  1.  -  CL 
PSI P=  P1*VL*< 1 .-PS1P)«0 


66666*XS^&- 


PSIP 


206  NS*2«NS 

DO  212  IS*  1.  NS 
NX  *  1  S 

233  VL  =  VL  4  DL 

TARG  *  (CL/VL  +  VL/CL  -  l./CVL*CLl>  /2. 

PHI  =  AT ANF ( SQRTF ( l «  -  TARG4TARG  1/TARG) 

IF  (PHI  )  214.215*215 

214  PHI=  PI  4  PHI 

215  EXPAR  =  VL*VL*XS0 

1 F (25. —EXPAR 1  231.231.216 

216  PSI  =  PNS*EXPF <-EXPAR)*PHI*VL 
GO  TO  (217. 218). NSGO 

217  TAU  =  PNS*EXPF (-EXPAR )*(l.+2.*EXPAR  -EXPF ( EXPAR/2. 1 1 *S I NF (PH  1 ) 

218  TPS l  =  TPS 1  +  PSI 
TTAU  =  TTAU  4  TAU 
PNS  *  6.  -  PNS 

212  CONTINUE 

231  PSI  =  PSIP  4  <TPSI*2.-PS I  )*XSQ*DL/(3.*PI 1 
GO  TO  <  219.2201  .  NSGO 

219  TAU  =  ABSF ( ( TT AU—  TAU/2 . 1 *DL/ (3.*PI *8MTIL 1 ) 

220  RETURN 
END 
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SHORT  RUN  SUBROUTINE 
SUBROUTINE  FLUX  (  PSI.TAU1 
COMMON  NK.NS.P1*  CU.BMTIL 
220  FORMAT  15F10.6I 
NS  •  0 
NK  «  0 
RAO I  •  0* 

8 MO  «  0« 

IF  <TAU1  213.206*206 
213  READ  220 «  XR.RADUR*  SF.SJ.SHK 
RA01  *  RADJR  -  CU 

IF  <  ABSF  fRADl I  *  0«0|  >  209.203.202 
206  BMO  «  XR  -  BMTlL 

IF ( ABSF (0MO>  —  0*0051  201 .201. 202 

202  PAUSE 

TYPE  220.RA01 .BMO. XR. RADJR »CU 
PS  I  «  1.0  • 

TAU  *  0.0 
GO  TO  210 
201  PSI  -  SF 
TAU«  SJ 
GO  TO  210 
2 05  PSI  ■  SHK 
TAU  ■  l  • 

210  RETURN 
END 
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MAUI  PROGRAM  LISTING 


COMPRESSIBLE  FLOW  MIXING  PROFILE  AND  ACOUSTICS  PROGRAM  -  VIDYA-63-4 


SUBROUTINE  ORDERS <N * V, 1 PH) 
DIMENSION  V(lOO>«IPH(|OOI 
LN  *  N 
DO  I  l»l«LN 

I  I phi i )  *  I 

2  K  >  0 

00  4  J*2*LN 
IF  mj-U-VIJI)  3*4.4 

3  VT  «  VIJ-ll 
VIJ-1 »  *  VIJI 
V<J)  «  VT 

IJ  •  IPHIJ-II 
IPH(J-I)  •  IPHIJI 
IPHIJI  »  IT 
K  *  1 

4  CONTINUE 

IF  CO  6*6*2 
6  RETURN 
END 
LIST 

label 


MATRIX  INDEX  AND  MODULAR  OVERRIDE  SUBROUTINE 
SUBROUTINE  MOOSCMFRSTI" 

COMMON  P|«PI2«PRADL*PRASQ*FRE0J. IMH.BM»XE.NEST.NUM*LEAV«PHI .alpo.r 
1A.S1G.C01M* ISAP«BMINT*RADC*SF,SH«SJR.SJC.M00SG.NM00M.M00P 
MODP* (MOOSv>  I )/25 
MFRST  «  I 

>F  (MODP-INMDOM-l  1)216*218*217 
217  NMD0M*NMD0M#2 
216  MFRST«MFRST  ♦  MOOP 
MODX«MOOSG-23*MOOP 
RETURN 
END 
LIST 


*  LABEL 

*  LIST 

*  XEQ 

*  CAROS  COLUMN 

*  LIST 

DIMENSION  XE150) «  ALP0(50)«  LEAVC50I.  NUM(50)*  KN 

1  (50 )  *  1  SAP (50 )«X1TA(S0) . GMW ( 50 ) *  PTE ( 50 ) . TTE (50 ) .ENTLP « 50 ) *PGCE (50 1 • 
2CPE (50  )  «GAMK ( 30 ) «RHOE(50 ) «EM (50  >  *UE (50  > * ACK ( 50 J  * EMAS ( 50 ) *EH( 50 ) «SF 
3(50  |*SH(SO)»  S JR (SO) *SJC(50) « GK ( 50 > ,GKON (50 > . SPL ( 25 > *08 (25 ) *FREQI 12 
45 )  *  I SSH (20 ) *NRE  J (  ),kAFP(25).1SPL(85)«OSPS(85)*ETA(2.50)*UN(20«20 
5)*  IFEX(20.20) *UZ (20*201 . DXBMO 1 20 * 20  )*UCONV (20*20 )*EV IS (20*20  )*OSFL 
5(2Q»20)«X0LD(20*20J *UN" (  «  «  >«UZK  «  «  )»UN2(  *  •  )>UZ2< 

7  *  .  )  *OSPL (25*25*25 )  • 

IA*SIG.CDIM* IS6P*BHINT*RADC«SF.SM.SJR*SvtC*M0DSG*NMD0M«M0DP 
COMMON  P I  «P 12, PRADL «PR ASQ «FREOU »  JHH*BM* XE« NEST. NUM«LEAV. PH| « ALPO.R 


.!  FOAM  AT  (6110* 2F 10.4) 

2  FORMAT  (2FlO.A» 

3  FORMAT  l2F|0tA.31 JO) 

4  FORMAT (  80H1 

isle  flow  mixing  profile  program// /osh 


COMPRESS! 


3  FORMAT (18H 
1 

6  FORMAT (1  OHO 

I13.I2H 

2 

8  FORMAT  I  1 3H0 

I  ( 12t2H)»F7«4» 
22.2H ) *  I 3.7H 

9  FORMAT C23H0 


INPUT  DATA///) 

0PFB«5.17H  IPEU.3.14H 

0P1 1 1 *3 . I 3H  IPEII'A//) 

NCI-I3.J3H  NCF*I3*12H  LPHU 

I SYM* 12*1 AM  NEST-I3.13M  INOX-13  //19H 

RA0C-F6.3.9H  CD1M«F6.3.8H 

XE<  12.2HUF8.2t-  7H  ALPO 

7H  LEAVI  I2.2HUI3.  7H  NUMI 12 .2H I ■ I 3.7M  KN( ! 
I5>A»<  I2.2HU13) 

0ALP(I2«1H, 1 2 * 2H ) »F7»A »7H  RA( 12. 1M, 


I  12,.2H)»F7.4.7H  DALPl  J2.1H.I2.2HUF7.4.TH 

10  FORMAT  I 76H1 

1  CASE  NUMBER  14* 

11  FORMAT (  61  HI 

1N01TI0NS///  93M 
2CITV  MACH  MOMENTUM 

3  TEMPERATURE 

4  FLUX//  9SM 

5/SEC  LB/SO  FT 

6  I2.1IH  0PF8.2.7H 

7E9.3.  6H  1PE9.3)) 

12  FORMAT (11 OHO  RADIUS 

1MENTUM  ENTHALPY  RADIAL 

2 

3  SHEAR  SHEAR 


RA( I2« 1H« 12«2H)bF7«4) 
OPERATING  CONDITIONS  FOR 


CONTOUR  STATIC 

ENTHALPY/  93H 

NUMBER  Fl 

DEG  RANK I ME 
LB/SO  FT//(19M 
F8.2.  6H  F6.3.  5H 


EXIT  CO 
VELO 

FLUX 


ANGLE  VELOCITY 
TANGENTIAL  total 
flux 

SHEAR  /ll OH 


ARADIANS  FT/SEC 
5  LB/SQ  FT  LB/SO 

13  FORMAT (  9H 
I5E12.3/I 

14  FORMAT  1 57H 

1  =  1  PE  I  I  «4//l 

15  FORMAT (6 12  ) 

16  FORMAT  <5Fl 0*4  ) 

17  FORMAT  (69H0  EXIT  AREA  OF  C 

IONTOUR  NUMBER  12. 2H  ■  F10.4  ) 

18  FORMAT  (94H0  CONTOUR  TOTAL  T 

1EMPERATURE  TOTAL  PRESSURE/  91H 

2  DEGREES  RANKINE  LB 

3/SO  FT//UBH  12.  3 OH 

4  FB.2.  24H  F8.2/I1 

19  FORMAT  <  89H1  CONTOUR  MOLECULAR  WEIGHT 

1  LAMBDA  ONE  LAMBDA  TWO  /89H 

2  FT2/S2-DEGR  IFT/S-DEGR 1SQ//20X13* IB 

3XF6. 3* 12X1  PE 10.3.10XE10.3  ) 

20  FORMAT  <  9XF 1 0*  3.5XFS  *3 .5XF6.A . 4XF6.4 . IPSE  12*3  ) 

23  FORMAT  <  40H  SOUND  POWER  SPECTRUM  ///29X  H  F 

1REQUENCY  FREQUENCY  SOUND  POWER  S  X 

2  H  INDEX  RB  REF  10-13“  /////// 

3  13  1PC10.3  QPF10.B  I 

300  FORMAT  <40H|  PSUEUDO  AXIAL  LOCATION-C  SUB  M«F»*3.25M  FT  C 

1ASE  NUMBER-  13  //I 

302  FORMAT  (1H  ) 


FT/SEC  DEG  R 
LB/SO  FT///) 

1PE10.3.5H 


LB/SQ  FT 


LB/ FT  SEC 


TEMP  MC 
/UOH 

FLUX 

FT 

LB/SO  FT 


0PF5.3.5H 


F6.I.4H 


F6.1.1P 


TOTAL 


IFT/S-DEGR 1S0//20XI3. IB 
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303  FORMAT (7BHt  SOUND  POWER 

1  INDEX  A NO  MASS- FLUX//  113»  -  PSUEOO  AXIAL  LOCATION  S 

SOUND  POWER  LEVEL  INDEX  MASS  FLUX  VELOCITY  TO  THE  8TH 

20 UNO  POWER  LEVEL  1N0EX  FREQUENCY  INOEX  VELOCITY  TO  THE  jBTH 

3—  INTEGRAL/1 06H  FEET  I FT  LB/SO 

4  SEC ISO  FT**10/SEC*«8///J 

304  FORMAT  I  76MI  COMPRESSIBLE  FLO 

IW  MIXING  PROFILE  PROORAM//  - 

234h  CASE  NUMBER— 12/ 

3  4 OH  PSUEOO  AXIAL  LOCAT 

4I0N,«  F9.5.3H  FT//  TOM  RADIUS 

5  ANGLE  TOTAL  PRESSURE/  BOH 

6  FEET  RAO IANS  INCHES  HG  GAUGE! 

303  FORMAT (23N  F9«3*  X3H  F6.3.  19H 

1  F7.3! 

307  F0RMAT(213.3F10*5.4F10*2  ! 

309  FORMATION  XI !2.2H)»F6*3«  9H  SPL I < I2.2H)* IPEl 1 *4 .BM  FM(I2 

1 .2H)*lPcIl*4.  OH  UOI t!2«2H»a|PEl!«4//) 

329  FORMAT  < 10X*F10*2.F10*2I 

330  FORMAT  (110! 

DSPF (  A  .  B I  ■  A  *  B  #SORTF(TU2RO*TU2RO*TUSRO*TAU! 

PI *3. 141 5927 
PI £*6*283 1853 
NT2  «  2 


INPUT  OATA  PE AD IN  AND  PRINTOUT 


READ  INPUT  TAPE  NTS  * 15.NT3.NT4 .NT5.NT6.NT7. NTS 
REAO  INPUT  TAPE  NT2. 1 .NCI «NCF. INDX.NEST.LPMI . ISYM.RADC. COIN 
READ  INPUT  TAPE  NT 2 *3*  t XE (X ) * ALP0 ( K 1  *ISAP(K >*LEAV(K) .NuM(K),KN«! 
1*K»1«NEST! 

FLPH1 «LPHI 

WRITE  OUTPUT  TAPE  NT3.4  - 

WRJ.E  OUTPUT  TAPE  NT3.5 .NC l .NCF .LPHI « ISYM.NEST. 1NDX.RADC.COIM 
S£MAS*0*  .• 

SEMOM-O* 

SENLP« 0.00001 
AR£AC«P1 *RADC«RAOC 
EAN02*0* 

DO  319  tC*2«NEST 
NUMK*  NUMOO 

READ  INPUT  TAPE  NT2.2* (DALP <N«K 1 ,RA<N.tC»«N«l «NUMX| 

1 SAPE* I SAP (K ! 

FLEAV  ■  LEAVIKl 

WRITE  OUTPUT  TAPE  NT3.8.K.XE <K ) «  K.ALPOdCl.K.LEAVdCl.K.NUMHC! 
l.K.KN(K).K<ISAP(X! 

NUMK*NUM(K> 

WRITE  OUTPUT  TAPE  NTS. 9. (N.K.DALPIN.X 1 .N.K.RAfN.K) «N»1 .NUMKI 


EXIT  AREA  CALCULATION 


GO  TO  1315.316.3171. ISAPE 
315  XITACKI  *  .25MRAU.KI  4  RAIN U 
DO  313  N  «  2.  NUMX 
313  X1TAOC  )*X I TA (K  »+  .2S*<RA(N.K> 
X  IT  A  OO  »  XITA«Kl#FLEAV 
GO  TO  319 


♦  RA(NUMK.K))«CALPU.K> 


RAIN-1.K1 >»DALP(N.«! 
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V  .V.V.V.V. 


* 


316  XlTAK)«PI*flA(l.K)»Rft(l  «K> 

GO  TO  319 

317  £1  •  RA(3.K)  *RA(3«K)  +flAU>K)*RAI|«K)  *2(*AAU<X|*OA(3«K}«COSFIM 
ILP(l.Kl) 

S2  -  RA<1«K)*>RA<1  .K)  +PA(2.K)«RA(2.K)~2.*RA<1 *K)  *«A(Z.K)*COSFCDAU 
1PI2.K 1 1 
X1TAIK }*Sl*S2 

WRITE  OUTPUT  TAPE  NT3, 17.K*XltA(K) 

319  EANOZ-EANOZ  ♦  XITAIK  > 

c 

c  EXIT  CONDITIONS 

'  00*999  NC  •  NC1.NCF 

READ  INPUT  TAPE  NT2.307.OCTO.NOCT.XINT.XOR.PRADL.FREQB.CMW.PS.CP 
READ  INPUT  1  APE  NT2  »  1 6  *(PTE  (K )  «GMW (K 1  .TTE(K) «EXA(K« 1  l.ETA  <K»2  1  »•“’* 
II. NEST  1 

.WRITE  OUTPUT  TAPE  NT3, IO.NC.OCTD.NOCT.XlNTvXOP.PRAOU.FREOa.PS.CHW. 
1CP 

DO  119  K-l.NEST 

119  WRITE  OUTPUT  TAPE  NT3. 18. K rPTEIK 1 .TTE «K ) 

IF (ETAf 1.1  1)120.120*121 

121  DO  120  K-  1 .NEST 

WRITE  OUTPUT  TAPE  NT3. 19*K.GMW(K) *ETA «*t 1 .ETAtK«2t 

120  CONTINUE 

ENTLO  •  <ETA<  l « 1 l+ETAI 1 .2  )/2«*TTE ( 1 ) )*TT£< 1 ) 

LBPAS- I 

00  113  K-l.NEST  . 

IFtETAU.l  nill.111.112 

112  ENTLPIKl  -  (ETAlK.l  >+ETA<K.2)/2.*TTEtKn»TTEtK» 

PGCE  CK  1  -49690 ./GMWIK) 

PGC  ■  PGCE IK 1 

IF  ( PTE (K l-PS 1  321.321.322  # 

322  TEK-TTE(K>*(PS/PTE1K)  1**0.206  #<» 

IIS  CPE (K)-ETA(K.l 1+ETA1K .2 )*TEK 
CP  »  CPE (K 1 

STeNl -1ETA  IK.  1 1+ETA1K  .2  )/2*i»TEK  )#TEK 

LBPAS-2 

GO  TO  378 

1 1  1  PGC  *  49690./  CMW 
PGCE I l 1 -PGC 
ENTCPIK1-  TTEIK1«CP 
CPM.CP 

ENTLO  «  CP»TTEM1  * 

378  GAM  -  1 •/ 1 1 PGC/CP1 
GAMKIKl-GAM 

ACH2  «  2.  *  I  (PTE  <K  )/PS)*«(l.  -l./GAHi>U  1/1GAM  -  1.) 

GO  TO  (407.408 1.L8PAS 

407  TEK  »  TTE (K  )  /  (1.+  ( GAM-l • )*ACN2/2. 1 

408  PKCE(K)«PS/(PGC*TEK1 
SV2  -  GAM  *  PGC  4  TEK 
GO  TO  (l 16.381 l.LBPAS 

381  STEN2=ENTLP(K)/( 1 .♦ACH2/2.*GAM/RH0E(K)«PS 
DSTEN.STEN2-STENI 
TeK»OSTEN/CPE (K i+tek 
IF  ( ABSF ( OSTEN— 1 0» i  116.116*119 
321  CPE(K)  ■  ETAlK.l  1-.ETA  (K.2  )*TTE  (K) 

CP  «  CPE (K 1 
LBPAS-  1 
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CO  TO  378 

116  EMIKI  a  ACH2  *  GAM  «  PS 

UEIKI  >  SQRTFI ACH2  #  SV2  I 
SVEIK)  a  SQRTFISV2I 
ACHIK)  a  UEIK1/SVEIK1 
EMASIK)  «  RHOE IK  1 *UE ( K ) 

.ENTLPtK)*ENTLRtKJ  -*ENTLO 
113  EHIK1  a  EMASIK >«£NTLP(K I 
CAM  a  SAMKIll 
PGCaPGCE(l) 

A5Pa  PS  #GAM  #GAM  «RGC  «$ORTF|CAN  #PGC  #  TE(1> 
UDJ  V  *  UEUNOXI 
T01-V  •  TTEllNOXl 
PHOOV»RHOE ( INOX } 

SHOIV  *  RHOEIINDX)  *  UEtlNDXI  ♦  UCtINOXI 
WRITE  OUTPUT  TAPE  NT3.lt .  ( K » T EK  •  UEIKI 

*1  )«<■  1  •NESTI 
WRITE  OUTPUT  NT3  «IA«A5R 
00  122  K>2,NEST 
REMAS»RHO€ IK  1  # IU€ (K  )-UE 111! 

semas«semas  +REMAS*XITA|K> 

SEMOM»  SEMOM+REMAS*XI T  A IK )  #  IUE  IK )-U£  < 1 >  > 
SENLP*SENLP+REMAS«ENTLP  IK  >#X  1TA  IK  » 

END  EXIT  CALCULATION 

FREQUENCY  NEST  INITIALIZATION  LOOP 


IF  INOCTI  313.513.514 

513  NOCT  a  10 

514  IFIOCTOI  515.515.516 

515  OCTO  a  1 

516  SPLIDaO 
DB I  1)  •  O 
FREQ! 1 1  la  FREOS 
NOCFO  a  10#  NOCT 
OCT  a  NOCT 

NOCT  a  OCT#OCTO 

FREOJaEXPF 10. 693 1 5/OCTOI 

00  561  JFRE0«2 . NOCT 

FPE01 I JFREQ I aFREQ 1  I JFREQ-l J#FREQJ 

DBIJFREOlaO, 

561  SPLIUFREOI  a  0.0 

Not  Operational _ 

SOUND  POWER  INITIALIZATION 

00  151  N« 1 .NOCT 
DO  151  K*2 .NMBM 
00  151  M* 1 . M I  XT 

151  OSPLIN.K.M)  a  -l.  _ 


VELOCITY  INITIALIZATION  LOOP 


C 

00  511  Mai, 20 
00  511  I ■ 1 .40 
SI  1  UNI  I .Ml  a  o. 

c 


»*TEU» 


•  ACHIKl.EMIK > .EH  IK 
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0 


m 

a 

.3 


■,a 

a 


41 


H 


Ti 

7 

a 


* 

'1 


-.1 

i 

_* 

t 

■j 

7 


t  INI T  1*1,  CONDITIONS 

t 

ft*U.  *  44  /<7»  LPMll  ♦  1 

*riN?»o 
I  N«€J*0 
4f TF*» • 

?eh*»i» 

ISC?  •  0 
NMPCM*  I 
I  SGI  I  •  0 
Oft  • 
moACK-  •  0 
NHJHX  •  11 

pru-o. 

^  mat  «  l 

’  SET!  •  0. 

rweo  •  o.oi 
xolo  •  o. 
oucnv  •  o* 

ST HP  a  0. 

X«0. 

RAOL-RAOC 

C  1N0EX  REGISTER  DETERMINATION 

C 

|E  (XI  NT I  94.94.9S 
-94  XlNT  «  COIN  /  4.  *  0.01 

95  |F(XOR>  517.517. SIS 

517  xon  «  io. 

518  BKRAT  «  XOR  *  RAOC  /  XlNT 

M<8M  «  OCTO  *  LOGF  (0HRATI  /O. 69315  ♦  2» 

BMINT  «  0.075  •  XIMT  - 

BH  «  C. 

0S1G2.RADC  /  10. 

C 

c  B  M  LOOP 

C 

DO  185  KA-1.NM8M 
I HP A 5  ■  1 

|F  (KA-ll  99.99.100 
1 00  MI1PAS  •  1 
C.O  TO  118 
99  MtlPAS»2 
t | M  WACO  «  2. *BM 

UMC.-DS1GZ4SET2 
00  7 27  I  I • I .NOCFO 
1C  t  O'.PM  I  I  **°» 

U".4C>I*  *  DSIG/8.0 


TAPE  NT3.200.8H.NC 
TAPE  NT6.300.BH.NC 


T  MAS  I  «  0. 

1 H I T  -  0. 
1NIPT  «  O. 

(  T  (  Nt  *  I 
Wi»  I  -O. 

Ml  «*0 

Wl<  4  1 4  out  put 
*Hitr  output 
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WRITE  OUTPUT  TAPE  NT7.300.BM.NC 
l_»Nt3«S 

LINE6*0 

LIXE7-8 

WRITE  OUTPUT  TAPE  NT3 
WRITE  OUTPUT  TAPE  NT6.12 
WRITE  OUTPUT  TAPE  NT7. 

FLOW  FIELD  MATRIX 


990  SIG*0.001*DSIG 
I BP AS  ■  I BP AS* l 
BUG  a  SIG 

RADIUS  LOOP 

DO  286  M*l,NMJMX 

OPHI  *  PI2/FL0ATF  <LPHI#IS» 

DPMI  «  DPHI  *SET2 

ISSY  a  IP|L«(M>1)  ♦  t  -  ISVM' 

ISSM(M)  a  ISSY 

DPHIPaOPHI/0«O 

PHI aO« 


ANGLE  LOOP 


DO  152  I  * l « I SSY 
IM»0 
MM*0 
MAT  a  1 
IFEX  a  0 
NVPAS  •  1 
DERAO  a  OSIG 
OEPMI  a  OPH1 


NO  MESH  SET  ON  NMDON 
NMOOM  a  8 


BYPASS  REGISTER  -FLUX  SUMMATION 

GO  TO  053. 205*5051  .IBPAS 
505  GO  TO  (204.288I.M8PAS 

FLUX  LOOP 

153  SET  a  0. 

SHEAR  a  0* 

SHEAC  ■  0* 

SHEAT  a  0» 

TM  a  EM< I ) 

H  «  EH  O  » 

CALL  F'_MA  <NC»  i .m.kai 
DO  971 K«£  «n£ST 
K«2+NEST-IK 

IF  (SF(K»  -  .999J  All. All. *12 
*12  KNK  a  | 


138. 
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CO  TO  413 
411  KNK  «  KN1K1 

413  TM  »  TM  ♦  (EM1K1  -  EM1KNK  »)*SF CIO 
Ha  H  +  (EH(K)  -  EH < KNK 11*  SHCK1 
GK (K )«EMAS <K>  *  SH<K» 


c 

c 

c 


X  ZERO  SHEAR  AND  ITERATION  BYPASS 


GO  TO  <328a349>aM3PAS  * 

328  SHEAR  •  SHEAR  -MEMIKI  -  EM ( KNK  1 1  * 

SHEAC  a  SHEAC  +  1EM1K1  -  EM1KNK1)  *  SdClKl 
SHERT  «  S0RTF1 SHEAR  *  SHEAR  ♦  SHEAC  *  SHEACI 
349  lF<SF;K)-0a999T9?«97«324 
97  SFT»ai-T+SF<Ki 


If*  1 

324  lF(SFT-OaOO!)272aS23«S23 
C  • 

c  POTENTIAL  CORE  BYPASS  * 

C 

272  UNU  .M1-UE1K)  ' 

TNTLPaENTLP<«> 

RHO  •  RHOE1KI 
GAM  a  GAMK(K) 

SHEAR  »  Oa  * 

SHEAC  ■  Oa 

TS  •  TTE (K  )/ ( 1  a  + < GAM- 1 a  1 /2 a * ACHM* ACHM 1 
ACHM  «  ACHfKI 
GO  TO  362 

S23  GO  TO  1 102*334  I »MBPAS 


END  BYPASS 


VELOCITY  ITERATION 


334  U3  ■  500 

IFEXII «MI«1 
GO  TO  98 


102  U3  »  ONUaMl 
98  U2  ■  Oa 


U1  *  lOOCOOa 
GO  TO  103 

104  AP»TM/2aO+CPM*PS/POC 

UNdaMl  H/(2a0*AP>  ♦  SORTF 


{  H*H/ 1 4  a 0* AP*AP I 


107 


UD  a  ON(IaM)  -  U3 

(ABSF  tOD- 1 • 1  l 05* 1 05 « 1 OT  * 

02  «  <U2*U3-U1  *0N<  1  «M)/  (02+U3-01— ON11  ♦  “ »  » 


01  a  03 
03  •  02 
02  «  ON t I «M> 

GO  TO  108 

103  RHO  a  TM/ (03  *  U3» 

TS  *  PS/< l700a*RHO) 

GO  TO  109 

108  RHO  a  TM  ✓  <ON(I«M>  •  ONtlaNU 
TS  »  PS/<PGC  *  RHO 1 


♦ 


<tm#entu>  /ap»i 


*-*<W 


■-*v- 
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109  IF(£TA(1 *11}  1 04. 104. 1 1 4 

CP  AW  PGC  LOOP 

114  TOTLK  •  0. 

CPM«0» 

CP«0* 

PGC  •  0*.' 

00  521  NK  -  2»  NEST  * 

GKON  INK)  «  GK  I NK  )  ✓  ( RHO*UN  (  I  «M  )  ) 

PGC  ■  P6C?«NK)#GK0NINK  >4PGC 

CP  ■  GKON  (NK)  •  (ETA INK* t )  4  ETA  INK.  2)*TS  ♦  ETA  (NK.  3)/<TS* 
ITS))  4  CP 

CPMbGK0N(NK)4{£TA(NK.1 )4ETA(NK.2)/2.«TS)  4CPN 
321  TOTOC  ■  TOTLK  4  GKON  (NK) 

.  PGC  ■  PGC  4  (1.  -  TOTLK)  #  PGCE  (I) 

CP  •  CP  4  (1.  -  TOTLK)  *  (ETA(l.l)  4  ETA  ()«?)  •  TS 
CPU  *  CPM  4  (1.  —TOTLK ) • ( ET A ( 1 « 1 )  4  ETA(1.2)*TS) 

GO  TO  104 

LOCAL  PROPERTIES 

103  GAM  «  CP  ✓  (CP-PGC) 

TNTlP  ■  K/(MHO*UNII«N)) 

SV  ■  SQRTF  (GAM  *  PGC  •  TS) 

TT-TNTLP/CP 

ACHM  »  UN(I.M)  /  SV 

382  TPG»  PS*(1.4(GAM-».)/2.«ACMM  ♦  ACHM)  **  (GAM/(GAM-).))-PS 
TP  «  TPG  4  PS 
HG1N  •  TPG/70.727 

WRITE  OUTPUT  TAPE  NT6. 1 3«S IG »PMI  .UN (  I  .M.TS.TM.H* SHEAR. SH6PC  .  SHERT 
IE (SENSE  SWITCH  3)  213.214  • 

213  WRITE  OUTPUT  TAPE  NT6.303 . S IG.PH1 «HG1N 

214  GO  TO  (4)4.413) .M8PAS 
413  U?(I.M)  ■  UN(I.M) 


oouuoo  Uuu 


X  BYPASS  AND  MESH  BYPASS 

/• 

414  60  TO  <S62*17£.206.X07)*  MAT 

INTERROGATION  LOOP 

204  MXPLN  ■  1 
UMAX  «0« 

255  UMIN  »  UN  tl«M) 

MOM  «  RAC0/0SI6  4  1* 

>101  a  RACO  ✓  <  SlG*DPH!)  4-1  • 

MM  IN  *•  MAXOF  (M-MOM.l) 

MMAX  a  MINOF  < M-f MOM » MRA DL  1 
1M1N  *  MAXOF(  t-IOI.'l) 

1MAX  *  MJNOF  ( 14101* 1SSY1 
DO  101  MP  a  MM IN*  MMAX 

11  *  ( ( I M IN  -1 >*(MP-1 > )/(M-I )  41 

12  ■  (  (  1 MAX  -1 )*{MP-I ) >/(M-l )  4  1 
DO  181  IP  *  11*12 
60  TO  (883*8821* MXPLN 

882  UMAX  »  MAX IF  <UZ  ( l P«MP >  «UMAX 1 
UMIN  a  MIN1F(UZ(IP.MP) .UMIN) 

60  TO  881 

883  UMAX  •  MAX1F  I UN( IP»MP ) «UMAX) 

UMIN  a  MIN1F  (UN( IP.MP ) .UMIN) 

881  IF  (NMDOM-B)  333.181*161 
_  This -  section,  now  ,  bvnaaaed 

INNER  MESM  INTERROGATION 

333  1FX  a  ( IF£X0(MP«1P)  1/100  ♦  I 
IND1  «  IFEXO(MP.IP) 

CO  TO  (187.192*193. 194). 1FX 

192  M0DS6  a  1 NO I 
CALL  MOD (MFRST ) 

DO  195  MPP=MFRST»8.NMDOM 
GO  TO  (884.695) .MXPLN 

885  UMAX  a  MAXlFtUZl (MODX. IP. MPP ) .UMAX ) 

UMIN*  MINjrlUZl (MODX* IP .MPP) .UMIN) 

GO  TO  195 

884  UMAX  a  MAX1F (UNI (MODX  « I P. MPP ) .UMAX ) 

UMIN  a  MIN1F(UNI (MODX  « IP. MPP. UMIN) 

195  CONTINUE 
60  TO  181 

193  MODSG  ■  1 NO I 
CALL  MOD (MFRST) 

00  196  lPP«MFRST.a.NMDOM 
GO  TO  (886. 867). MXPLN 
887  UMAX  a  MAX1F(UZ1 (MODX. IPP.MP), UMAX) 

UMIN  >  M1N1F1UZ1 (M0DX.1PP.MP1.UMIN) 

60  TO  196 

886  UMAX  a  MAX1F  (UNI (MODX  , I PP .MP ) .UMAX) 

UMIN*MINIF(UN1 (MODX  .IPP.MP  ) .UMIN) 

196  CONTINUE 

GO  TO  181  * 

194  M0DS6  »  IND1 
CALL  MOO (MFRST) 


I 


00  198  MPPaMrRST.8.NM00K 
00  198  1PP*  mfrst  .  6«  NMOOM 
60  TO  <668*8891 .MXPLN 
889  UMAX  a  MAXIF<UZ2<M0DX. IPP »MPP > .UMAX  I 
UM1N*  M1N1F 1UZ2 (MODX? JPP.MPP )  ,U*INI 
GOTO  198 

886  UMAX  =  MAX1F  <UN2 (MOOX  . IPP.MPPl .UMAX 

UM1N  «  MIN1F  4UN2  < MODXMODX  • ! PP « MPP ) . UM 1 N 1 

198  CONTINUE 
181  CONTINUE 


.  CONVECTIVE  VELOCITY  CALCULATION 

187  GO  TO  <281<2821<MXPLN 

281  UCONN  «  < UMAX  t-  UMlNl  /2.0 
CMACH  a  UCONN/ SVEUJ 

OXBMx  (CBl*(l»4C62*<l  »*"EXPF  <-XOLD/  <  2  »*RAOC  >  I  1«CMACH  1*  (UMAX+UM INI 1/ 
1  4UMAX-UM1N) 

IF  (KA-3)  875*873.873 

873  IF  4X0L01  874.874.875 

874  MXPLN  *  2 
GO  TO  255 

282  uCOn P»  (UMAX  «•  UMINI/2. 

CMAP  =  UCONP/SVE  <11 

OX8M04  1  . M  1  a  ICB1*  (UMAX-tUMlN)  >/(UMAX-OMl N 1 
X  DETERMINATION 


X  «  <OXOM  ■♦OXeMCCl*  M11/2.«8M 
GO  TO  876 

875  XcXOt_D  (  I  «  M  )  ♦  ( DX8M0  <  1  .Ml  40X8X1 

876  TAU-SKfcRTI I .M1/0X8M 
OXQMOU.Ml  a  DXBM 


F2.0*(BM-BMOLD» 


EOOY  VISCOSITY  AND  CONVECTIVE  CALCULATIONS 


XRA  s  X/<XOLD+O.OOQ1 I 
IF  <  KA-2  1  34S.345.346 

345  UCONV  < I .Ml  *  UCONN 

346  UCOVR  »  UCONN  /  UCONV U .Ml 
OLNUC  «  LOGF  <UCOVR)/lOGF<XRAI 
UCONV (I. M|  *  UCONN 

EVISN  ■  RH0  *  UCONN  *  8M  /<2.0*0XBMI 
OEVIS  a  <EVISN  -  EV 1 S< 1 »M  11  /  IX-XOLOl 
EVIS  < I .M  1  a  EVISN 

ATTENUATION  CAi_CULAT  ION 


I F < PAOL  -  SIG1  199.199. 201 
199  RAOP  a  RAOL 
GO  TO  202 

201  RADP  «  S1G 

202  EFUNC  «  1.0  -  RADP*E ANOZ  /  <«ADL  « 
XEFEC  a  CX  *  BM  *  EFuNC  »  CMACH 
EVISM  a  EVISN  +  DEVIL  XEFEC. '2.0 
FPEOO  «  freo 

FREQ  »  CFREO  ♦  UCONN  /  X 


ARE  AC  1 


-li£- 


no  nflft 


AT<RK  *  FREO  *  FREQ  *  EVlSM  «  t A./3.+  {GAM- 1 •  1/PRADL I  /(RMOEIU*SVE 
ATTN*EXPF  *-ATKRK*XEFEC I 


FREQUENCY  NEST  I NDECES 

IFREO  «  OCTO  *  LOGF  {FREQ  /  FREOB1  /  0.69315  +  1.0 
NFREO  *  10  *  LOGF  (FREO/FREGBl/0.69315+1.0 
NFREO  *  10  #  LOGF  {FREO/FREGB  1  /O. 69315  +1.0 

GO  TO  562 

GO  TO  215  ^ 

•  _ This  section  now  bypassed 

MATRIX  SUBDIVISION 

205  IF  {NMOOM  -8)  355055.125 
356  1FUFEXU.M1  -2001  352032.125 
.  352  IF  (lFEX(l.M)l  354054.353 

353  1NREJ  «  INREJ  +  1 

NREJ(INREJ)  •  M0DFUFEX{1  .Ml-l.lOO)  +1 

354  IF  (M-l)  125.125.129 
129  IF  ll-ll  125.125*133 
133  INDI  *  I 

II*  Ul-l)*«M-2»»/{M-n  +  l 
SFSGR  «  SFTIl.Ml-SPTm.M-l) 

I F ( ABSF  4  SF  SGR I  —0  .11  152.152. 131' 

135  INDI  •  1 NO I  +1 

152  SFSGP  *  SFT { 1 «M l-SFT ( 1 — 1 »M) 

IF (ABSF(SFSGP 1-0.1 1125.125. 158 
158  INDI  •  INDI  42 
NVPAS  *  2 
GO  TO  123 

337  GO  TO  <125. 167.168.1691. INDI 

167  1FIN1  *  IFINl+l 
IFEX  1 1 «M l  »  IF1N1 
MODSG  -  IF INI 

SIG  •  S I G+ I SET 1 -7 .5 1 »DS I GP 
MAT  «  2 

DERAO  «  DS1GP  *OEL 
CALL  MODS {MFR5T 1 
'  IF  {NMDOM-81  336.562.562 
336  DO  203  MM*MFRST .8.NMD0M 
GO  TO  153 

172  UNHMODX.I  »MM)*UN(  I  .Ml 
GO  TO  (416.4171.MBPAS 
4  i  7  U2UM00X.  I  .MM)  *  UNIJ.M  1 
416  SIG  *  SIG  ♦  DERAO 
GO  TO  123 
203  CONTINUE 
GO  TO  562 

168  IFIN  *  IF  IN  ♦  l 

-  JFEX< I .Ml  *  JFIN4 100 
MODSG«lFINl 

PHl«PHl  -MSETl-7.51*0PHIP 
DEPHI  *  DPH1P  «OEL 
MAI  •  3 

CALL  MODSIMFRSTl 
IF {NMDOM-81  229.562.562 


-llo- 


oo  u  o  o .  uuu 


229  00  209 ' t M«MFRST  «8*NMD0M 
GO  TO  153 

206  UNI (M00X*IM.M)«UNtI*M> 

GO  TO  (418*419)* MBP AS 

419  UZ 1  ( MODX * IM*M  )  *UN ( 1 *M) 

418  PHI  «  PHI  +  OEPH1 
GO  TO  125 
209  CONTINUE 
GO  TO  S62 
169  IFIN2-1FIN2+1 

,  IFEXd  *M)«IF1N2+200 
M0i)SG»lFlN2 

SIG  «  SIG  ♦  (SET1~7*5)«0SIGP 
PHI  >  PHI  +(SETI-7«5)*DPMIP 
OERAD  *  DS1GP  *DEL 
•  0EPH1  *  OPHIP  #OEL  ' 

MAT*4 

CALL  MODS < MFRST 1 
!F(NM00M-8)  223*562 *562 
223  DO  211  MM«MFRST.8*NMD0M 
00  212  IM«MFRST*8,NMD0R 
GO  70  153 

207  UN2(M0DF(M0DX  «1M*MM)  •  \jn\  1  .Ml 

GO  TO  (501 *5021* MflPAS 

502  UZ2(M00X, IM*MM)-UN<1«»1 
501  PHI  ■  PHI  ♦  DEPHl 
GO  TO  125 
212  CONTINUE 

SIG  »  SIG  ♦  DERAO 
211  CONTINUE 
MAT  -  1 

_ _ GO  TO  362  _ _ _ 


SOUND  POWER  FACTOR 


125  7U2R0»TAU/RHQ*TAU ' 

DSFLU«“)=  PSFLU«M)  +  05PF  <  DERAD  *  DEPHl  | 
GO  TO  (562*344)*  NVPAS . 


RATIO  CALCULATIONS 


562  RR  »  SIG/RADC 
SR1«SH£AR/SHD!V 
SR2  «  SHEAC/SHOIV 
SR3«  SHEAT/SHOIV 
RHOR  «  RHO/RHODV 
URAT  ■  UN ( I « M )/UDI V 
TTRbTT/TDJV 
TR«TSFTDIV 
8 MR  •  8M/X 

•WRITE  OUTPUT  TAPE  NT3 . 20* RR . PH I *RMOR* URAT « TR* TTR.SRl «SR2*SR3 

INTERMEDIATE  8YPAS  REGISTER 

* 

GO  TO  ( IS2.541 «2tS)« IBPAS 
541  GO  TO  (344*152)  .NVPAS 
C 


ft  o  o 


SOUND  POWER  CALCULATION 

>15  FRED  *  NFREO  -  NFREO 
I F  <  ABSF  <  FRED  >  22*  ♦  22G  ♦ 

pR»  FRE0IKA.N1XI*UC©NVI I »MJ/UCONN*X/XOLOt I «Mt 
jFRo  l 0.»LOGF (FREO/FREO I 1 1FREO l >/ .69315  +i# 

I DIF  •  JFR-NFREO 
JF<IOIF»  358*362*358 
>58  TF{OSPL<NFR£Q*KA*MlX+I >  3(2«3(2«>^3 
IF IMIX-20 1363*364* 364 

363  NIX  «  NIX  ♦! 

OSp”«rSQ^«.«.*.-OSPU<NFBE0.«».-.K.*05L*C«^U«'L«n/*» 

364  0SPS<NFRE0»«DSPS<NFREG)+DSL»CA/Pl2*Ft.PrtJ/AJ» 

DSL-O. 

.  GO  TO  1356. 506). LFINL  _  **iiroMN*ATTN 

356  OOSPL  -  DSFLU.M>40XBM#*-7/U.-0LNUC)*UC0NN*ATTN 

226  OOSPL  »  0SFLU«M)*DXBN*4-7/«.-0LNUC»*UC0NN»ATTN 
IF(l-l)  273*273*276 

275  IF<1-1SSYI  276*273*273 
273  D0SPL-0DSPL/2. 

276  DSL  *  DSL+OOSPL 

U8*UN«I.NI**8 

GO  TO  166 

MASS  MOMENTUM  AND  ENERGY  INTEGRATIONS 

344  TMAS!*SI6*«H0*UNU«Ml  *DESIG*DEPM! 

TMt *TMASI*UN<1 *H)  * 

TNTPI*TNTLP*TMASl 
GO  TO  <284*278  )  «NVPAS 
284  l F  <  I 1  >  277*277*283 
283  IFtt-lSSYl  278.277*277 

277  TMASl =TMASI/2» 

TMl«TMI/2. 

TNTPI«TNTPl/2o 

278  TMAST  •  TMAST  4  TMASI 
TMlT  «  TMIT  ♦  TMl 
TNTPT  *  TNTPT  ♦  TNTPi 
IX-X/C01M42* 

233  GO  TO  < J66*203*209*212)*NAT 


TERMINATION  OF  -ANGLE  LOOP 


166  XOLDU  «M>  *  X 
L I NE3-L I NE3+ l 
LINE6-LINE6+I 
L I NET  «L I NE7+ 1 
IF<LINE3-A0)  241*241*244 

241  WRITE  OUTPUT  TAPE  NT3«300*BM«NC 

244  IFILINE6— 40*  242*242*245  • 

242  WRITE  OUTPUT  TAPE  NT6*300«BM*NC 

245  IF  (LINE7-40I  243*243,152 


-1L<- 


%% 


o  r»  non 


243  *rite  output  tape  nt?«3oo*bm.nc 

152  Phl«Pnl-*OPHl 

TERMINATION  OF  RAOIUS  LOOP 
S1G»$J6+DSIG-BuG 
236  B U-G = 0  » 

GO  TO  4998*998. 2881* IBP AS 
288  RADL  «  S1G 
MmADL  «  M 

bmolD  «  a« 

•GO  TO  <  326.325  1  *M8PAS 

325  BM  ■  BMINT 
'  GO  TO  327 

326  BM»-F«£OJ  *  8M 
„  LF1NL«2 

GO  TO  357 


ACOUSTIC  BANO  SEARCH  AND  AOOITION 


CALL  OROERS  (NOCFO.DSPS. ISPL1 

JK»1 

JMl«t 

jeno*o 

JDFSG*  1 0/0CTt> 

NOJ*NOCFO-t 

nojp*nocfo 

459  DO  453  JM«JMl*NOJ  ' 

1 K  <  D5PS  <  JM )  453.453*458 

456  jP*jm+1 

OO  451  J*  JP*  NO JP 
GO  10  (*57**t5c.  >  *J< 

457  14-  (DSPS  (0)1  *.5-.. 454*455  .  * 

454  N£  S t  G=NO  JP 
GO  TO  456 

455  ISP10*J5-PL(J)-1SPL<  JM) 

IF (ABSF ( ISPLDI-JOFSGI  452*452*451 

452  DSPS  (  J>'<  )  =DSPS  ( JM  l'rOSPS(J) 

I SPL ( JM ) - ( DSPS ( J ) *FlOATF ( l SPL ( J ) ) 4DSPS  < JM  >  4FL0ATF . 1 SPL ( J Mil/ (OSPS 
1  (  JM!-*0SPS(  J)  »  ’ 

1 SPL ( J  >  4-99 

DSPS ( J  5  »-l • 

JEN0=1 

451  CONTINUE 
JS 1 G* JMI 

IrtJEND)  374*374*453 

453  CONT INUE 
GO  TO  374 

456  N0J-NFS1G-1 

nojp*nfsig 

JM  1*2 
JK«2 

GO  TO  459 

374 

INNER  INDEX  UN  TERMINATION 


327  NMJ«ix=  (2O.*BM/R4DC-M0.  I/SET2  +1* 
J  F ( \Kj4X-2 0  I  600 . 600  *  60 1 


601  S£T2  «  S£72*2. 

I SET  «  tSET  +  I 
MRA0L* (MRADL+l 1/2 
NM00J  «  SET2 
60  TO  327 
NHOOJ  *  2  ♦*  1SET 

600  i SET I  »  LOOT  t  BMStCDIM  •  0*131 » 

SETU'-2**isET! 

IF  (NMOOM-B)  531*332*532 
531  NMDOM»  SET  1 2 
332  QEU  »  NMOOM 
SET t  *  1 SET1 

TERMINATION  OF  BM  LOOP 

WRITE  OUTPUT  TAPE  NT3.303 
*  00  141  1FRE0* 1 *NOCT 

WRITE  OUTPUT  TAPE  NT3.5.BM.SPSI 1FREQI * IFRE0.U81 
141  CONTINUE 

185  continue  This  aectlon  notf  operational 

band  assignment  and  shift 

383  00  366  IF»1.N0CT 
KSIG  •  0 

mixt*2o 

00  366  M«  1*MIXT 

CALL  ♦  ORDERS  «NMBM.OSPL*KAFP > 

00  386  KA«  2.NMBMP  . .  . 

MIXT>MIXM(MI 

NORAT  •  l 

IF  (0SPL(1F.KA«M)  391*395*393 

391  IF(IF-l)  395*393*665 

665  IF (DSPL I  IF— l «<A«M )  393*394*394 

393  IFUF-2J  395.393*666 

666  1F(0SPL< 1F-2.KA.M)  393*396*396 

396  NDRAT-2 

394  IF  ( IF- (NOCT-l ) )  667*395.393 

667  IFIDSPLUF+l  .KA.M)  397.398*398 

397  IF ( IF- (NOCT-2 ) I  668*395.393 
666  IFIDSPLUF-te.KA.M)  393*399*399 

395  GO  TO  392 

398  GO  TO  <661.6621*  NORAT 

661  DSPLUF.KA.M)  •  <DSPL<  IF-l.KA.M)  *  DSPLUF  +  1.KA.MU/2* 

662  OSPLUF.KA.M)  ■  IDSPH  1F-2.KA.M)  ♦  DSPLOF+l  .KA.MD/2* 

399  GO  TO  (663*664),  NDRAT 

663  OSPLUF.KA.M)  »  <  DSPL  <  IF-l  .KA.M  >  +DSPL  UF+2.  KA.M  )  )/2. 

664  DSPLUF.KA.M)  «  <DSPL<  IF-2.KA.M)  +  DSPL < IF+2 »KA.  M)  )/2* 

KAO  «  KAFP(KA-M  )-KAFP(KA>  • 

IF  (ABSF  (KAD)-l)  385.383*386 

385  KAPP  »  KAFP<KA+l) 

KSIG  •  l 

IF  ( DSPL ( 1 F+ 1 .KAPP . M ) )  387.387.388 

387  D5PLUF+1  .KAPP.M)  «  DSPL  (IF.  KAPP.  M  ) 

388  IFCOSPLUF-i  .KAPP.M)  369.389.391 

369  OSPLUF-1  .KAPP.M)  *  DSPlUF.KAPP  »1*M| 

392  0SPLUF.KAPP.-1  »M)  «  <DSPL<  IF.KAPP-l  *M)  +•  OSPlUF.kAPP.,1)  )/2 

386  CONTINUE 


IF(KSlG-n  390.383.383 
390  DO  384  IF-1.N0CT 
00  384  M«1.20 

00  384  KA-2.NM8M 

SPL  (IF)*  DSPLUF.KA.M)  +5PL(IF) 

SAFO  *SPL (IF >/0. 0000000000001  +1« 

OB  ( I F  )  *  1 0*l0GF  C  SAFO  )//0. 268 
364  WRITE  OUTPUT  TAPE  NT3.23. IF.FREi (IF;.DB(  :g» 
CALL  COPY < NT6. 1 i 
.CALL  C0PYINT7.il 
CALL  COPY  (NTS . 1 1 
END  FILE  NT4 
999  CONTINUE  • 

REWIND  NT4 * 

*  CALL  EXIT 


SUBROUTINE  FlMAINC.I.M.KAI 

COMMON  PI .P12.PRADL*PRAS0.FfiEQJ.IMH.BM»XE»NEST.NUM.LEAV.PMl.ALP0* 
1RA.S1 G. CD  I M< I  SAP  .BMlNT »  RADC » SF « SH.SJR.SJC 
DIMENSION  MMU  <  B  )  *KAL I8> .SF (30 ) «  SH(30) «SJQ 130 ) *  SJC (30) »SFP<  *  l.S 
l HP (  ,  >.SJCP<  «  J.SJRPl  •  » 

1F(NC-1 >99.20.50  ‘ 

20  !F(I-1 )  99.21*25 

21  IFfM-l)  99.22.28  > 

22  IAL-1 
1FIKA-D99.23.2T 

22-  FR|_0*L0GF  (FREOJ) 

RAB1-RADC/BMINT 

SET2«t. 

MMU  <  D  » 1 
IFfi-l 
MAXM-0 
MIT*1 

24  KAL (MIT) *LOGF (RAB 1* (SET2-.5 J/FRLO+l .9999 
IF (KAL (MIT >-NM8M >25.40.26 

25  M1T-MIT+1 
SET2-SET2+SET2 
GO  TO  24 

26  KALI  MIT >*NM6M 
GO  TO  40 

27  REWIND  NT5 
I  JU*NMDO J 
MLO*I 

GO  TO  34 

28  MIT  «  (M-l >*NMDOJ 

IF (Ml T  — MAXM )  30.30.29 

29  MAXM*MIT 
MML(ISET>*M 
IFC-l 

GO  TO  40 

30  MLO*MLO+MJU 

IF { MLO-MML ( I AL  >  32.  32.31 

31  lJU*MJU/2 
MuU-'MwO-I  JU-»<M-1  ) 

MJU=I JU*(3-H)+ (MML( IAL1-1 >/2 
I AL* I ALT l 

32  MIT *MJU— 1 

CALL  FILE<NT5.M1T) 

34  MJU* IJU 
IFC*0 

GO  TO  38 

35  IF( IFC  >99.36 

36  DO  37  MAT*2,IJU 

37  READ  TAPE  NT5.  MIT 

38  READ  TAPE  NTS .M J T .MAT . ( (SFP (KAP.K ) . SHPIKAP »K ) .SJCP (KAP. K > » SJRP (<AP 
l,K).KAP*MIT,MAT).K»2.NEST  > 

33  DO  39  K»2 . NEST 
SF  C  K 1  *  SFP ( K A .  K  ) 

SUR(0«SJRP<KA,K» 

SJC(K)=>SJCP<KA«K) 

39  SH(K)*SHP(KA.K> 

WRITE  TAPE  NT4 .  (SF  (K  ) . SH (K  )  .SJR( < ) • S JC (K ) .K*2«NEST 1 
RETURN 


*0  M2»N 
12«I 
1L  - 1  SET 

41  IF (MODF (M2«2)  99*42*44 
.  42  IF  <MODF (  12*2 )  99*43*44 

43  1L-IL+1 
M2-M2/2 
12-12/2 
GO  TO  41 

44  KALl-KAUlLl 
MlT-KALl-KA+1 

CALL  FLUX  t  Ml T ) * SFP «  SHP  «  S JCP * S  JRP 

WRITE  TAPE  NTS. KA, KALI.  t(SFP(KAP.K).SHP(KAP,Kl.SJCP<KAP,K»,SJRPl«A 
1P.K) »KAP«KA .KALI 1.K-2.NEST1 
IF ( 1— I SSV )  33.49*99 
49  END  FILE  NTS 
‘GO  TO  33 

90  READ  TAPE  WT4  «  (SFIKJ.SHOC) *SJR  OC 1 «  SJC <K ) .K-2.NEST) 

RETURN 
99  HALT 
END 
C 

c 

SUBROUT l NE  FLUX (MIT*  SFP . SHP . S JCP . S JRP  > 

COMMON  P I  .  P 1 2 .  PR  AOL  «  PRASO  .  FREQ  J  .  I MH .  0M .  XE  .  NEST  .  NUN  .LEA  V  •  PH  1  .  ALPO. 
1RA.SIG.CD1M.1SAP 

DIMENSION  SICt  ..SIR*  ««VI<  )«T1(  )*CMX(  ).VA0(  l.TAOl  )*SA 

1R0(  ) « SACO <  ).XE<30»  *NUM(30>,LEAV«301. 1SAPI30) •AlP0(30!«D£LPC80« 

230 ) «  RA  <80*301 . SFP (  .301.SHPI  *30 ) .SJCPJ.  «301*£JRP<  *30» 

DO  193  K*2*NEST 
NODE  a  2 
B0=8M 

00  201  18-1. NIT 
S1CUB1-0.0 
SIR ( 16 ) -0«0 
Vl (IB)-O.O 
TUI81-S.0 
CMX (IB) -BO+XE IK  1 
201  B0«B0#FREQJ 
NUMK-  NUMIK) 

leaf«=leav<k> 

1SAPE-ISAP<K) 

30  PCAL-PHI-ALPOIKJ 

31  ABPA=ABSF<PHAL) 

IF ( ABPA-Pl >34.34.32 

32  PHAL=PHAL-SI6NF(PI2*PHAL1 

33  GO  TO  31 

34  COSPA-COSF (ABPA) 

RADO-SQRTF  <  <RA  (NUMK  «K )-S l G  >*  IRA  (NUMK « K )-S I G  )+2»*RA  1NUMK.K  ) *S  1G* 

1 U.-COSPA)  1 

IF  (RADO-0.0005  *  COIN)  35*35.36 
RADO-O.O 

35  NODE  -  l 

GO  TO  (4A.36.44). ISAPE 

36  C0ST0- (SIG-RA 1NUMK.K  >  * COSPA 1/RAOO 
IF (ABSF(C0ST0>-1 .0 >38.37.37 

37  Tf;0-;.°i-E!GNP(Pl,C0ST0n/2. 


,150. 


c 

c 

c 


S INTO*  0.0 
GO  TO  39 

38  SINTQ  «  SIGNF  tSQSTF  tl«  -COSTO*COSTO» .PhALI 
40  TMO« AT ANF ( S iNTO/COSTO ) 

IF (COSTO)  151.152.152 

151  THO*THO+Pl 
GO  TO  154 

152  I F  < THO )  153. 15*«154 

153  TH0=TH0+P12 

154  GO  TO  (39.143.39). ISAPE 

CL*RAD0/»A(1«KJ 

CRA-RACl.K) 

CALL  CFLUX (MlT.Vl  »TI  .SIC. Cl. . CPA 1 

00  21 0.  1B«  1  .'HIT 

SlR< I8)»SIC< IB»*COSTO 

.  SlC(IB)*SIC(18)*51NT0 

SJC<1B)*S1R(IB)*S1NT0 
SIR< IB)«SI»( 1B)*C05T0 
GO  TO  143 

39  00  203  I8»l. HIT 
RADX  *  RADO/CMX (IB) 

VAOU0l«l.-exPF  (-flAOX*fiADX) 

GO  TO  (Al.A21.lMH 

41  TAO(  IB)  *  l  .-EXPF  {- (RAOXaPRADL  )*(  RADX*PRADL  )  ) 

Vz  si  « (O.S8623*ERROSF  (PAOX )  +  *ADX  *  (VAO(  !B)-1*0»*CM 

SARO  < IB  >  *SA  *COSTO 
SAC0(1B)«SA  *S I NTO 

203  CONTINUE 

leaf  integration 

44  DO  134  J*1*L£AF 
oo  134  n*:.nuhk 
PHAL  ■  PHAl.  -  OALP (N.ICI 
ABPA  *  ABSF  (PHAL) 

IF  (ABPA  -  PI)  50.50*48 
PHAC  «  PHAL.  —  SI  GNF  (P 12  .PHAL.  ) 

GO  TO  47 

RAlTi  SGRTF^ (PA(N.K)— S!G}*(WA(N«K>'»S16)+2.*PA<N»K>*SIG*C  1 


46 

47 

48 

49 

50 


1 COSPA ))  _ 

IF  (RAD  -0.0005  *  SIGJ  Pi »  51.  52 

51  NODE  «  1 

GO  TO  132  _ 

52  COST  *  (SIG-RA(N.O*COSPA)/«AO 
1F(ABSF(C0ST  1-1.0)54.53.53 

53  TH  »(P1-S1GNF (PI .COST  11/2* 

SINT  •  0.0 

54  S?NT°.6S1GN?(S0RTFU.  -COST  *  COST)  .PHAL) 
56  TH*ATANF ( SINT/COST ) 

] F (COST  )  161.1 62*152 

161  YHsTH+PI  * 

GO  TO  59 

162  I F ( TH )  163.59.59 

163  TH.TH+PI2 

59  GO  TO  (60.62) .NODE 


60  N00E  a  2 
OTH  *0.0 
GO  TO  121 

62  DTH  »  TH  -  THO 

63  ABOTH  a  ABSF(DTH) 

IF(ABOTH  -  OTHM) 121.121.65 

63  IF (ABOTH  -  »I >66.66.1©! 

C 

C  0  -  360  CORRECTION 
C 

101  THO  a  THO  +  SlGNFtPl2.DTH> 

106  GO  TO  62 

C 

C  INITIATION  OF  AUXILIARY  INTEGRATION  • 

C 

66  LQ*AB0TH/DTHM+1* 

*  Q  »  LQ 
68  OTHaDTH/O 

RCRC  a  RADO#COSTO  -  RAD*COST 

RCRC  =  RCRC  +  SIGNF(0. 00000001 «  RCRCJ 

70  ABLE* (RA00*S1NT0-RA0*SINT  J/RCRC 

71  BKR=.«A00*SINT0-ASLE*RA00*C0ST0 

c 

C  AUXILIAPY  integration 

c 

75  00  92  L-l.LO 

76  THaTHO  ♦  pTH 

77  COST  *  COSFITHl  ' 

79  SINT *S I NF {TH ) 

80  RAD*BKR/(S1nT-ABLE*C0ST > 

00  204  18*1. MIT 
RADX-RAD/CMXt 18) 

81  VA  a  1..-  EXPF  (-RADX  *  RAOX ) 

83  SA  *(0.88623  *  ERRORF(RADX)  +  RAOX  *  ( VA- T.O ) )*CN 
©A  SAR  a  SA  *  COST 
SAC  *  SA  *  SINT 

86  Vl  (  IBJ*Vi  <  I6)  +  <VA+VA0U8>  >*DTM 

89  SIC( I8)*SIC< 18  J  +  (SAC+SACO( IB>  >#DTH 
83  SlR(I8)=SIR<lB)+<  SAR+SARO ( IB  > )#DTH 

GO  TO  (90*91 >  »1MH 

90  TA«-EXPF(-RADX*PRAOL>* (RADX*PRADL>  1  +  1  • 

TI ( I  6 ) *T  I  ( IB )  +  ( T A+T AO ( IB) >»DTH 

TAO ( IB  >  «TA 

91  VAC ( 18 ) *VA ' 

SACO ( I B  )*SAC 
SA«0( IB  >*SAR 

204  CONTINUE 

92  THO  a  TH 

93  GO  TO  132 
C 

C  MAIN  LINE  INTEGRATION 

c 

121  00  203  18*1. MIT 
RAOX*RAO/CMX( IB  > 

VA  a-EXPF  (-RAOX  *  RAOX)  +  1. 

122  SA  *(0.68623  •*  ERRORF  (RAOX )  +  RAOX  *  (VA-l.O)) 

SAR  -  SA  *  COST 
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124  SAC  •  SA  *  SINT 

VI < 18>*VI <  JB)+ <VA+VAO< IB) )*OTH 

126  StR< 18) -SIR ( IB>+(SAR+SARO< IB) )*OTH 
S1C( IB)«SIC( I8)+(5AC+SACO< 18) )*OTH 
GO  TO  <127,128 )«IKH 

127  TA*-EXPF{-(RADX*PRADL)*(RADX*PRADL>  >+1.0 
Tl < IBl-T! < l8)+(TA+TA0l IB) )*DTH 

TA0{ 18)«TA 

128  VAO i  IB ) *»VA 
SACOt IB )*SAC 
SARO  < IB  )«SAR 

205  CONTINUE 
THO  »  TH 
132  RAOO  ■  RAO 
COSTO  *  COST 
134  S INTO  ■  SINT 
143  SFP  < I8iK)*VI < IB)  * 

CO  TO  (191*192)  IMH 

192  TI < IB  )*VI < IB) 

191  SHP  ( 1 8 • <) *T I  (  18)  ' 

SJCP(IB«K)«SJC<I8) 

193  5 JRP <1B*K)«$IR<18) 

RETUR> > 

END 


SUBROUTINE  CFLUX  <MIT.VS*TS.SSC*C_ 

COMMON  PI  *PJ2*PMADL*PRASQ*FR£QJ*  IMH«:;M*XE 

DIMENSION  ARR  140)  *VS<30  ) •  T5 (30 )'«SSC  (33  )  *  XC(20)  •  SFL<40) 

ISJ-1 

DL1-20, 

DU2«4. 

EXPOC" < 1 • -CL ) * < 1 • -CL ) 

CLl-U/CL 
CL2-CL+CL1 
CL1 *CL-CL1 
VL»ABSF<CL-1« ) 

OL "MAX 1 F < CL/PL l «  < CL- 1 • ) /OLE ) 

VLO-VL 

8*0, 16666667 
A»0. 33333333 
SFO»0« 

SFL (1 ) «0# 

INCREMENTS  AND  AREAS 
IF (CL- l • )  203*204*204 
204  NS«1 »/DL+*9999 
VN-NS 
DL- t  */VN 
ARR ( 1 )*0» 

CO  TO  206 

203  NS=CL/OL+. 99999 

IF (NS)  231*231*224 
224  VN“NS 

OL-CL/VN 
ARR  ( 1  )»M*EXPOC 
206  NS-NS+NS+i 

DC  215  IS-2*NS 
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VU*VL+OL 

VLS«VL#VL 

CPHI2*CL1/VL+VL/CU 

S PH I 2- SORTF <  4 . -CPH 1 2«CPHI 2 1 

PHI«ATANF(SPH12/CPKI2» 

IF<PHl)  214*215.215 

214  PHlsPHI+Pl 
CTHE2-CL2-VLS/CL 

THE- AT ANF ( SORTF  <4  *-CTHE2*CTKE2 ) /CTHE2 ) 

IF( THE)  217*275.275 
217  THE-THE+PI 

275  sfl<is)-<b*sfo+a*sphi2)*du  +sfu: is-a i 

SFO-SPHI2 

C«A 

A-B 

B«C 

215  APPUS»»VtS*PHI+TH£-VL/a.3*SPHI2 

C  *  FLUX  INTEGRATIONS  AT  SEQUENT  l  AL  B  PtANES 
BO-BM 

00  233  1S-1.M1T 
Vl-O. 

T 1  *0* 

SIC-O. 
bd-bo+xeiio 
IF(BO)  260.230 *281 
260  XSQ=a  000000. 

GO  TO  296 
281  Y. 3Q-CRA  /BO 
XSQ-XSO*XSQ 
201  BO=BO*FREQJ 
IS-ISJ+l 
ISO-1' 

0CJ-0I«*FL0ATF  (  1SJ  J 

VL-ABSF  <CL— 1 . 1  . 

VL-VL+Ol.J/2 
272  VUS=VU*VL 

EXPAR-VUS*XSO 
EXPV2-EXPF (-EXPAR/2* 1 
EXPV«EXPV2*EXPV2 
ARD-ARR(  1,S)-ARR(  ISO) 

VI«VI+ARO*EXPV 
GO  TO  <270.271 I.1MM 

270  T I *T I +ARD*EXPF ( -EXP AR-PRASO  > 

271  SIC  -SIC  +EXPV*U  .♦EXPAR+EXPAR-EXPV2 /*  ( SFL  (  IS  1-SFl.t  ISO  1  } 
I  SO- IS 

VL-VL+Dt-J 

IS=XMlNOF< IS+ISJ.NS) 

IF(ISO-NS)  272.296.296 
296  VS( IB  )=V1/PI*XSQ 
TS« IBI-TI/PI4XS0 
SSC ( IB  > -SI C/P I*XSQ 
J  F  < ARR  <1  )5  233.233.298 

298  VS  4  IB ) « 1 •— EXPF (-EXPOC*XSQ  )4VS(IB> 

CO  TO  (299.295). 1HH 

295  TS ( IB ) -VS< IB) 

GO  TO  233 

299  TSUB)  •  l.-EXPF<-EXP0C*XS3 
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*PRAS0)+TS<18» 


233  CONTINUE 
RETURN 


* 


LIST  OF  SYMBOLS 

area 

area  of  circle  having  a  radius  equal  to  R„  „ 

c,e 

nozzle  exit  area 
acoustic  constant 

constants 

frequency  constant 
length  constant 

porosity  parameter  for  attenuation 

enthalpy  flux  function  (Eq.  (7}J 
concentration  flux  function  (Eq.  (7)) 

Mach  number 
convective  Mach  number 

momentum  flux  function  (Eq.  (5)) 

pressure  (static) 
acoustic  power 

Prandtl  number 

geometric  radius  (circular  nozzle  exit  radius) 

flow- fie Id  coordinate 

largest  radius  of  the  mixed  flow  field 

gas  constant 

temperature 
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bounding  velocity 
convective  velocity 

*  i* 

volume 

attenuation  length  measure 

acoustic  velocity 
momentum  spreading  parameter 

enthalpy  spreading  parameter 

coefficient- 
specific  heat 

mean  specific  heat 
frequency 

frequency  at  point  -of  maximum  noise  production 
frequency  determined  by  Strouhal  relation 
enthalpy 

relative  enthalpy  (Eq.  (A. 17)) 
local  concentration 

eddy  heat-condution  coefficient  (Eq.  (A. 3)) 

radial  coordinate  (polar) (for  elemental  jet  flows) (Fig.  10) 

variable  of  Integration  (Fig.  10) 

polar  coordinate  (radial) (Fig.  10) 

temperature  (variable  of  integration) 

local  mean  velocity,  (x-directed) 

fluctuating  component  of  the  x-directed  velocity 

local  mean  velocity  (r-directed) 
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fluctuating  component  of  the  r-directed  velocity 
watts 

axial  coordinate 

space  correlated  mean  eddy  volume 

attenuation  energy  ratio,  coordinate  angle  (Fig.  10) 
Kirchhoff  attenuation  parameter 

ratio  of  specific  heats 
specific  heat  coefficients 
eddy  viscosity  coefficients 

density 
shear  stress 

angle  coordinate  (Fig.  10) 
velocity  parameter  (Eq.  (8b) ) 
eddy  characteristic  dimension 

Subscripts 


adjacent  gas  flows 
exit 

gas  (concentration)  index 
tensor  components 
ambient 
maximum 

exit  element  index 

relative  (enthalpy)  (Eq.  (A. 17)) 

vector  directions 

Superscripts 


stagnation 


mean 
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Figure  2.12b  -  Typical  input.  (Concluded) 
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Figure  2.13  -  Typical  Output 


Distance  Radius  Velocity  Shear  Enthalpy  Mach  Momentum  Shear 
Ratio  Ratio  Ratio  Ratio  Flux  _  Number  Flux  Flux 

45.1204  0.0452  CL'o'oolD  <£L.00CKD  1.000  CS.000Q>  Cf.6bdg>  CgT.ObbjE 

.0500  .2690  .0510  1774139  .9920  0  0 

45.1204  -.0122  .0409  0.0000  • 

Frequency  Sound  Power 

1  1.250E+02  O.OOOE-OO 

2  1.574E+02  O.OOOE-OO  . 

3  1.984E+02  C.OOOE-OO 

4  2.499E+02  1.281E-01 

5  3.149E+02  2.896E-01 

6  3 . 968E+02  O.OOOE-OO 

7  4.999E+02  5.445E-01 

8  6. 299E+02  8.029E-01 

9  7.937E+02  O.OOOE-OO 

10  9-999E+02  8.989E-01 

11  1. 259E+03  8 . 248E-01 

12  1 . 587E+03  7..191E-01 

13  1. 999E+03  6.402E-01 

14  2.3 19E+03  5.794E-01 

15  3.174E+C3  5.332E-01 

16  3.999S+03  O.OOOE-OO 

17  5.039E+03  5.984E-01 

18  6.349E+03  4.734E-C1 

19  7.999E+03  4.394E-01 

20  1.007E+04  4.044E-01 

21  1. 269E+04  3.727E-01 

22  1.599E+04  3.616E-01 

23  2.015E+04  3.807E-01 

Frequency  Sound  Power  DB  24  2.533E+04  4.152E-01 


1  1.250E+02  O.OOOE-OO  0.0000 

2  1.574E+02  4.912E-02  116.8336 

3  1 . 984E+02  O.OOOE-OO  0.0000 

4  2.499E+02  1.281E-01  120.9943 

5  3 . 149E+02  2.896E-01  124.5344 

3 . 968E+02  O.OOOE-OO  0.0000 
4.999E+02  5.445E-01  127.2737 

6. 29SE+02  8.029E-01  128.9594 
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3-0  Investigation  of  the  Effects  of  Flow  Temperature  and  Pressure  on  Jet 

Soise  Generation  and  Suppression 

3-1  Introduction 

In  order  to  apply  jet  noise  suppression  concepts  to  afterburning  engines, 
an  understanding  of  jet  noise  generation  and  suppression  at  high  flew  tempera- 
t ures  must  be  developed  beyond  the  present  state-of-the-art.  The  validity  of 
the  correlation  between  free  jet  mean -flow  velocity  profiles  and  noise  genera¬ 
tion  and  suppression  by  heated  jets  has  been  demonstrated  in  Section  1  of  this 
report,  but  the  flow  temperature  range  investigated  does  not  extend  to  levels 
achieved  by  afterburners.  This  section  describes  the  results  of  acoustical 
investigations,  carried  out  in  addition  to  the  original  scope  of  the  res*:erch 
program,  for  the  purpose  of  providing  information  on  noise  generation  by  high 
tarperature  jets,  far  both  simple  circular  nozzles  and  suppressor  configura¬ 
tions.  The  results  of  this  study  provide  basic  information  on  noise  generation 

* 

by  high  temperature  jets,  and  give  further  insight  to  the  problem  cf  evaluating* 
noise  suppressor  performance  by  means  of  cold  flow  scale  model  tests. 

Experimental  investigations  have  been  carried  ,ut  to  evaluate  the  'jet 
noise  generated  by  a  conical  nozzle  over  a  range  of  flow  temperatures  and 
pressure  ratios,  extending  up  to  a  temperature  of  3300*R  and  pressuie  ratio 
3.u.  The  jet  noise  spectra  and  overall  sound  power  levels  are  presented,  and 
the  effect  of  jet  density  on  noise  generation  is  described.  The  results  of 
the  high  temperature  scale  model  conical  nozzle  tests,  as  well  as  concentric 
jet-  and  full  scale  engine  data  have  been  vised  to  establish  an  empirical  rela¬ 
tionship  between  nozzle  energy  flux  and  total  sound  power  generated,  applicable 
to  concentric  jet  flows  of  unequal  velocities  and  temperatures  as  well  as  to 
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conics!  nozzles  at  all  pressure  ratios  end  temperatures. 

Scale  model  acoustical  evaluations  of  suppressor-nozzle  noise  reduction 
performance  at  high  temperatures  have  been  accomplished,  vith  results  indicating 
that  jet  noise  suppressor  action  is  temperature  as  well  as  pressure  ratio 
dependent.  Data  obtained  indicates  that  noise  reduction  may  be  more  difficult 
to  achieve  at  high  flow  temperature,  ana  that  suppressor  scale  model  evalua¬ 
tions  must  be  carried  out  vith  heated  flow  or  suitable  substitute  gas.  Further 
experimental  actii  .ty  in  the  area  of  jet  noise  generation  and  suppression- at 
flow  temperatures  in  the  afterburner  range  is  recormnended  for  tbe  purpose  of 
defining  generalised  temperature  scaling  laws  applicable  to  all  nozzle  config¬ 
urations  . 

3.2  Noise  Generation  by  High  Temperature  Jets 

Experimental  investigation”  on  noise  generation  by  jets  operating  at 
temperatures  up  to  3300°R  have  been  carried  out  utilizing  a  free-field  acoustic 
test  facility  (Figure  3.1),  cleccri  ed  in  Section  1.3.1  of  this  report.  Tests 
were  conducted  vith  a  conical  nozzle  3.85  inches  in  diameter  at  the  exit, 
operating  at  pressure  ratios  of  1.7,  2.0,  2.5,  3.0,  and  3«l»,  varying  the  flow 
total  temperature  from  560°R  to  3300°R  for  each  pressure  ratio  condition.  Sound 
pressure  spectra  were  measured  along  a  circular  arc  at  35  feet  from  the  nozzle 
exit,  and  at  various  angles  from  the  discharge  direction.  Sound  power  speerta 
were  computed  cy  integration  of  the  sound  pressures  over  the  hemispherical 
radiation  surface. 

3,2.1  Instrumentation 

The  flew  conditions  which  were  controlled  and  monitored  during  the 
experimental  investigation  were  total  pressure  and  total  temperature  at  the 
inlet  to  the  nozzle.  An  8  inch  diameter  flow  section  preceding  the  nozzle 
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inlet  contained  a  static  pressure  tap  and  three  IR  vs  IR-Rh  thermocouples 
which  were  continuously  monitored  during  the  tests.  The  area  ratio  of 
approximately  k’l  from  the  measuring  section  to  the  nozzle  throat  provided  an 
effective  plenum  upsream  of  the  nozzle  inlet. 

Acoustic  measurements  were  made  using  a  Bruel  and  Kjaer  type  I4I33  condenser 
microphone  system  which  was  calibrated  by  means  of  a  B  &  K  pistonphone.  The 
microphone  system  frequency  response  was  flat  to  liO  kcps,  covering  the  entire 
frequency  range  of  measurements  without  need  of  response  correction.  Microphone 
output  was  filtered  through  a  B  &  K  1/3  octave  spectrometer,  over  a  range  of 
band  center  frequencies  from  200  to  32,000  cp3.  Graphic  representation  of  the 
sound  pressure  spectra  was  obtained  by  recording  spectrometer  output  with  a 
B  &  K* level  recorder.  Sound  power  spectra  and  overall  PWL  were  computed  from 
the  pressure  level  spectra  obtained  at  a  distance  of  1$  feet  from  the  nozzle 
and  at  angles  of  20,  30,  I4O,  $0,  60,  75,  90,  110,  130,  and  150°  from  the  flow 
direction* 

3,2,2  Test  Results 

The  measured  total  acoustic  power  levels  and  sound  power  spectra  are 
presented  in  Figures  3.2  and  3.3a  through  3.3d,  for  the  range  of  nozzle  flew 
conditions  investigated.  In  Figure  3*2  the  overall  sound  power  levels  are 
plotted  against  flow  velocity  for  the  purpose  of  comparison  with  the  Lighthill 
eighth-power  parameter:  k°aAU®c0~'’  ,  The  flow  velocity  used  is  a  calculated 
value,  assuming  that  the  flow  fully  expands  isentropicly  from  the  measured 
total  pressure  to  the  ambient  static  pressure  of  30  in,  Hg,  As  would  be 
expected,  excellent  agreement  with  the  Lighthill  equation  is  obtained  at 
moderate  pressure  ratios  (up  to  2,0)  and  relatively  low  temperatures  (up  to 
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approximately  1?00°R).  The  limits  of  applicability  of  the  Lighthill  parameter 
are  clearly  visible  in  the  figure,  and  a  flow  temperature  influence,  apart 
from  velocity,  is  apparent.  Thus  it  is  seen  that  for  estimates  of  jet  noise 
total  power  at  high  flow  pressure  ratio  and/or  temperature  the  Lighthill 
relationship  is  inadequate.  Figure  3.2  shows  that  at  low  flew  temperature 
(560°R)  sound  power  generation  increases  a3  U2®  over  the  range  of  pressure 
ratios  from  2.0  to  3.ii»  while  at  a  flew  temperature  of  330O°R  the  overall 
sound  power  is  proportional  to  approximately  U^.  Examination  of  the  increase 
in  sound  power  resulting  from  velocity  increase  due  to  flow  temperature  rise 
at  constant  pressure  ratio  shows  that  at  a  pressure  ratio  of  3.0  sound  power 
increase  is  proportional  to  U2  from  1100  to  3300°R. 

This  variation  of  the  velocity  exponent  from  2  to  20  suggests  that  the 
Lighthill  parameter  might  take  the  form:  k(?j/6a)nAU®c“^  .  In  Figure  3.U 

the  effect  of  jet  density  on  sound  power  generation  is  shown.  Over  the  entire 
range  of  flow  conditions  investigated  the  sound  power  is  observed  to  increase 
in  'proportion  to  jet  density  squared,  considering  the  data  in  a  rather  gross 
manner.  A  detailed  observation,  however,  indicates  that  for  a  given  flow 
temperature  the  sound  power  varies  with  (£j/?a)n  ,  where  n  is  temperature 
dependent.  The  magnitude  of  this  temperature  dependence  is  demonstrated  by 
the  substantial  change  in  slope  of  the  sound  power  ^s  jet  density  curve, 
progressing  from  temperature  of  5&3  to  3300°R.  The  jet  density  ratio  exponent 
n  is  plotted  as  a  function  of  flow  total  temperature  in  Figure  3.5.  Values 
of  n  vary  from  *12  to  -6,  and  it  is  observed  that  the  relationship  between 
n  and  flow  temperature  is  approximately:  n  ,  For  low  pressure  ratio 

jet  flow  at  moderate  temperatures  the  flow  density  effect  on  jet  noise  generation 
is  negligible.  However,  the  results  of  the  investigations  described  herein 
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extend  knowledge  of  the  behavior  of  jet  noise  generation  to  a  broader  range  of 
flow  conditions,  including  high  temperature  and  high  pressure  ratio  flow  which 
is  likely  to  be  experienced  in  advanced  propulsion  systems. 

3,3  Correlation  of  Jet  Noise  with  Nozzle  Energy  Flux 

Results  of  the  scale  model  investigations  on  noise  generation  by  simple 
circular  jets  operating  at  temperatures  up  to  33DO°R  have  been  used  to  establish 
a  general  empirical  relationship  between  total  acoustic  power  and  jet  stream 
power. 

Figure  3.6  shows  normalised  power  level,  (PWL  -  lOLog^gt),  plotted  as  a 
function  of  noaxle  total  energy  flux.  The  scale  model  conical  noxxle  data 
presented  is  that  described  in  Section  3.2,  and  turbojet  «*gine  noise  data  is 
included  to  demonstrate  the  generality  of  the  relationship.  For  simple  conical 
noxxles  the  energy  flux,  €  ,  is  given  by: 

€  -  ht(W/A)  (1) 

where  hj.  is  the  flow  total  enthalpy  per  pound  of  fluid. 

A  tentative  method  for  ineludiiig  noise  generation  by  concentric  flows  in 
this  correlation  has  also  been  developed,  and  three  concentric  flow  examples 
are  included  in  Figure  3.6.  A  table  of  nozzle  flow  conditions  is  also  given 
in  the  figure.  The  dual  stream  concentric  flow  nozzle  utilized  in  obtaining  the 
noise  data  is  shown  in  Figure  3,7.  Flow  parameter  variations  tested  included  a 
heated  high  velocity  outer  flow  with  cold  inner  flew,  and  a  case  with  nearly 
equal  flow  velocities  but  with  substantial  temperature  difference  between  the 
two  streams.  The  confirmation  of  the  validity  of  the  correlation  for  concentric 
flow  requires  additional  experimental  investigation,  but  results  obtained  thus 
far  are  extremely  encouraging.  The  correlation  developed  is  as  follows: 
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For  a  single  stream,  Eq.  (6)  reduces  to  the  simple  form  of  Eq.  (1).  The 
significance  of  Eq.  (6)  is  that  for  consideration  of  concentric  flows  it  it 
necessary  to  sum  the  energy  flux  for  each  of  the  jets.  The  inner  jet  flux  is 
defined  herein  as  the  relative  flux  with  respect  to  the  surrounding  jetj  thus 
the  inner  flux  is  based  upon  the  static  enthalpy  plus  the  relative  velocity 
head. 

The  efficiency  of  conversion  of  jet  stream  power  to  acoustical  power  is 
seen  in  Figure  3.6  to  be  uniquely  related  to  nozzle  energy  flux.  Good  correlation 
between  model  and  full  scale  engine  data  has  been  demonstrated,  thus  indicating 
that  the  nozzle  energy  flux  parameter  is  quite  general  and  valid  In  application 
to  the  problem  of  jet  noise  estimation, 

3.1t  Effects  of  Flow  Temperature  on  Jet  Noise  Suppressor  Action 

The  value  of  cold  flow  scale  model  investigations  in  the  development  of 
jet  noise  suppressors  is  well  recognized,  both  from  the  standpoint  of  time 
requirements  and  hardware  expense  and  complexity.  Up  to  the  present  time,  jet 
noise  suppressor  investigations  have  been  carried  out  under  the  assumption  that 
noise  reduction  achieved  in  cold  flow  tests  can  be  directly  related  to  expected 
full  scale  suppressor  noise  reduction  by  Strouhal  number  scaling.  Furthermore, 
experience  has  shown  that  jet  noise  suppressors  in  general  perform  poorly  at 
lew  pressure  ratios,  with  noise  reduction  performance  (compared  to  an  equal 
area  conical  nozzle)  improving  as  pressure  ratio  increases.  In  order  to  apply 
scale  model  techniques  to  the  development  of  a  suppressor  for  afterburning 
turbojet  engines,  it  was  deemed  necessary  to  investigate  the  effects  of  hi$x 
temperature  flow  on  suppressor  performance  in  order  to  determine  the  validity 
of  scaling  laws  arid  the  Mach  number  effect  on  suppression. 
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Scale  model  suppressor  noise  investigations  were  carried  out  utilizing 
the  nozzles  shown  in  Figures  3.8,  3.9.,  and  3.10.  These  suppressor  nozzles  are 
of  approximately  the  same  flow  area  as  the  conical  nozzle  described  previously 
in  this  section,  and  are  described  in  detail  in  section  1  of  this  report. 
Measurements  of  total  acoustic  power  were  made  in  the  same  manner  as  for  the 
conical  nozzle,  though  over  a  somewhat  more  limited  range  of  flow  pressure 
ratio  and  temperature. 

Comparisons  of  total  acoustic  power  generation  have  been  made  for  the 
suppressors  arid  conical  nozzle.  For  convenience,  the  conical  nozzle  data  is 
replotted  in  terms  of  PWL  vs  >ressure  ratio,  for  various  temperatures,  in 
Figure  3.11.  Figure  3.12  presents  corresponding  data  for  the  three  suppressor 
models:  (l)  8-Iobe,  (2)  19-tube,  and  (3)  l8-segn»ent  with  plug.  In  order  to 
determine  the  effect  of  flow  temperature  on  suppression  of  jet  noise,  the  data 
of  Figures  3.31  and  3,12  has  been  cross  plotted  in  Figure  3.13,  which  shows 
overall  FWL  as  a  function  of  flow  temperature  for  the  four  nozzles  operating 
at  three  pressure  ratios:  1,8,  2, It,  and  3.0.  It  can  be  readily  observed  that 
for  each  of  the  suppressor  nozzle3  the  noise  reduction  (the  difference  between 
the  conical  nozzle  and  suppressor  nozzle  curves)  at  constant  pressure  ratio  is 
temperature  dependent.  Furthermore,  this  temperature  dependence  appears  to  be 
different  for  each  suppressor.  Note  that  at  pressure  ratio  1.8  little  noise 
reduction  is  shewn  for  cold  Flow,  and  that  as  temperature  is  increased  the 
suppression  capability  increases  and  then  finally  declines.  At  flow  temperatures 
near  1000°F,  only  slight  Mach  number  (pressure  ratio)  effect  is  observed,  and 
for  pressure  ratio  2«U  little  temperature  effect  is  evidenced.  It  is  fortuitous 
that  turbojet  engine  suppressor  development  up  to  the  present  time  has  been  for 
engines  operating  at  pressure  ratios  near  2 and  temperatures  near  1000°F, 
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since  noise  reduction  at  pressure  i*atio  2.1*  shows  little  temperature  dependence. 
However,  at  low  pressure  ratio  (1.8)  cold  flow  test  results  would  not  provide 
a  satisfactory  estimate  of  sound  power  reduction  experienced  at  increased 
temperatures. 

While  no  general  means  of  correlating  suppressor  performance  for  all  types 
of  suppressors,  for  various  temperatures  and  pressures,  has  been  developed, 
the  results  of  the  investigation  have  several  important  implications.  First, 
evidence  has  been  presented  to  show  that  jet  noise  reduction  is  not  simply 
pressure  ratio  dependent  for  a  fixed  configuration.  Furthermore,  the  variation 
in  noise  reduction  performance  as  a  function  of  flow  temperature  for  different 
suppressor  configurations  indicates  that  model  suppressor  evaluation  must  be 
ultimately  performed  experimentally  at  flow  temperatures  and  pressures 
representing  full  scale  conditions*  Finally,  the  need  for  further  research  into 
normalization  of  suppressor  performance  for  temperature  and  pressure  effects 
is  emphasized. 

3.5  Conclusions  and  Recommendations 

Results  of  scale  model  experimental  investigations  on  noise  generation  by 
jets  from  conical  and  suppressor  nozzles  have  been  presented  for  flow  temperatures 
ranging  up  to  3300°R.  The  following  important  conclusions  are  summarized: 

1)  iioise  generation  by  conical  jets  is  not  adequately  described  by  the 
Light hill  equation  except  for  low  pressure  ratio  and  moderate  temperature  flow. 

2)  A  temperature  dependent  flow  density  effect  on  noise  generation  by 

conical  jets  has  been  found  to  occur.  Total  sound  power  varies  with  , 

where  n  is  temperature  dependent, 

3)  An  empirical  relationship  between  sound  power  and  nozzle  energy  flux 
has  been  developed  for  conical  nozzles,  including  concentric  flows,  for  the 
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m^-£cm»acJ-fc  v<ac  nm: 


i  i 


range  of  pressure  ratios  from  U,  to  JA.  -  **«**”’  “  35°°^ 

U  let  noice  suppressor  acousticsl  performance  is  tesperatar.  as  well  as 

pressure  ratio  dependent,  so  that  scale  node!  oppressor  development  s»st 
include  final  Rental  evaluation  under  flop  conditions  representing  full 

scale  engine  temperature  as  well  as  pressure  ratio. 

»  is  recommended  that  further  research  effort  be  directed  toward  lending 
the  generality  of  the  energy  flux  except  to  c^lex  nestle  designs,  includirg 
further  verification  of  applicability  or  -centric  flows,  let  noise  sup^«si* 
research  should  bs  oontdnued  for  the  purpose  of  developing  general  laws  for 
scaling  temperature  and  pressure  ratio  effects  on  noise  reduction  performance 
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Figure  3,1  High  Temperature  Test  Facility  for  Jet  Noise  Evaluation 
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Correlation 


Figure  3.7  Concentric  Flow  Hotzle 
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Figure  3.13  -  Effect  of  Flow  Temperature  and  Pressure  Ratio  on  Suppressor 
Acoustic  Performance r 


4.0  SXPERIMSIfEAL  STiffY  ON  TRANSMISSION  AND  DISSIPATION  OF  SOUND  THROUGH  A 
TURBULENT  JET  WAKE 


4.1  Introduction 

A  major  assumption  of  the  jet  noise  suppression  theory  relating  the 
aerodynamic  and  acoustic  properties  of  jet  flow  is  that  ail  of  the  stream¬ 
generated  noise  radiates,  without  significant  loss,  into  the  far  field.  For 
complex  suppressor -nozzle  designs,  the  validity  of  this  assumption  can  he 
questioned,  since  sound  sources  associated  with  the  inner  flow  elements  are 
surrounded  by  high  velocity,  turbulent  flow.  There  exists  some  experimental 
evidence  that  acoustic  energy  is  dissipated  in  the  highly  turbulent  jet  mixing 
region.  In  Reference  1,  jet  noise  cata  is  presented  showing  that  sound  levels 
measured  Id  the  far  field  of  a  rectangular  array  of  nozzles  is  higher  when  the 
measuring  location  is  oriented  along  the  long  side  of  the  multiple  tube  config¬ 
uration  as  compared  to  the  short  end  of  the  array  (Figure  4.l).  This  effect, 
referred  to  as  shielding,  indicates  that  perhaps  the  sound  generated  by  inner¬ 
most  elements  of  a  suppressor  iu  unable  to  propagate  through  the  turbulent 
regions  created  by  surrounding  suppressor  elements.  Although  theoretical 
justification  for  significant  absorption  of  audio  frequency  sound  in  a  jet  wake 
is  lacking,  it  is  important  that  the  possibility  be  investigated  since  positive 
evidence  of  sound  dissipation  would  require  modification  of  the  correlation 
between  acoustic  and  aerodynamic  properties  of  the  jet  wake.  Furthermore, 
absorption  of  acoustic,  energy  in  the  jet  wake  could  explain  in  part  the  increase 
in  noise  reduction  by  a  suppressor  as  pressure  ratio  is  increased. 

Experiments  lave  t>een  conducted  in  an  effort  to  determine  whether  or 
not  dissipation  of  sound  does  occur,  or  whether  the  noise  reduction  effect 


referred  to  as  shielding  is  a  directivity  rather  than  power  level  change  effect. 


4.2  Experimental  Procedure 

The  determination  of  the  magnitude  of  sound  energy  absorption  in 
turbulent  flow,  as  related  to  jet  noise  suppression,  requires  that  the  power 
level  of  a  source  be  determined  after  propagation  through  a  turbulent  jet  wake 
compared  to  the  source  power  level  in  the  absence  of  the  turbulent  jet. 

Ideally,  n  discrete  frequency  sound  source,  which  could  easily  be  distinguished 
from  oroad-band  jet  noise,  could  be  located  within  a  jet  wake,  and  sound  power 
level  change  could  he  determined  in  a  reverberant  room  for  various  jet  veloci¬ 
ties.  However,  the  problem  of  locating  such  a  sound  source  within  a  jet  wake 
without  altering  its  acoustic  output  is  quite  difficult. 

An  alternative  method  for  evaluating  sound  absorption  by  the  jet  wake 
involves  placement  of  the  sound  source  outside  the  jet  boundary  while  posi¬ 
tioning  the  acoustic  pickup  within  the  wake.  While  this  scheme  eliminates  the 
source  variation  problem  mentioned  above,  it  presents  the  difficult  problem  of 
extracting  the  pure  tone  source  signal  from  the  jet  noise  and  jet  pressure 
fluctuations  especially  for  high  Mach  flow.  Even  more  difficult  would  he  the 
evaluation  of  convective  effects  on  the  sound  propagating  through  the  jet  flow. 
Since  the  sound  path  through  the  wake  would  depend  on  flow  velocity,  comparison 
of  levels  measured  with  and  without  flow  would  he  exceedingly  difficult. 

In  a  third  possible  approach  to  the  experimental  investigation,  the 
high  intensity  pure  tone  sound  source  can  he  located  outside  the  jet  wake,  and 
a  survey  of  sound  pressure  level  is  made  by  traversing  a  microphone  along  the 
Jet  boundary,  but  not  in  the  flow,  on  the  opposite  side  of  the  wake.  Free 


field  acoustic  environment  is  required,  and  effects  of  changing  sound  path  length 

due  to  convection  of  the  signal  in  the  jet  flow  can  he  minimized  by  locating 
« 

the  sound  source  sufficiently  f«-**  from  the  jet  wake.  Because  of  the  decreasing 
flow  velocity  in  the  jet  in  the  downs  bream  axial  direction,  a  focusing  of  the 
sound  propagating  through  the  flow  would  occur  as  described  schematically  in 
Figure  4.2.  The  sound  waves  entering  the  wake  in  the  high  flow  velocity  region 
near  the  nozzle  exit  are  convected  further  downstream  than  waves  entering  in 
the  lower  velocity  regions  of  the  jet,  thus  forming  the  shadow  zone  near  the 
nozzle  exit  and  somewhat  fccuoing  the  rays  at  a  downstream  location.  In  order 
to  evaluate  the  possibility  of  sound  dissipation,  the  intensity  of  the  trans¬ 
mitted  sound  must  be  measured  over  a  suitably  long  region  of  the  jet  so  that  a 
valid  average  sound  pressure  can  be  determined,  considering  both  the  shadow 
zone  and  region  of  apparent  amplification. 

Experimental  investigations  were  conducted  following  the  three 
previously  described  procedures,  but  concentrating  on  the  first  and  third 
methods.  The  second  scheme,  with  the  sound  pickup  within  the  jet  wake, 
presented  considerable  difficulty  in  extracting  the  pure  tone  sound  from  the 
jet  pressure  fluctuation  and  thus  it  was  abandoned. 

The  acoustic  signal  for  which  attenuation  effects  were  to  be  eval¬ 
uated  vac  produced  by  a  Harunan  generator  (Figure  4-3).  The  device  consists  of 
a  3 /B  inch  diameter  convergent  nozzle  which  discharges  at  pressure  ratio  of  4.5 
into  a  l/2  inch  diameter  cavity.  The  generator  produces  a  high  intensity 
discrete  tone.  The  frequency  of  the  whistle  is  a  function  of  the  cavity  depth 
which  can  be  varied  by  a  piston  in  the  cavity.  The  acoustic  power  level 
produced  is  on  the  order  of  150  db  (re  10~°  watts)  over  a  frequency  range  of 
1200  to  9000  cps. 
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Experiments  on  the  dissipation  of  sound  with  the  source  within  the 
jet  were  carried  out  using  the  apparatus  shown  in  Figure  h.k.  The  Hartman 
generator  was  located  in  a  cylinder  within  an  annular  jet  so  that  the  jet  flow 
would  not  affect  its  acoustic  output  Measurements  of  reverberation  sound  pres¬ 
sure  level  were  made  in  a  reverberant  room,  for  various  source  frequencies  and 
nozzle  pressure  ratios.  Average  pure  tone  sound  levels  were  determined  hy 
slowly  moving  the  microphone  within  the  room  and  reading  its  output,  separated 
from  jet  noise  by  means  of  a  narrow  band  filter,  with  an  acoustic  integrator 
(Ref.  2).  Static  pressure  was  measured  within  the  Hartman  generator  enclosure, 
and  the  maximum  variation  from  ambient  was  2.5"  H^O  at  nozzle  pressure  ratio 
P/P0 . 2.2.  Tills  erriall  change  in  pressure  could  not  affect  the  whistle  perfor¬ 
mance.  Results  of  the  invest  Igation  are  shown  in  Figure  b.5,  where  reverbera¬ 
tion  sound  pressure  level  for  the  whistle  tone  is  plotted  against  nozzle  pres¬ 
sure  ratio  for  various  pure  tone  frequencies.  A  small  decrease  in  sound  level 
is  observed  as  the  jet  flow  begins,  but  this  reduction  does  not  appear  to 
increase  as  flow  pressure  ratio  is  raised.  If  the  sound  energy  is  dissipated 
in  the  turbulent  flow,  then  the  tore  turb'ilent  higher  pressure  ratio  flow 
should  provide  increased  noise  reduction.  Since  this  is  not  the  case,  a  plau¬ 
sible  explanation  for  the  slight  decrease  in  sound  pressure  level  is  that  the 
supply  pressure  to  the  Hartman  generator  decreases  slightly  as  nozzle  airflow 
is  initiated.  Since  rather  small  supply  pressure  variation  can  affect  the 
whistle  output,  this  is  believed  to  be  the  source  of  the  noise  reduction.  In 
the  succeeding  experiment,  supply  pressure  was  monitored  much  more  closely  so 
as  to  preclude  such  an  occurrence. 
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Further  experiments  were  conducted  to  determine  if  sound  in  the 
audible  frequency  range  is  attenuated  during  propagation  through  the  turbulent 
flow  exhausting  from  a  nozzle.  As  in  the  previously  described  work,  a  high 
intensity,  pure  tone  sound  source  (Hartman  Generator)  provided  the  acoustic 
signal,  and  a  3~l/2  inch  diameter  conical  nozzle  was  utilized  to  create  the 
turbulent  wake.  The  experimental  set  up  was  as  Shown  in  Figure  4.6.  The  sound 
source  was  located  k-l/2  nozzle  diameters  downstream  of  the  nozzle  exit,  »nri 
several  diameters  to  the  side  of  the  stream.  On  the  opposite  side  at  the  Jet 
wake  a  microphone  was  traversed  over  a  distance  of  approximately  4  nozzle 
diameters.  Thus,  by  positioning  the  turbulent  flow  between  sound  source  and 
acoustic  pickup  its  effect  on  the  acoustic  signal  could  he  determined.  In 
order  to  measure  only  the  direct  sound  and  eliminate  reverberation,  the 
Acoustics  Laboratory  anechoic  room  was  utilized.  The  room  is  capable  of 
handling  the  airflow  from  the  3-1/2  inch  diameter  nozzle  operatic  at  pressure 
ratio  up  to  2.4.  To  test  whether  attenuation  tffirJ&gh  the  stream  could  be 
measured  in  this  facility,  a  mechanical  sound  barrier  6”  in  diameter  was  posi¬ 
tioned  between  sound  source  and  microphone  in  the  approximate  location  of  the 
jet  wake.  The  attenuation  measured  using  the  mechanical  barrier  in  place  of 
the  jet  wake  was  13  db;  sufficient  to  indicate  that  reflected  and  refracted 
sound  is  at  a  level  low  enough  so  as  not  to  obscure  attenuation  throug .  the 
stream,  which  might  be  of  lesser  magnitude. 

Source  frequencies  of  2.5,  5-0,  and  9-0  keps  were  used,  and  nozzle 
pressure  ratios  up  to  2.4  were  tested.  So  as  to  obtain  a  usable  signal  to 
noise  ratio,  a  Panoramic  Sonic  Analyzer  with  25  cps  filter  band  width  was  used 
as  a  signal  filter  for  the  microphone  output. 
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A  plot  of  sound  pressure  level  along  the  microphone  traverse  path, 
for  three  nozzle  pressure  ratios,  is  shown  in  Figure  4.7.  Note  that  "shielding" 
occurs  in  the  region  near  the  nozzle  exit,  and  amplification  is  experienced 
further  downstream  due  to  convection  of  the  signal  in  the  jet  wake-.  Results 
indicate  that  the  peak  amplitude  measured  is  independent  of  flow  velocity,  but 
that  the  flow  has  a  pronounced  influence  on  directivity.  Thus,  if  dissipation 
of  sound  does  occur,  it  is  insignificant  in  magnitude  in  the  audible  frequency 
range.  The  results  indicate  that  the  "shielding"  of  noise  from  one  jet  by  the 
flew  from  another  as  reported  by  other  researches  is  a  directivity  effect,  and 
that  the  decrease  in  noise  experienced  at  one  position  is  accompanied  by  an 
increase  where  the  sound  is  convected  downstream. 

4.3  Conclusions 

The  results  of  the  experiments  described  indicate  that  little,  if  any, 
sound  is  absorbed  by  the  turbulent  jet  wake,  in  agreement  with  the  work  in 
Ref.  3-  Thus,  no  modification  of  the  correlation  of  aerodynamic  and  acoustic 
jet  properties  is  required  from  this  consideration. 

The  shielding  effect  of  a  jet  wake  has  been  shown  to  be  a  result  of 
directivity  rather  than  sound  power  change.  The  problem  of  measuring  the  sound 
power  of  a  jet  in  the  presence  of  a  shielding  wake,  requiring  exceedingly 
precise  acoustical  measurements,  has  been  bypassed  by  substitution  of  a  pure 
tone  sound  source  for  the  jet  noise  source. 
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,  v*  O  -v1  ■ 


Hoiss  Source  within  Jet  Stream 


fifura  U.T 


I 

! 

* 

j 

I 

i 

I 

i 

1 

I 

1 

1 

i 

a. 

i 


j, 

« 

i 

i 

I 

I 

I 


5.0  *N  EVAUJATION  OF  SELECTIVE  HATER  INJECTION  AS  i  JgT  NOISE  SUPPRESSION  TECHNIQUE. 


5-1  Introduction 

Conventional  techniques  for  jet  noise  suppression  basically  rely  on  sub¬ 
division  of  the  main  nozzle  flow  into  a  number  of  separate  streams  in  such  a 
manner  as  to  produce  favorable  jet  mixing  through  interaction  of  the  elemental 
jets.  The  velocity  profile  theory  indicates  that  the  sound  power  generated  by 
a  free  jet  depends  only  on  the  mean-flow  velocity  distribution  in  the  mixing 
region  of  the  wake.  If  indeed  the  noise  generated  at  a  given  axial  location 
along  the  jet  is  primarily  dependent  on  the  velocity  profile  at  that  position 
and  sufficiently  independent  of  past  history  or  the  flow  and  the  means  employed 
to  achieve  the  velocity  distribution,  then  a  mixer -type  suppressor  could  effec¬ 
tively  he  replaced  by  another  scheme  which  would  alter  the  mean-flow  velocity 
profiles  without  first  subdividing  the  main  jet.  By  injecting  water  sprays 
Into  a  heated  jet,  the  temperature  and  velocity  of  the  flew  can  be  lowered  i9  ** 
selected  regions  of  the  wake,  altering  the  velocity  profiles.  While  it  is 
known  that  spraying  large  quantities  of  water  into  the  jet  wake  will  reduce  the 
sound  power,  generated  (Ref.  l)  the  quantity  of  water  required  is  such  that 
this  means  of  suppression  cannot  be  considered  as  practical  for  in-flight  use. 
However,  by  proper  utilization  of  limited  water  flow,  it  is  possible  to  alter 
the  mean-flow  velocity  profiles  rather  than  accomplish  a  gross  velocity  reduc¬ 
tion  for  the  entire  jet.  It  is  the  suppression  effect  of  this  means  of  water 
injection  which  has  been  investigated  and  is  reported  here  in  fulfillment  of 
the  proposal  work  statement  (Ref.  2). 
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The  objective  of  the  investigation  was  to  determine  feasibility  of  the 
selective  water  injection  concept  for  jet  noise  suppression.  As  previously 
stated,  the  noise  reduction  achieved  by  non-uniform  introduction  of  vat-ei 
sprays  should  be  greater  than  reduction  by  thorough  mixing  and  evaporation  of 
the  water  in  the  jet  wake,  if  the  scheme  is  to  be  successful.  Thus,  as  a 
standard  for  comparison,  a  theoretical  calculation  of  overall  power  level 
reduction  was  made  for  a  jet  with  completely  mixed  and  evaporated  water.  The 
results  are  shown  in  Figure  5-l»  and  are  based  on  the  assumption  that  acoustic 
power  depends  on  the  eighth  power  of  the  flow  velocity.  The  change  in  flow 
velocity  was  determined  from  heat  balance  calculation  for  the  air  and  water 
mixture. 

5-2  Experimental  Procedure 

The  experimental  approach  adopted  was  to  measure  noise  reduction  for  a 
heated  jet  with  water  injection  as  compared  to  an  unsuppressed  jet.  Rather 
than  make  a  complete  survey  of  the  sound  field  for  sound  power  determination, 
it  was  decided  to  measure  the  more  significant  acoustic  parameter  of  sound 
pressure  level  reduction  at  the  angle  of  maximum  jet  noise  radiation  (40*  from 
the  jet  axis).  On  the  basis  of  noise  directivity  studies  reported  in  References 
3  and  4,  it  is  seen  that  the  one  point  sound  pressure  spectrum  measurement  at 
the  angle  of  maximum  noise  can  be  used  to  estimate  the  sound  power  level 
spectrum. 

The  noise  measurement  system  consisted  of  a  Bruel  and  Kjaer  type  4133 
microphone  (with  flat  frequency  response  to  40  keps)  located  15  feet  from  the 
jet  nozzle  at  an  angle  of  40*  from  the  jet  axis.  Microphone  output  was  frequency 


-  228 


analyzed  by  means  of  a  Bruel  and  Kjaer  spectrometer,  and  recorded  on  a  B  4  K 
level  recorder. 

The  test  facility  used  was  an  outdoor  free  field  test  stand,  with  no 
significant  reflecting  surfaces  except  the  ground,  and  capable  of  continuous 
heated  flow.  The  test  nozzle  employed  was  a  4  inch  diameter  convergent  design, 
and  is  shown  with  water  spray  tube  configurations  in  Figure  5*2.  Flow  pressure 
and  temperature  were  measured,  as  well  as  water  and  air  mass  flow.  The  airflow 
measurements  were  used  to  monitor  the  effective  nozzle  area,  thus  providing  a 
rough  indication  of  water  evaporation  upstream  of  the  nozzle  exit.  A  number  of 
spray  tube  configurations  were  employed,  with  water  inlet  locations  at  two 
diameters  upstream  of  the  nozzle  exit  for  initial-  tests  and  twelve  diameters 
upstream  for  subsequent  runs.  The  tubes  were  immersed  radially  into  the  air- 
stream,  and  wete  equally  spaced  around  the  circumference  of  the  air  duct. 

Eight  tubes,  with  an  orifice  0.050  inches  dia.  in  the  end  of  each,  were 
used  in  the  configuration  with  tube  location  2  diameters  upstream.  Immersion 
depths  of  0  (tube  ends  flush  with  nozzle  wall),  2  inches,  and  alternate  tubes 
at  0  and  1  inch  were  tested.  Waier  flow  was  varied  from  0  to  3-5  gallons/minute 
for  airflow  conditions  of:  (a)  nozzle  pressure  ratio  «  1.8,  total  temperature 
-  l400aF,  and  (b)  pressure  ratio  *  2.7,  total  temperature  *•  1240°F.  Sound 
pressure  level  reduction  at  the  position  15  feet  from  the  nozzle  and  at  an 
angle  of  40*  from  the  jet  axis  was  determined  by  comparing  runs  with  and  with¬ 
out  water  flow  for  a  fixed  airflow  condition  and  spray  tube  geometry. 
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For  the  tests  with  water  injection  at  12  nozzle  diameters  upstream  of  the 
exit,  6  tubes  were  used.  Water  was  introduced  through  a  single  l/8"  diameter 
hole  in  the  end  of  each  tube,  and  in  an  alternate  scheme  water  was  sprayed 
through  8  holes  0.045"  dia.  located  along  the  length  of  each  tube.  Water  flow 
rates  up  to  27$  of  the  airflow  (by  weight)  were  achieved  for  flow  temperatures 
up  to  1340*F  and  over  a  range  of  pressure  ratios  from  1.35  to  2.15. 

5-3  Discussion  of  Results 

Initial  investigations  of  noise  reduction  by  means  of  selective  water 
injection,  utilizing  the  spray  tubes  located  2  diameters  upstream  of  the  nozzle 
exit,  indicated  that  little  evaporation  of  water  occurred  in  the  flow  region 
near  the  nozzle  exit.  Even  at  the  maximum  water  iluw  rate  (13$  of  the  airflow) 
only  4  db  sound  pressure  level  reduction  was  achieved,  as  shown  in  Figure  5.3. 

Based  on  the  theoretical  curves  in  Figure  5*1»  for  complete  mixing  of  the 
water  with  the  heated  airstreaa,  a  6  db  sound  power  level  reduction  would  be 
expected.  While  it  appears  that  this  amount  of  suppression  at  the  source  might 
be  obtainable  in  practice  by  complete  mixing  of  water  and  air,  the  configura¬ 
tion  tested  did  not  achieve  noise  reduction  as  intended  by  alteration  of  flow 
velocity  profiles.  No  change  in  airflow  occurred  as  water  flow  was  increased, 
indicating  that  only  a  limited  amount  of  evaporation  could  have  taken  place  in 
the  nozzle,  since  formation  of  steam  would  effectively  decrease  the  flow  area 
available  for  the  heated  air. 

Examination  of  the  sound  pressure  reduction  spectra  (Figure  5.3)  reveals 
that  at  the  highest  pressure  ratio  the  noise  reduction  was  shifted  to  a  lower 
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frequency.  This  is  explained  by  consideration  of  the  increased  jet  velocity, 
since  the  water  droplets  were  transported  further  downstream  (toward  the  loca¬ 
tion  of  lower  frequency  sound  source)  before  evaporation  and  resultant  cooling 
and  deceleration  of  the  flow. 

In  order  to  achieve  more  complete  vaporization  of  the  water  injected  into 
the  flow,  further  tests  were  conducted  using  the  spray  tubes  located  12  uozzle 
diameters  upstream  of  the  exit-  Comparison  of  suppression  for  the  two  spray 
tube  designs  (one  hole  in  the  end  of  each  tube  compared  to  8  smaller  holes 
spaced  along  the  tube)  showed  no  significant  difference  between  the  coarse  and 
fine  sprays.  Bather,  noise  reduction  was  found  to  depend  essentially  on  water/ 
air  ratio  for  a  given  airstream  pressure  and  temperature.  Overall  sound  pres¬ 
sure  level  reduction  as  a  function  of  vater/alr  ratio  is  shewn  in  Figure  5-4, 
and  it  is  apparent  that  suppression  is  a  linear  function  of  water/air  ratio 
over  the  range  tested.  Comparing  the  curves  for  equal  pressure  ratios  hut 
varying  temperatures  it  is  evident  that  noise  reducti  increases  as  the  flow 
temperature  rises,  and  a  maximum  reduction  is  achieved,  after  which  the  addi¬ 
tional  temperature  rise  adversely  affects  suppression.  Interpolated  data  from 
Figure  5-4  for  water/air  ratio  of  15$  is  replotted  in  Figure  5*5  bo  indicate 
the  effects  of  water  residence  time  in  the  airstream  (flow  velocity  effect). 

It  is  the  combined  effects  of  increase  in  flow  temperature  causing:  l)  more 
rapid  vaporization  of  waxer,  2)  decreasing  droplet  residence  time,  and  3)  pro¬ 
viding  less  noise  reduction  capability  (as  shown  in  Figure  5-l)  which  cause  the 
maximization  of  suppression  noted  in  Figure  5-4. 
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The  calculation  adverse  effect  of  temperature  results  from  the  fact  that 
for  a  given  water  air  ratio  the  flow  temperature  change  accompanying  complete 
mixing  and  evaporation  is  essentially  independent  of  initial  temperature  and 
thus,  percentage-wise,  the  coding  effect  decreases  with  increase  in  Initial 
temperature. 

The  effects  of  water  injection  on  the  jet  noise  sound  pressure  level 
spectrum  are  shown  in  Figure  5.6.  Seduction  is  achieved  over  a  wide  range  of 
frequencies,  and  it  is  significant  that  even  at  low  pressure  ratio  suppression 
is  accomplished  in  the  high  frequency  region  (for  the  test  results  shown,  the 
sound  pressure  level  spectrum  peak  occurs  at  1250  cps)-. 

5-4  Conclusions 

Results  of  the  selective  water  injection  investigation  indicate  that  the 
desired  alteration  of  jet  mixing  cannot  he  achieved  by  this  means  in  order  to 
affect  noise  reduction.  Although  the  mean-flow  velocity  profiles  can  be  altered 
tl irou&h  use  of  selective  sprays,  the  jet  mixing  is  not  greatly  influenced  and 
the  suppression  accomplished  is  less  than  would  be  possible  through  uniform 
mixing  of  water  and  air. 

The  adverse  effect  of  high  temperature  and  flow  velocity  on  noise  reduc¬ 
tion  capability  of  water  injection  and  the  large  watery air  ratios  required 
indicate  that  this  scheme  might  not  be  practical  for  application  to  existing 
aircraft  and  engines. 

Effects  on  nozzle  aerodynamic  performance  have  not  been  measured,  but 
losses  would  be  small  if  water  evaporation  and  resultant  cooling  of  the  flow 
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occurs  primarily  external  to  the  nozzle,  as  is  indicated  by  test  results. 

Jet  noise  suppression  by  means  of  water  injection  is  unusual  in  that  high 
frequency  noise  reduction  can  be  achieved,  and  furthermore,  suppression  can  be 
accomplished  for  low  pressure  ratio  flow. 
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Figure  5.1  -  Calculated  Reduction  of  Sound  Power  Level  for  C&ispleta 
Evaporation  and  Mixing  of  Water  viih  Jet  Exhaust. 
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Figure  $.2  -  Schematic  of  Water  Injection  Kottle 
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7irur«  5.3  -  Sound  Pressure  Level  Reduction  of  Water  Injection  Koacle  Measured  at  Uo  from  Jet  Axis 


Injection  Nozzle  Overall  Sound  Pressure  Level  Reduction 
;le  of  Maximum  Noise  (^0*  from  Jot  Axis) 
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Figure  5»5  -  Vater  Infection  Nozzle  Overall  S?L  Reduction 
at  Angle  of  Kaxiaiuai  Noise 
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Figure  5,6  -  Water  Injection  Nozzle  Noise  Reduction  at  Angle  of  Maximum  Noise 


